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In this new text the author presents an ex- 
position of the analytic geometry of three- 
dimensional space. The material covers the 
standard topics of space analytic geometry but 
provides a treatment of the subject which per- 
mits immediate generalization to n dimensions. 
This treatment ties the subject to modern math- 
ematics, and, in particular, to modern algebra. 
The use of the theory of vector spaces and ma- 
trices permits a major simplification in the 
proofs and in the exposition in general. Thus 
the aim of the book is to provide a modern and 
simpler treatment of the subject matter which 
permits easy generalization and fits the sub- 
ject into its proper place in modern mathe- 
matics. 

The early part of the text is simplified by the 
use of the concepts of inner and scalar prod- 
uct. A brief but adequate chapter on the theory 
of matrices provides a full, clear exposition 
of the principal axis transformation in the 
n-dimensional case. An additional feature is 
provided by the chapter on spherical coordi- 
nates where the mathematics of the approxima- 
tions used in actual physical measurements of ' 
direction and distance in rotated coordinate 
systems is presented. This material is r</^ usu- 
ally found in texts on solid analytic geometry. 

Chapters I and II contain a treatment of the 
equations of lines and planes. After a prelimi- 
nary study of the linear operations in n-dimen- 
sional vectors and inner products are interpreted 
geometrically, and from a consideration of 
scalar products and axis translations, the para- 
metric equations of a line are obtained. The 
vector approach then yields a very simple 
derivation of the normal form of an equation 
of a plane, and the standard forms of plane and 
line equations are rather immediate conse- 
quences. 

Chapter III presents classical elementary sur- 
face and curve theory. Chapter IV contains the 
usual treatment of spheres, and Chapter V the 
classical descriptions of quadric surfaces in 
standard position. This chapter ends with a 
rather novel classification of quadrics accord- 
ing to certain invariants. 

Chapter VI is an exposition of that part of the 
theory of matrices needed for a complete de- 
velopment of the so-called principal axis trans- 
formation. A full account of the orthogonal 
reduction of a real quadratic form in n varia- 
bles is given, and the theory is applied in 
Chapter VII to the three-dimensional case of 
quadric surfaces. A discussion of the sym- 
metries of "dric surfaces is included. 

Chapter Vi spherical coordinates, contains 

a discussion of some practical aspects of the 
theory of rotations and translations of axes in 
space. 

The final chapter offers a brief presentation of 
the elements of projectiye geometry. The ef- 
fect on the theory of linear transformations 
of the use of homogeneous coordinates is given 
.and the chapter contains a rigorous uidtrix 
proof of the invariance of the cross ratio under 
protective transformations. 
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PREFACE 

In a recent text on college algebra the author gave a brief pre- 
senlation of what seems to him to be the best basis for a modern 
course on plane analytic geometry, that is, the algebraic vector 
approach. The present text contains an extension of this 
approach, yielding an exposition in full for the three-dimensional 
case, and thereby ties up the study of space analytic geometry 
with the theory of vectors and matrices. 

Chapters 1 and 2 contain a treatment of the equations of lines 
and planes. After a preliminary study of the linear operations 
on n-dimensional vectors, rectangular coordinates are introduced, 
three-dimensional vectors and inner products are interpreted 
geometrically, and, from a consideration of scalar products and 
axis translations, the parametric equations of a line are obtained. 
The vector approach then yields a very simple derivation of the 
normal form of an equation of a plane, and the standard forms of 
plane and line equations are rather immediate consequences. 

Chapter 3 contains an exposition of classical elementary sur- 
face and curve theory. Chapter 4 contains the usual treatment 
of spheres, and Chapter 5 gives the classical descriptions of 
quadric surfaces in standard position. The latter chapter ends 
with a rather novel classification of quadrics according to certain 
invariants. 

Chapter 6 is an exposition of that part of the theory of matrices 
which is needed for a complete development of the so-called 
principal axis transformation. A full account of the orthogonal 
reduction of a real quadratic form in n variables is given, and the 
theory is applied in Chapter 7 to the three-dimensional case of 
quadric surfaces. This latter chapter ends with a discussion of 
the symmetries of quadric surfaces. 

The first seven chapters of this text provide an exposition of 
the basic topics of solid analytic geometry, the material being 
adequate for a one-quarter course on the subject. The remain- 
ing chapters contain additional material for longer courses or out- 
side reading. Chapter 8, on spherical coordinates, contains a dis- 
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cussion of some practical aspects of the theory of rotations and 
translations of axes in space. It is quite clear that rectangular 
coordinates are not as practical for actual measurements as are 
the coordinates of range, azimuth, and elevation, and in this 
chapter methods are developed for actual computation of the 
changes in these measurable coordinates after translations or 
rotations of axes. The chapter ends with a discussion of gno- 
monic charts. 

Our final chapter contains a brief presentation of the elements 
of protective geometry. The effect on the theory of linear trans- 
formations of the use of homogeneous coordinates is given, and 
the chapter contains a rigorous matrix proof of the invariance of 
the cross ratio under projective transformations. It is hoped that 
the use of modern algebraic techniques in this and in the earlier 
chapters of the present text will serve to make the subject of 
solid analytic geometry fit better in the teaching of modern 
mathematics than it has in the past. 

ADRIAN ALBERT 

CHICAGO, ILL. 
March, 1949 
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CHAPTER 1 
COORDINATES AND LINES 

1. Vectors. A sequence P = (xi, . . . , x n ) of n numbers Xi is 
called an n-dimensional vector. The elements x\, . . . , x n are 
called the coordinates of P and # is the ith coordinate. We shall 
limit our attention to real vectors, i.e., to vectors whose coordinates 
are real numbers. 

The vector whose coordinates are all zero is called the zero 
vector and will be designated by 0. A real vector may be inter- 
preted as a representation, relative to a fixed coordinate system 

with O as origin, of a point in ?i-dimensional real Euclidean space. 


It may also be interpreted as the line segment OP directed from 

O to P. These interpretations have little intuitive significance 
except for the cases n ^ 3, and we shall carry out the details in 
this text for the case n = 3. 

The sum P + Q of two vectors P = (xi, . . . , x n ) and Q 
= (yi, . . , 2/n) is the vector (xi + yi, . . . , x n + y n ) whose 
ith coordinate is the sum X{ + ?/ t of the ith coordinate of P and 
the ith coordinate of Q. We leave the verification of the follow- 
ing simple results to the reader: 

Lemma 1. Addition of vectors is commutative, that is, P + Q 
= Q + P for all vectors P and Q. 

Lemma 2. Addition of vectors is associative, that is, (P + Q) 
+ R = p -f (Q + R) for all vectors P, Q, R. 

Lemma 3. The zero vector has the property that P + = P 
for all vectors P. 

Lemma 4. Let P = (xi, . . . , x n ). Then the vector P 
= ( xi, . . . , Xn) has the property that P + ( P) =0. 

Lemma 5. // P and Q are any vectors the equation P + X = Q 
has the solution X = Q + (~P). We call this vector the difference 
of Q and P and write X = Q P. Then the ith coordinate of 
Q P is the difference of the ith coordinate of Q and the ith coordinate 
off. 

1 
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EXERCISE 

Verify the five lemmas. 

2. Scalar multiplication. If a is a number and P = (xi, 
. . . , x n ) is a vector, we define the scalar product of P by a to be 

aP = 



Evidently IP = P, (-l)P = -P, OP = 0. The reader should 
verify that 

a(bP) = (ab)P, (a + 6)P = aP + bP, a(P + Q) 

= aP + aQ 

for all scalars a and b and all vectors P and Q. 
A sum 

P = 



of scalar products a/P/ of vectors P/ by scalars a ; , is called a linear 
combination of Pi, . . . , P m . We shall say that PI, . . . , P m 
are linearly independent vectors if it is true that a linear combina- 
tion aiPi + ' ' + flwPm if and only if ai, . . . , a m are all 
zero. If PI, . . . , Pm are not linearly independent, we shall 
say that Pi, . . . , P m are linearly dependent. 

Let Ei be the vector whose ith coordinate is 1 and whose other 
coordinates are all zero. Then 

P = (Xlj . . . , X n ) = X i E l + + XnE n . 

Thus every vector is a linear combination of EI, . . . , E n . If 

XiEi + ' * * + XnE n = P = 0, then (xi, . . . , x n ) = 0, that is, 

Xi = #2 = = #n = 0. It follows that Ei y . . . , E n are 
linearly independent. 

EXERCISES 

1. Show that if P = (#1, . . . , x n ) and Q = (?/i, ...,?/) are not 
zero then P and Q are linearly dependent if and only if Q is a scalar 
multiple of P. 

2. Show that if PI, . . . , P m are linearly independent and P w +i is 
another vector then Pi, . . . , P w , P m +i are linearly dependent if and 
only if Pm+i is a linear combination of Pi, . . . , P m . 

3. Compute the following linear combinations of Pi = (1, 1, 2, 3), 
P 2 = (0,l, -1,2),P 3 = (-2,1, -1,2). 

(a) 2P, + P, + P, (b) P l + 3P, - 2P 3 (c) 3P, + 2P 2 - 4P 3 
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4. Use Exercises 1 and 2 in determining which of the following sets of 
three vectors are linearly independent sets. 

(a) (1,-1,2),(1 ( 1,0),(0, -1,1) 

(6) (2,1,1), (1,- 1,1), (5, 4, 2) 

(c) (1,0, -2), (2, -1,2), (4, -3, 10) 

(d) (1, -1, !),(-!, 2, !),(-!, 2, 2) 

(e) (1,0, -1, 1), (0, -1, 1, -1), (4, -1, -3,4) 
(/) (5, 1, -2, -0), (1, 1, 0, -2), (2, -1, -1, 0) 
(17) (1,0, 0,0), (1,1, 1,1), (3, 1,1,1) 

(h) (1, 1, -1, 2), (2, 2, -2, 3), (3, 3, -2, 6) 

5. Prove that any three two-dimensional vectors are linearly depend- 
ent. 

6. Prove that any four three-dimensional vectors are linearly de- 
pendent. 

3. Inner products. If P = (xi, . . . , x n ) and Q = (7/1, . . . , 
t/ n ) are any two vectors, we shall call the number 

(1) P ' Q = *i?/i + + x,#n 

the inner product of P and Q. Evidently, P Q = Q P. 
The norm of a vector P is defined to be the inner product 

(2) P P = Xi 2 + ' ' ' + X n 2 . 

If P is any real vector, the number P P ^ and has a nonnega- 
tive square root 



(3) t = 

which we shall call the length of P. 

A vector P is called a unit vector if P P = 1. Thus a real 
unit vector is a vector whose length (and whose norm) is 1 . 

Lemma 6. Every real nonzero vector is a scalar multiple of 
exactly two unit vectors. These are the vectors U = t"" J P and U, 
where t is the length of P. Then if P = tU, where t > and U is 
a unit vector, the number t is the length of P. 

For proof we first let P = tU where U = (ui, . . . , u n ) is 
a unit vector. Then P -P = (toi) 2 + + (tu n )* = 



+ . . + u f) = t\ and = V?" 7 * 5 ; * = V^ 7 ?, if * ^ 0. 
Conversely, let f/ = fT, where t = \/P ' P. Then 17 17 
= (r^i) 2 + + (r 1 ^) 2 - t~*(xS + + x n 2 ) = 1 and 
U is a unit vector. The vector U= t~ l P is clearly also a 
unit vector. 
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EXERCISES 

1. Give the norms and lengths of the following vectors: 

(a) (2,2, -1) (d) (1, -1,1, -1) 

(6) (1,1,0) (e) (1, -1,2, 1) 

(c) (1, -4,8) (/) (3,2, -1,1,1) 

2. Give the unit vector (\/P P)~~ 1 P for each vector of Exercise 1. 

4. The angle between two vectors. If P = (xi, . . . , x n ) 
and Q = (t/i, . . . , t/ n ) are any two real nonzero vectors, the 
difference 

(4) '(P-PXQ-Q) - (P-Q) 2 ^0. 

For (P P)(Q - Q) = (zi 2 + + x n *)(yi* + + yS} is 
the sum of Xi 2 ?/i 2 + x 2 2 2/2 2 + + x n 2 y n 2 and all expressions 
of the form (xT/y) 2 + (#/2/;) 2 for 1 ^ i < j ^ n. The square 
(P - Q) 2 - (zij/i + +x n 7/ n ) 2 is the sum of x^yS + xfy? 
+ ' + # n 2 y n 2 and all products of the form 2o: t ?/ t x ; 7// for 
1 ^ i < j ' ~ ^- The difference then is the sum of all expressions 
of the form (>,?/,) 2 + (z/j/,-) 2 - ^x^XfUi = (xy/ x/y t ) 2 for 1 <; i 
< j ^ HJ and must be nonnegative. 

The numbers P P, Q Q, and (P Q) 2 are all positive, and we 
have shown that 



It follows that there exists an angle 6 between and 180 such 
that 

p . o 

(5) cos 6 = ^L^-^^. 

^ ' 



We define this angle 6 to be the angle between the vectors P and Q. 
Two vectors are said to be orthogonal (i.e., perpendicular) if 
cos = 0. Then P and Q are orthogonal if and only if their inner 
product 

(6) P Q = Zl^ + + XnVn = 0. 



Thus, if P and Q are any vectors, we multiply corresponding 
coordinates and add the products. The sum so obtained is zero 
if and only if P and Q are orthogonal. 
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EXERCISES 

1. Compute P Q for each of the following vector pairs P, Q : 

(a) (1, 1, -1), (1, 0, 1) (/) (1, -1, 1, 1), (1, 1, 1, 0) 

(6) (1, 2, 3), (-1, 1, -1) (?) (2, 3, -1, 6), (3, -2, 6, 1) 

(c) (1, 1,2), (0, -1,1) (K) (1,2,3,4), (2,- 1, - 1,1) 

(d) (-1,0,1), (2, 1,1) .0) (4, -6,1,2), (1,2, -1,2) 

(e) (-1, 3, 2), (1, 1, -1) (j) (3, 1, -1, 1), (0, 1, 1, 0) 

2. Which pairs are orthogonal? 

3. Compute cos 6 for each nonorthogonal pair. 

5. Directed lines. Directed lines are frequently used in the 
geometry of three-dimensional Euclidean space, i.e., in ordinary 
solid analytic geometry. Every pair of distinct points P and Q 
in space determines a line passing through P and Q. We shall 
use the notation PQ for this line and shall prefix the word ray 
when we mean the ray PQ, which is the half line from P through Q. 

Let us assume that a unit of measurement has been prescribed 
and that we have measured the length of the line segment joining 
P to Q in terms of this unit. The result is a real number that is 
positive if P and Q are different points and is zero only when P 
and Q coincide. We shall use the notation |PQ| for this measure- 
ment of length. 

When P and Q are points on a directed line, we shall use the 
symbol PQ for the signed length of the segment joining P to Q 
and directed from P toward Q. Then PQ = |PQ| if the direction 
from P to Q is the positive direction on the line, and PQ = |PQ| 
if the direction from P to Q is opposite to the positive direction on 
the line. See Fig. 1 in which PQ > and RQ < 0, and note that 
in all cases PQ = -QP. 

1 t 1 

P Q R 

FIG. 1. 

1 1 1 

P R Q 

FIG. 2. 

If P, Q, R are on a directed line, it should be clear from Figs. 1 
and 2 that PQ + QR = PR. This equation may be generalized 
to the case of any finite number of points on a directed line and 
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the generalized equation is 



6. Orthogonal projections. A theorem of solid geometry 
states that through a given point P there exists precisely one 
plane perpendicular to a given line L. This plane intersects L 
in a point P' such that the line PP f is perpendicular to L. We 
shall call P' the orthogonal projection of P on L. 

If P and Q are any two points, we shall designate by PQ the 
line segment which joins P to Q and which is directed from P to Q. 
Project P and Q orthogonally on a directed line L, and obtain 
projections P' and Q 1 '. Then we define the orthogonal projection 

of PQ on L to be the signed length P'Q' . It follows that the 
orthogonal projection of QP is the negative Q'P' of the orthogonal 

projection of PQ. 

A directed broken line joining two points P and Q is the 
geometric configuration consisting of the directed line seg- 

ments PPi, PiP2, , PnQ for any finite number n of points 

Pi, ... , P n . Let us use the notation PPi P n Q for 
such a configuration and define the orthogonal projection of 



PPi PnQ to be the sum P'Pi' + PJP*' + . + 

By formula (7) this sum is equal to P'Q' and we have proved the 

following: 

Theorem 1. The orthogonal projection of any directed line 

-4 

segment PQ on a directed line L is equal to the orthogonal projection 
on L of any directed broken line from P to Q. 

As we have said, a ray PQ is a half line that begins at the point 
P, passes through Q, and extends indefinitely. If PQ and PR 
are two rays from the same point P, there is a unique angle 6 





Q 
FIG. 3. 
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between them such that 0^0^ 180. We define this angle to 
be the angle between the ray PQ and the ray PR. 

Let P and Q be any points and P' and Q' be their respective 
projections on a directed line L, as in Fig. 4. Then Q, P, P' are 
three of the vertices of a parallelogram, and we find the fourth 
vertex R by drawing a line segment P'R such that \P'R\ = |PQ|. 




FIG. 4. 

Define the angle between the ray PQ and the line L to be the 
angle between the ray P'R and the ray from P' in the positive 
direction on L. By the standard ratio definition of the cosine 
of an angle we have 

cos = i^lf = p| 

We have proved the following : 

Theorem 2. Let 6 be the angle between PQ and a directed line L. 

Then the projection of PQ on L is equal to |PQ| cos 0. 

7. Rectangular coordinates in ordinary space. A rectangular 
Cartesian coordinate system in ordinary three-dimensional real 
Euclidean space is a certain one-to-one correspondence between 
the points of space and three-dimensional real vectors (x, y, z). 
The construction of the correspondence begins with the construe- 
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tion of three mutually perpendicular directed lines intersecting 
at a point called the origin of the coordinate system (see Fig. 5). 
The three lines are called coordinate axes. The first of them 
is a vertical line directed upward. It is called the z axis. The 
second line is a horizontal line in the plane of the book and is 
directed to the right. It is called the y axis. The remaining 
line is the x axis. It should be thought of as a line perpendicular 
to the plane of the book and directed toward the reader. The 



R 



Q 



FIG. 5. 

specification of a coordinate system will be completed as soon as 
a unit of measurement, which will be used for all measurements 
of lengths of lines, is given. This is usually done by specifying 
a unit point U on the x axis such that OU is the unit of length. 

Let P be any point in space, x be the projection of OP on the 
x axis, y be the projection of OP on the y axis, and z be the pro- 
jection of OP on the z axis. Then the vector (x, y, z) is uniquely 
determined by P and we write P = (x, t/, z). 

Conversely, if (x, y, z) is given, we can draw a plane perpendicu- 
lar to the x axis and through a point C on the x axis such that 
OC = x. All points P on this plane have the property that the 

projection of OP on the x axis is x. We may construct another 
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plane perpendicular to the y axis through a point B on the y axis 
such that OB = T/, and finally a plane perpendicular to the z axis 
and through a point R on the z axis such that OR = z. The 
intersection of these three mutually perpendicular planes is the 
unique point which is such that P = (#, ?/, 2). 

This completes our description of a rectangular coordinate 
system in ordinary space. Let us now observe some simple 
properties of a coordinate system. It should be clear that each 
pair of axes defines a plane. We call the three mutually perpen- 
dicular planes that are so defined the coordinate planes. The 
x, y plane is that determined by the x and y axes and z is the 
perpendicular distance from the x y y plane to P = (x, y, z). 
The ?/, z and z, x planes are defined similarly, and x and y are also 
perpendicular distances from coordinate planes to P. 

The coordinates x, y, z are any real numbers and therefore may 
be positive, negative, or zero. The reader should verify the 
following statements and answer the questions: 

1. The x axis is the set of all points P such that y = z = 0. 
What are the corresponding equations for the y axis? The z axis? 

2. The x, y plane is the set of all points such that z = 0. 
What is the corresponding equation for the y, z plane? The z, x 
plane? A plane parallel to the x, y plane and three units above 
it? A plane parallel to the y, z plane and three units behind it? 

3. The coordinate planes divide all of space into octants. 
The forward, right-hand, upper octant is that where x ^ 0, 
y ^ 0, z ^ 0. We call this the first octant, label it 7, and draw- 
most of our diagrams as if the points being studied are all in the 
first octant. This will restrict only the generality of the diagrams 
but not that of the mathematical arguments. 

4. The octants labeled II to VIII are those where the corre- 
sponding sets of signs are given as follows :(, + , +), ( , --,4-), 

(+, -, +), (+, +, -), (-, +, -), (-, ~, -), (+, -, -)- 
Describe the octants in terms of the words forward and backward, 
right and left, and upper and lower. 

EXERCISES 

1. Verify the statements, and answer the questions listed above. 

2. How can an equation be used to describe the set of all points with 
equal x and y coordinates? Give the equation and the nature of the 
geometric configuration. 
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3. Let a given point P be such that the projection OP' of OP on the x 

axis is equal to cos a, where a is the angle between the positive x axis 

* 

and OP. What is the projection of OP' on OP? 

8. The length of a vector in space. A directed line segment 

OP is usually called a vector and this is the inspiration of the 
name we have given to number sequences (x, y, z) . We shall 
now show that if P = (x, y y z), then the length of OP is given by 



\OP\ = V? + y 2 + z\ 

that is, our formula for the length of a vector is an actual length 
in the three-dimensional case. 

Use Fig. 5 and observe that \OA\ is the diagonal of a rectangle 
whose sides are x and y. Then \OA\ 2 = x 2 + y 2 . Also, AP = z 
and \OP\ is the hypotenuse of a right triangle whose legs arc \OA\ 
and \AP\. By the theorem of Pythagoras |OP| 2 = |OA| 2 + z 2 
= x 2 + y 2 + z 2 . This yields the formula above. 

ORAL EXERCISE 

Give the lengths of the following vectors: 

(a) (-1,2,3) (d) (-2,2,1) (0) (0,0, -3) 

(6) (1,1,2) (e) (3,6,0) (A) (3,0, -4) 

(c) (1, -2, -2) (/) (-1,1,1) (i) (5,7, -1) 

9. Lines through the origin. Every line OP through the origin 
is composed of two rays. One of them is the ray OP and the 
other is the ray OQ, where Q = P. For if P (x, y, z} is on 
a line and P = ( x, y, z), it should be evident geometri- 
cally that P is also on the line. 

Every ray OP contains a point U such that \OU\ = 1. Then 
U = (X, ju> v) is a unit vector and 

(8) X 2 + M 2 + * = 1. 

There is exactly one other unit vector on the line OP and this is 
the point U = ( X, ju, ?) 

Define a to be the angle between the positive x axis and the 
ray OP, to be the angle between the positive y axis and the ray 
OP, and 7 to be the angle between the positive z axis and the 
ray OP. By Theorem 2 and our definition of coordinates 

(9) x = \OP\ cos a, y = |OP| cos 0, 2 = |OP| cos 7. 
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We call a, /?, 7 the direction angles of the ray OP, and it should 
be clear that any two of them and the sign of the third uniquely 
determine this ray. In fact, the set of all rays making an angle 
with the positive x axis form a cone of lines, the set of all rays 
making an angle /? with the positive y axis form a cone of lines, 
and the ray OP is one of the two intersections of these two cones. 
The dependence of 7 on a and /3 is expressed by formula (11) 
below. 

Since the point U is on the ray OP and |OC7| = 1, we may apply 
formula (9) to see that 

(10) X = cos a, JJL ~ cos |8, v cos 7. 
By formula (8) 

(11) cos 2 a + cos 2 )8 + cos 2 7 = 1. 

This relation expresses the dependence of the three angles a, /J, 7. 
The numbers cos a, cos /3, and cos 7 are called the direction 
cosines of the ray OP, and we have seen that they satisfy formula 

(11) and are the coordinates of the unique unit vector on this ray. 
If Q = (xij 7/1, zi) is on the ray OP, then x\ = \OQ\ cos a, 

yi = \QQ\ cos j8, 21 = \OQ\ cos 7. By formula (9) we have 

(12) Q = fP, 



* |OP| 

If Q is on the line OP but not on the ray OP, then Q is on the 
ray OP and Q = tP. We have proved the following geometric 
interpretation of our algebraically defined operation of scalar 
product : 

Theorem 3. A point Q is on a line through the origin and 
another point P = (x, y, z) if and only if Q is a scalar multiple 
t(x, y, z) of P. Then |t| = IOQHOPI- 1 and t art <Q accord- 
ing as Q is or is not on the ray from through P. 

EXERCISES 

1. Find the direction cosines of the ray OP for each of the following 
vectors P: 

(a) (-1, -1,1) (d) (1,0, -1) 

(6) (1,2,2) (e) (-1,1,2) 

(c) (-2,1,1) (/) (2, -1,3) 



12 
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2. What previous exercise is identical with Exercise 1 except for the 
actual numbers involved? 

3. What are the direction cosines of the coordinate axes? 

10. The angle between two vectors in space. Let P and Q be 

any two points distinct from the origin and 9 be the angle between 
the ray OP and the ray OQ. Let U = (A, M> v) be the unit vector 




B 



X 

FIG. 6. 

on the ray OP and V = (Ao, MO, ^o) be the unit vector on the ray 

OQ. Project OU orthogonally on the line OQ directed from O 
toward Q. The projection is \OU\ cos 6 = cos 6. But this is the 

same as the projection on OQ of the broken line OABU of Fig. 6. 

The projection of AB is the same as the projection of OC and, 
since the cosine of the angle between the ray OC and the ray OQ 

is Ao, this projection is AA . Similarly, the projection of OA on 


OQ is MMO and that of BU on OQ is WQ. Then 

(13) cos 6 AAo + MMO + VVQ. 

If P = (x, y y z) and Q = (x , ?/o, Zo), then formula (13) becomes 
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/ i 4 \ / XX *T" 2/2/o i zzo 

(14) cos e = .^^==^== ya 



For |OP| = Vz 2 + y 2 + z\ \OQ\ = Vzo 2 + 2/o 2 + zo 2 , and * 
= |OP|X, y = |OP| M , 2 = |OPh * = |OQ|X , 2/0 = |OQ| M o, 20 
= |OQ|z/ . We also see that a line through and P = (x, 2/, z) 
and a line through and Q = (z , 2/o, Zo) are perpendicular if and 
only if xx Q + yz Q + zz = 0. 

The exercises of Sec. 4 are exercises on the material of this 
section. 

ORAL EXERCISE 

Show that the triangles whose vertices are P, Q, and the origin are 
right triangles in the following cases: 

(a) P = (2, -1,4), Q = (3,2, -1) 

(6) P = (1, -1,0),Q = (1,1,6) 

(c) P - (2,1, -2), Q = (3,4,5) 

(d) P = (1,1, -2),Q = (1,1,1) 

11. Translation of axes. The correspondence between points 
P and vectors #, ?/, depends on the use of a fixed coordinate sys- 
tem. If the coordinate system is altered, so is the correspondence. 

It is desirable to investigate the effect of a translation of axes 
on this correspondence. Such a change is the result of setting 
up an x f , y', z' coordinate system in which the (new) x f axis is 
parallel to the x axis, the y' axis is parallel to the y axis, and the 
z' axis is parallel to the z axis. Thus the x', y f ', z 1 coordinate 
system may be conceived of as having been obtained by a motion 
called a translation in which axes are moved parallel to themselves 
and the main effect is a change in origin. 

Every point P in space will now have two sets of coordinates. 
We will call the first set the x, y, z coordinates and the second set 
the x' , y', z' coordinates. The origin 0' of the second coordinate 
system has a set of x, y, z coordinates that we shall designate by 
x , 2/0, 20, as in Fig. 7. 

Suppose now that P has coordinates x, y, z. This means that 
P is z units above the x, y plane. Note that z units above means 
z units below if z is negative. The x', y 1 plane is 2 units above 
the Xj y plane and it should be clear then that P is only z z 
units above the x' ', y' plane, that is, z' = z z . By similar 
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considerations we obtain the relations 

(15) x' = x - x 0) y' = y - yo, z' = z - z 

between the two sets of coordinates of P. 

Conversely, if x', y', z' and XQ, y, ZQ are given, we can use the 
relations 

(16) x = x f + xt>j y = y' + 2/0, z = z' + ZQ, 
between the two sets of coordinates. 

z 9 




FIG. 7. 

Formula (16) is sometimes used in making a translation of 
axes, which simplifies what we shall later call an equation of a 

surface. 

EXERCISES 

1. What translation of axes will simplify the following equations of 
surfaces? 

,. (^ - l) 2 , (y + 2) (z- 3)* _ 
W 4 -I e 9~~ ~ 

(6) 2(x + 2) 2 - (y - 3) 2 + (z + I) 2 = 5 

(c) x 2 - 2x + 3y 2 - 6?/ + 42 2 - 12 = 

(d) 2x 2 + 3y 2 - z 2 + Go; - 12z/ + z = 9 

(e) 3a; 2 - 3?/ 2 + Oz 2 = 6a; - 9y + 24z 
) 3x - 2/ 2 + s 2 = 6 
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2. Find the #', y', z' coordinates of each of the following points P if 
there is a translation of axes with origin at 0'. 
(a) P = (1, -1,1), 0'= (2,2, -1) 
(b)P = (0,0,0), 0'=(1, -1,1) 
( C )P = (-1, -1,1), 0'= (-1, -1,2) 
(d)P = (3,1,1),0' = (-2,2,1) 
(e)P = (l,2, -3),0' = (1,0, -1) 
(/) P = (2,4,0),0' = (0, -1,0) 

3. Let the vectors P in Exercise 2 be the a/, ?/', 2' coordinates of cor- 
responding points. Give their x, y, z coordinates. 

4. Find the x', y', z' equation of the surfaces given by the equa- 
tions of Exercise 1 after a translation of axes moving the origin to 
0' = (1, -1,2). 

12. Geometric addition of vectors. Let P\ (xi, y\, Zi) and 
P* = (#2, 2/2, 22) be two vectors. Then Pi#P 2 determines a plane 




FIG. 8. 

and, as in Fig. 8, we can determine a point Pa in this plane such 
that OP 3 is parallel to PiP 2 and in a corresponding direction, 
|OP 3 | = |PiP 2 |. The coordinates of P 2 , relative to an x' , y', z 1 
coordinate system with origin at Pi, will be (x 2 Xi, y<z y\, 
z 2 ~ Zi). But P 3 is located with respect to the x, y, z axes 
exactly where P 2 is located with respect to the x' ', y', z' axes. It 
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follows that 

(17) P 3 = P 2 - Pi = (xi - xi, ?/ 2 - j/!, z 2 - zi). 

This shows that to construct the vector difference P 3 = P 2 Pi 
we construct OP 3 parallel to PiP2 and having the same length. 

Since P 3 = P 2 Pi, we have P 2 = P 3 + PI. This shows that 
if Pi and P 3 are given, their sum, Pi + P 3 , is the end point of the 
diagonal of the plane parallelogram determined by the directed 
line segments OP\ and OP 3 . 

13. The length of a line segment. If Pi = (xi, y\, Zi) and 
P = (#2, 2/2, 2) are points in space, the line segment PiP 2 has the 
same length |PiP 2 | as 0P 3 , where P 3 = P 2 - PI (or P, = P, 
- P 2 ). It follows that 



(18) |PiP 2 | = 



EXERCISES 



1. Compute the lengths of the line segments joining tho following; 
pairs of points: 

(a) (0,0,0), (2, -2,1) (e) (1, 2, 3), (-1, -1, 1) 

(b) (1, -1, 1), (3, -3, 2) (/) (-1, 1, 1), (0, 1, 2) 

(c) (-1, -3,4), (-2,1, -4) (g) (1,2, 1), (-1,0, -2) 

(d) (1, -1, 1), (0, -1, 1) (h) (1, 2, -3), (3, -2, 1) 

2. Prove that the triangles formed by line segments joining the follow- 
ing triples of points are isosceles: 

(a) (-1, -3,4), (-2,1, -4), (3, -11,5) 
(6) (2, -1,2), (1,2,0), (4,0, -1) 
(c) (0,0,0), (1, -1,2), (-1, -2,1) 

3. Show that P, Q, / are collinear by using the distance formula in 
the following cases: 

(a) P = (1, -1, 2), Q = (0, 1, 1), R = (2, -3, 3) 

(6) P = (-1, 3, 0), Q = (3, -5, 4), R = (-2, 5, -1) 

(c) P = (1, 1, -2), Q = (-1, 0, -4), R = (5, 3, 2) 

14. Direction numbers. If L is any line in space, we may 
construct a line L' through parallel to L. The direction in 
space of L' is determined by that of L, and conversely. But the 
direction in space of L' is determined by any nonzero vector 
a, 6, c on L'. We shall therefore call any such vector a set of 
direction numbers of L. 

If PI = (xi, 3/1, 2i) and P 2 = (#2, 2/2, 22) are on L, then PI Pa 
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= (x\ #2, j/i J/2, 2i ^2) is on L f and thus the vector 
PI P2 is a set of direction numbers of L. Conversely, if 
PI = (xi, y\j Zi) is on L and a, 6, c is a set of direction numbers of 
L, then we may take (a, 6, c) = P 2 Pi, where P2 = (x\ + a, 
2/i + b y Zi + c) is a point on L. 

Every line L in space may now be prescribed by a point PI on 
L and either a second point PS or a set of direction numbers 
a, &, c. We form 



+ fe 2 + c 2 \/aH- fc 2 + c 2 



and X, M, " are the direction cosines of a ray from parallel toL. 
We shall call X, AI, v a set of direction cosines of L. Then X, 
M, v are also a set of direction cosines of L. 

Let LI and L 2 be any two lines. Then the corresponding lines 
Li and L 2 ' through intersect at and define a single acute 
angle 0, which we shall define to be the angle between L\ and Z/2. 
Let ai, &i, ci be a set of direction numbers of LI and a 2 , 62, c 2 
direction numbers of 1/2. Then we may surely use the formula of 
Sec. 10 and see that 



n 

COS V = - -=-= :-:.-_ 



Moreover, two lines are perpendicular if and only if aia 2 + 6 t 6 2 
f CiCz = 0. 

EXERCISES 

1. Find a set of direction numbers for each of the lines defined by the 
pairs of points of Exercise 1 of Sec. 13, and give a set of direction cosines 
in each case. 

2. In each of the following cases the first vector is a point PI on a line 
and the second vector is a set of direction numbers of the line. Give a 
corresponding point P 2 on the line. 

(a) (1, -1, 2), (3, 2, -1) (c) (2, 4, -1), (-4, -2, 1) 

(6) (1, 1,2), (-1, -1,3) (d) (1, -2, -1), (3, -1,2) 

3. Show whether or not the line JPiP-2 is perpendicular to the line P$P\ 
in each of the following cases: 
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. A 

(a) P 1= = (1,-1,2),P 2 = (3, -2, 1),P, = (2, -1,3), P 4 = (1,1,3) 
(6) P t - (3, 2, - 1), P, = (1, 2, - 1), P, = (3, 2, - 1), P 4 = (1, 3, 0) 

(c) P l - (2, 4, -1), P 2 = (2, 3, 1), P 3 = (4, -1, 2), P 4 = (6, 3, 3) 

(d) P! = (1, 1, 2), P 2 = (-1, -1, 1), P 3 = (1, 1, -1), 

P< = (2, -1, 1) 

(e) P t = (3, -1, 1), P 2 = (1, 1, -2), P 3 = (1, 2, -3), 
P 4 = (5,3, -I)' 

(/) P, = (2,-l,0),P 2 = (1, 1,1), P 3 = (1,0, 1),P 4 = (-1,1,1) 
(g) P 1 = (1,2, -1),P 2 = (-1,2, 1),P,= (0, -1, -2), 
P 4 = (2, -1, -1) 

4. Compute the cosine of the angle between Pip2 and PsP4 in each of 
the cases of Exercise 3 where PiP 2 is not perpendicular to PJV 

5. Show that the triangles whose vertices are P, Q, R are right tri- 
angles in the following cases, and compute the cosines of the acute angles. 

(a) P = (1, -2, 3), Q = (3, -3, 7), R = (4, 0, 2) 

(b) P = (1, -1, 0), Q = (4, 3, 5), R = (3, 0, -2) 

(c) P = (0, 2, 1), Q = (0, 0, 1), B = (-1, 1, 1) 

(d) P = (1,0, -1),Q = (2,l,0),fi = (3, -1, -2) 

6. Show that the points P, Q, /, >S are the vertices of parallelograms, 
and determine which lines are the pairs of parallel sides. 

(a) P = (2, 6, 3), Q = (3, 2, 2), R = (0, 5, 4), 5 = (1, 1, 3) 
(6) P = (2, -2, 4), Q = (1, 0, 2), R = (2, -1, 3), S = (1, 1, 1) 

7. Compute the lengths of the sides of the parallelograms in Exercise 5. 

8. Show that the points of Exercise 3 of Sec. 13 are collinear by using 
direction rather than distance. 

15. Equations of a line. If two distinct points Pi = (xi, yi, Zi) 
and Pz (x%, y^ 22) are given, we may translate our coordinate 
axes to Pi, and the x f , y f , z' coordinates of P* will be (x<i x\, 
J/2 ~ 2/i, z 2 Zi). The x 1 \ y', z' coordinates of P = (Xj y, z) are 
# ~" #i> V ~~ Vi> z ~~ z i> an d Theorem 3 implies that P is on P\Pi 
if and only if 

(19) x - xi = t(x<i - zi), y - yi == t(y 2 - j/i), 

z - Zi = /(z 2 - Zi), 

where if we direct the line PjP2 from Pi toward P% then 

(20) t 

v 



The three equations of formula (19) are called a set of parametric 
equations of the line through PI and P2. They have the property 
that if we prescribe the parameter t we obtain the coordinates 
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x, y, z of a point P on the line which actually has the property 
(20). Conversely, if P is a point on the line, its coordinates will 
satisfy these equations when t is defined by formula (20). 

Equations (19) also permit us to find the coordinates of points 
on PiP2 such that P divides the line segment PiP2 in a prescribed 
ratio. For the midpoint of PiP2 we have t = ? and thus x 
= Xi + i(x 2 Xi) = -2(22 -*- Xi). Indeed the midpoint is given 
by the formulas 

(21) X 



--, , 

ILLUSTRATIVE EXAMPLES 

I. Find the midpoint P of the line segment PiP 2 if Pi = (6, 3, -4), 
P 2 = (-2, -1,2). 

Solution 
The coordinates of P are 

6- 2 3-1 , -4 + 2 

ff _ _____ _ _ * 1 n i _ _ 1 y _ _ ! _ __ _ 1 

fr 9 ^> y 2 2 

II. Let the line PiP2 of Example I be directed from PI toward P%, and 
define / by formula (20). Find the points P(f) corresponding to t = ^, 
I, -jy, and 2, and give a rough sketch showing the relative positions of 



Solution 

We use the formula P(f) = PI + J(7' 2 - PI) = tP* + (1 - OPi and 
so obtain 

P(T) = ^(-2, -1, 2) + 1(6, 3, -4) = (V, i 2), 
P( _l) = _(_ 2 , -1,2) +2(6,3, -4) = (14,7, -10), 
p(-\~) = -_f(-2 f -1,2) +|(6,3, -4) = (10,5, -7), 

/>(2) = 2(-2, -1,2) - (6,3, -4) = (-10, -5,8). 

The following is a sketch showing the relative positions of the six points. 
- 1 - 1 - 1 - 1 



FIG. 9. 

EXERCISES 

1. Give the midpoints of each of the following line segments PiP->: 

(a) (1, -1, 1), (3, -3, 2) (e) (1, 2, -3), (-3, -2, 1) 

(6) (1, -1, 1), (0, 1, -1) (/) (4, 1, 6), (2, -1, 4) 

(c) ( -2, 1, 1), (3, -2, 2) (g) (1, 0, 2), (2, - 1, 1) 

(d) (3, -1, 2), (I, -1, -2) (h) (1, 1, 1), (2, 2, 5) 
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2. Give a set of parametric equations of PiPa in each of the cases of 
Exercise 1. 

3. Work out the problem of Illustrative Example II in each of the 
cases of Exercise 1 for the values 

/ \ * 2 3 o 3 /JA j 1 o 3 5 

(a) t = ^f, -*, -2, -T (6) t = *, -2, T , T 

4. The origin and the points Pi, Po of Exercise 1 form a triangle in 
each case. Find the remaining midpoints of sides, and give sets of 
equations in parametric form of the medians. 

16. Equations in symmetric form. If a, 6, c is any set of direc- 
tion numbers of a line through a point PI = (rci, T/I, Zi), then 
P% = (#1 + a, y\ + b, Zi + c) is on the line and formula (19) 
becomes 

(22) x x\ = ta, y ?/i = /&, 2 21 = te. 
When />c ^ 0, we may write these equations as 



the common value of the three ratios being t. These new equa- 
tions are called a set of equations of the given line in symmetric 
form. A special case of this symmetric form is that of 



x - xi _ y - T/I __ z - 



where X, M> v are a set of direction cosines of the given line. If 
PI is given, the form given by (24) is unique apart from the 
replacement when desired of X, ju, v by X, ju, ~~ v. 

Two lines are parallel if and only if they have the same sets of 
direction numbers. It is then a simple matter to use formula (23) 
to find a set of equations of a line through a given point and 
parallel to a given line. 

If one of a, b, c is zero in formula (22), we can still use formula 
(23), provided that we agree that whenever a denominator in 
formula (23) is zero we shall delete the formal ratio and set the 
numerator equal to zero; for example, the equations 

x - x l _y - y, _ z - zi 
- - 
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for be 7* shall mean 

(26) * = * ^ = *-^A 

and the equations 

/97>> x - xi _ y - yi __ z - z l 

(Z7) _ o _ 

shall mean 

(28) x = xi, y = yi. 

These results are derived from our original equations of formula 
(22). We may then treat all problems by the use of formula (23) 
and shall convert formula (25) into formula (26) and formula (27) 
into formula (28) when the sets of equations assume these irregu- 
lar forms. 

ILLUSTRATIVE EXAMPLES 

I. Give a set of equations of the line joining (3, 1, 2) to (4, 1, 1). 

Solution 

x - 3 _ y+ 1 _ z - 2 
-1 ~ -2 ~ 3 

II. Give a set of equations of the line through (3, 1, 2) parallel to 
the line joining (1, 2, 4) to (4, 1, 2). 

Solution 

x - 3 __ ?M- J _ 2 _- 2 
3 -3 2" 

III. Give a set of equations of the line joining (1, 2, 4) to ( 1, 3, 4). 

Solution 
The formal set of equations is 

x *~ * _ I/ ~ 2 _ z 4 
_____ _-, ___ = __ 

and these equations must be replaced by 



IV. Show that the four points Pi = (1, -1,2), F 2 = (2, 1,3), 
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PS = (4, 1, 2), and 7^ = (5, 1, 1) are the vertices of a parallelo- 
gram, and indicate which lines are the parallel sides. 

REMARK: A translation of axes to Pi replaces our vectors by Pi Pi, 
P 2 - Pi, P 3 - Pi, P 4 - Pi. Then P 2 - Pi, P 3 - Pi, P 4 - Pi must 
be nonzero vectors, one of these vectors must be the sum of the remain- 
ing two and these latter two vectors must not have proportional coordi- 
nates. We have evidently used the parallelogram law for addition of 
vectors. 

Solution 

Compute P 2 - Pi = (1, 2, 1), P 3 - Pi = (3, 0, -4), and P 4 - Pi 
= (4, 2, -3) = (P 2 - Pi) + (P 3 - Pi). This shows that 1\ P 2 , P.,, 
P 4 are the vertices of a parallelogram in which PiP 2 is parallel to PzP* 
and PiP 3 to PaP 4 . 

EXERCISES 

1. Write a set of equations in parametric form for each of the lines 
joining the following pairs of points, using formula (23) with the first 
point given as Pi. 

(a) (2, -1, 3), (3, 1, -2) (/) (1, -1, -2), (2, -2, -3) 

(6) (2, 1, 2), (2, -1, -2) (0) (2, 1, 3), (2, -1, 2) 

(c) (3, -1, 2), (-1, 1, 1) (K) (4, 1, 0), (-1, 1, 2) 

(d) (0, 1, -1), (1, 1, -1) (i) (1, 2, -3), (2, 2, -4) 

(e) (1, 2, 3), (3, 2, 1) 0) (1,4, -2), (-1,3, -1) 

2. Write a set of equations of a line through the origin parallel to the 
line PiP? for each case of Exercise 1. 

3. Write a set of equations of a line through PI parallel to the line 
joining the origin to P 2 for each case of Exercise 1. 

4. Write a set of equations of a line through ( 1, 1,2) parallel to the 
line PiP 2 for each case of Exercise 1. 

5. Show that the following sets of four points are sets of vertices of 
parallelograms : 

(a) (2, -1, 3), (3, -2, 5), (4, 0, 2), (3, 1, 0) 
(6) (1, -1, 2), (2, 3, -4), (2, I, -1), (1, 1, -1) 

(c) (1,2, -3), (1,2, -1), (2, -1,3), (2, -1,5) 

(d) (0, 1, -2), (1, 3, -3), (3, 6, 1), (2, 4, 2) 

(e) (3, 2, -2), (4, 0, 9), (5, 1, 2), (2, 1, 5) 

6. Write a set of equations for each of the sides of the parallelograms 
of Exercise 5. 

7. Write a set of equations for each of the diagonals. 

8. Write a set of equations for each of the lines parallel to a pair of 
sides and bisecting the other pair of sides. 
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9. Compute the lengths of the sides and diagonals of the parallelo- 
grams of Exercise 5. 

17. Distance between a line and a point. We require a pro- 
cedure for computing the unsigned distance between a line and a 
point not on the line. Let the line L be given by a set of equations 

x - xi _ y - 7/1 __ z - zi 



and the point be PI = (# 2 , ?/ 2 , z 2 ). Then we require |PP 2 | where 
P is the foot of the perpendicular from P 2 to L. If is the angle 
between PJ\ and L, then |PP 2 | = \PiP\\ * |sin 0| and thus 

|PP 2 | 2 = |P 2 Pi| 2 (sin 2 0) = [(x l - xtf + (y, - ?/ 2 ) 2 + (zi - z 2 ) 2 ] 

(1 - cos 2 0). 
But 

= r--r [X(X 2 - Xi) + M(2/2 - Vl) + 1>(Z2 



It follows that |PP 2 | may be computed by the use of 

(29) |PP 2 | 2 = (xi - x 2 ) 2 + (2/1 - y 2 ) 2 + (zi - z>)* 
- [\(xi - x 2 ) + M(?/I ~ 2/2) + 



If direction numbers are given rather than direction cosines, it is 
necessary to compute direction cosines before formula (29) can 
be applied. 

EXERCISES 

Compute the distances between points and lines in the following cases: 

(a) P,- (0,0,0);^ = ^=^ 
(6) P,= (1,0,0);*"- = V ~ - = ~ 

(c) P 2 = (l, -1,0);- = 7 = 



(d) P 2 = (0, 1, 2); 3z = 2)/ = 62 + 12 

(e) P 2 = (1,1, -1);2* = 7/- l=2e-4 

( r\P -t 11 IV g " 1 _ +J _ z ~ 2 
(J) * 2 ( 1, i, i), 2 ~ 3 ~ l 

M P - re 7 2V x + 2 - ^il - 1- 2 

w r 2 w, /, ^;, 



CHAPTER 2 
PLANES 

1. The normal form. If we draw a ray from the origin O per- 
pendicular to a given plane, the ray is called the ray normal to 
the plane and the line of which it is a part the normal line. The 
ray will have direction cosines X, /z, v, which are uniquely deter- 
mined unless the plane passes through the origin. But in this 
case the selection of X, ju, v, rather than X, ju, v, will not 
affect our results. 

Let p be the distance from to the given plane so that p ^ 0. 
Then the normal will intersect the plane in the point P 
(pX, pp., pv). A translation of axes, which moves the origin to 
Po, replaces the coordinates x, y, z of the arbitrary point P by its 
#', y', z' coordinates 

(x - pX, y - p/z, * - pv). 

The ray through P and the vector whose transformed coordinates 
are (X, /z, v) is normal to the given plane. Then P is on the plane 
if and only if the vector just defined is orthogonal to PoP, that is, 
if and only if (X, M, v) (x pX, y pp, z pv) = 0. We com- 
pute this inner product and see that P is on the given plane if and 
only if \x + py + vz - p(X 2 + M 2 + ^ 2 ) = 0. Since X 2 + M 2 
+ u 2 = 1, we have proved that P (#, y, z) is on the given plane 
if and only if 

(1) \x + y + vz = p. 

We shall call formula (1) the equation of a plane in normal form. 
It is completely unique when p > 0. When p == 0, the equation 
\x M?/ vz = is equally valid and the selection of one set 
of direction cosines rather than their negatives may be regarded 
as implying that a positive direction on the line normal to the 
plane has been selected. This is of little importance. 

EXERCISE 

Give the equation in normal form of the two planes p units from the 
origin and normal to the line joining Pi to Pi in the following cases: 

24 
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(a) p = 5, P l = (-1, 2, 3), P 2 - (1, 2, -1) 
(6) p = 2, P! = (-2, 1, -1), P 2 = (1, 3, -2) 

(c) p = 1, Px = (0, 0, 0), P 2 = (2, 2, 1) 

(d) p = iP 1 = (2,l,4),P 2 = (0, -1,2) 
() p = 7, P! = (4, -1, 6), P, = (6, -3, 8) 

2. The general equation. Every linear equation 

(2) az + by + cz + d = 0, 

in which a, b, c, d are real numbers and a, b, c are not all zero, is 
an equation of a plane; for we can write this equation in an 
equivalent form 

(3) tax + eby + ecz = \d\, 

where e = 1 if d ^ and e = 1 if d > 0. The number 



(4) t = Va r +W+~c'- > 0, 

and we may divide formula (3) by t to obtain 



A point P (x, y, z) satisfies the equation of formula (2) if and 
only if it satisfies the equation of formula (5). But then formula 
(2) is an equation of the plane for which 

(6) p = M, x-, = = , 

i.e., the plane p units from the origin and having a normal ray 
whose direction cosines are the numbers X, /z, v of formula (6). 
It is important to observe that the original coefficients a, 6, c 
are a set of direction numbers of the normal line. 

If abed ? 0, the general equation may be written in the form 



This equation is called the equation of a plane in intercept form 
and the numbers 

/o\ dj,d d 

(8) = , /--^ g = ~ 

are called the intercepts of the given plane. They are indeed the 
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nonzero coordinates of the three intercept points 
(6,0,0), (0,/,0), (0,0,0) 

where the given plane cuts the coordinates axes. 

If a plane defined by formula (2) passes through a point 
Pi = (xi, y\, zi), then ax\ + byi + cz\ + d = and 

(9) a(x - xi) + b(y - y,} + c(z - zi) = 0. 

This equation involves the coordinates of PI and a set of direction 
numbers of the normal line. We shall call it the paint, direction 
number form of an equation of a plane. 

EXERCISES 

1. Use formula (9) to write an equation of a plane through Pi per- 
pendicular to the line joining P 2 to PS in each of the following cases: 

(a) P! = (-1, 2, 3), P 2 = (-3, 1, 2), P 3 = (-5, 4, 6) 

Ans. 2x 3y - 4z + 20 = 0. 

(6) P l = (1, -2, -1), P 2 = (4, -1, 2), P 3 = (-1, 4, -1) 

(c) P 1 = (2, -1,4),P 2 = (-3,1, -2),P 3 = (1,2, -4) 

(d) Pa = (-1, 2, -3), P 2 = (-1, 4, 3), P 3 = (0, -3, -4) 

(e) P, = (-2, 1, 3), P 2 = (4, -6, 7), P 3 = (3, -4, 5) 
(/) P 1 = (1, 4, -1), P 2 = (6, 7, 8), P 4 = (9, 10, 12) 

(g) P! = (3, 2, -2), P 2 = (-9, -8, 0), P 4 = (6, -8, 7) 

2. Find p, X, ju, v for each of the following planes: 

(a) 2x - 2y + z = 6 (e) 6y + 2x - 3z = 28 

(6) 3x - 6y + 7z = -1 (/) 12x + 4y + 62 = -49 

(c) 3 + 2y + 2z = -6 (g) x + T/ - 2z = 6 

(d) -3z + 6y + 22 = 14 (h) -x + y + 2z = -3 

3. Find the intercepts of each of the planes in Exercise 2. 

4. Find a set of values of X, /x, v for each of the following planes 
through the origin: 

(a) x + y = (d) 2x + y = 2z 

(b) 3z + 4y = (e) -x + y + 2z = 

(c) 3z - 4?/ = 62 (/) 3o: + &y + 2z = 

5. Write an equation of each of two planes parallel to a corresponding 
plane of Exercise 4 and one unit from the origin. 

6. Write a set of equations in symmetric form of a line perpendicular 
to the given plane and through the point ( 1, 2, 3) in each of the cases 
of Exercise 2. 
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3, Planes through three points. Let Pi = (xi, y\, Zi), P 2 
= (#2, yzj zz), and P 3 = (#3, 2/3, z- 6 ) be three distinct points. By 
Sec. 9 of Chap. 1 the three points are collinear if and only if 
(x* xi, 2/3 yi, 23 Zi) is a scalar multiple of (x 2 #j, 
2/2 2/1, 22 21). 

If a plane passes through Pi, its equation is a(x #1) + 
b(y 2/0 + c(z Zi) = 0. If it also passes through P 2 and 
P 3 , we have the relations 



(10) 

a(z 3 - a?i) + 6(2/3 - 2/1) + c(z 3 ~ zi) = 0. 

These two equations in a, 6, c do not have proportional coeffi- 
cients when Pi, P 2 , P 3 are not collinear, and they can be solved 
to yield a set of direction numbers (a, 6, c) = (0, 0, 0) of the 
normal to the given plane. Then formula (9) is an equation of 
the required plane. 

ILLUSTRATIVE EXAMPLE 

Find an equation of the plane through ( 1,2,3), ( 3,1,2), 
(-5,4,6). 

Solution 
Equations (10) become 

-2a - b - c - 
4a - 26 - 3c - 

Then -46 = oc, 8a = c, so that -46 = 40, b = -10a. Then 

a(z + 1) + 6(y - 2) + c(z - 3) 

= a(x + 1) - 10a(y - 2) + Sa(z - 3) = 

and x - 10t/ + 82 4- 1 + 20 - 24 = x - IQy + Sz - 3 = is an 
equation of the plane. Ans. x IQy + 82 = 3. 

EXERCISE 

Find an equation of the plane through Pi, P 2 , PS for each set of points 
listed in Exercise 1 of Sec. 2. 

4. Parallel planes. Parallel planes have the same normal line. 
Then every plane parallel to ax + by + cz + d = is obtained 
if we leave a, &, c fixed and vary d. 

If we put an equation of a plane in normal form \x + ny 
+ vz = p, we obtain all planes parallel to the given plane and 
on the same side of the origin by varying p ^ 0. The distance 
between the given plane and a second plane \x + Mj/ + vz = p\ 
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is then clearly \p p\\. If, however, we consider all planes 
parallel to the given plane and on the other side of the origin, we 
see that their equations are (\x + ny + vz) = p lw Then the 
distance between Xx + ny + vz = p and (\x + ny + vz) = p\ 
is p + pi. 

ILLUSTRATIVE EXAMPLES 

I. Find an equation of a plane which is parallel to the plane 3# 2y 
+ 4z = 7 and which contains the point P = ( 2, 3, 3). 

Solution 

The equation may be taken to have the form 3x 2y + 4z = k y 
where k is to be determined so that P is on the plane. Then 6 6 
+ 12 = = k and the required equation is 3x 2y + 4z = 0. 

II. Find the distance between the plane 3x 2y + 4z = 7 and the 
plane 3;r - 2y + 4z = 12. 

Solution 

The normal forms of these two planes are obtained by multiplying 
the equations by (9 + 4 + 16)~* = (29)""*. The distance is then 
(29)-*12 - (29)"*7 = 5(29)"*. 

III. Find the distance between the plane 3x 2y + 4z = 7 and 
the plane 60; 4?/ + 82 = 15. 

Solution 

The distance is 7(29-*) + Y(29)"* = ^(29)-* = \ \/29. 

IV. Find the equations in normal form of the two planes that are 
five units from the plane x + 2y -{ 2z + 3 = and parallel to it. 

Solution 

The normal form of an equation of the given plane is obtained by 
multiplying the given equation by (1 + 4 + 4)~* = - J and thus is 

x-2y - 2z 
3 

The two planes are then 6 units and 4 units from the origin where we 
interpret a negative distance as equivalent to the replacement of X, ju, v 
by X, n, v. The answers are then x 2y 2z = 18 for the 
value 6 and x 2y 2z = 12 for the value 4. 

EXERCISES 

1. Find an equation 6f a plane parallel to the given plane and through 
the given point in the following cases: 
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(a) 3x - 2y + z = 1, (-1, -2,1) 

(6) 2x + 2y-z = 0, (-1,1,1) 

(c) z + 2/ + 3z = 4, (-2,1,1) 

(d) 2x-y-2z = 7, (3, -1,3) 

(e) 3* + 4y - 62 = 183, (4, 5, 2) 
(/) I8x - Uy + I7z = 987, (1, 1, 0) 
(g) Ite - 7y + 6z = 111,968, (2, 3, 1) 

2. Find an equation of a plane parallel to the given plane and having 
x-intercept 3 in each of the cases of Exercise 1. 

3. Find an equation of a plane parallel to the given plane and two 
units further from the origin in each case of Exercise 1 except (e), (/), 
and (0). 

4. Find the distances between the following pairs of planes: 
(a) 2x + 2y - z = 3, 2x + 2y - z = 18 

(6) 2x - y - 2z = 6, -2x + y + 2z = 12 

(c) 3x + 2y + 6z = 14, 3z + 2y + 6z = 21 

(d) 6x - 2y + 32 = 28, 12x - 4y + 62 = 35 

(e) 18x - 6y - 9z = 35, -Go; + 2y + 82 = 21 

5. Find equations of pairs of planes four units from the first of the 
planes in Exercise 4. 

5. Distance from a plane to a point. If an equation ax + by 
+ cz + d = of a plane is given and a point P = (zi, yi, z\) is 
given, we require a formula for the distance d from the plane to P. 
We first observe that we may divide the equation of the plane by 
\/a 2 + 6 2 + c 2 to convert it to the normal form 

\x + ny + vz = p. 

Then the equation of the plane parallel to tue given plane and d 
units further from the origin is \x + py + vz = p + 5. Since 
P is on this plane, we have 

(11) 8 = Xzi + M2/i + *i ~ P- 

If 5 > 0, the point P and the origin are on opposite sides of the 
plane, and if 8 < 0, the point P and the origin are on the same 
side of the plane. 

If the given plane passes through the origin, the meaning of the 
sign of 5 is not clear, and we shall simply use the formula 



ILLUSTRATIVE EXAMPLES 

I. Find the distance from the plane 2x + 2y z = 3 to the point 
(2, -1,5). 
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SoliUion 

The distance is [2(2) + 2(-l) - (5) - 3] = -2. 

II. Find the distance from 2x + 2y z = to (2, 1, 5). 

Solution 
The distance is |i(4 - 2 - 5)| = 1. 

EXERCISES 

1. Find the distances from the following planes to the corresponding 
points: 

(a) 2x + 2y - z = 3, (5, 3, -2) 
(6) -2* + y + 2z = -2, (2, 1, -2) 

(c) 3z + 2y + 6z = 14, (1, -2, 1) 

(d) 12x-4y + 6z = 35, (-1, -1,1) 

(e) 18* - 9y - 9z = 14, (1, 1, 3) 

(/) -6* + 2 + 3z = 0, (-1, 2, -3) 

2. Find the z coordinate of a point P = ( 3, 2, 2) if the distance from 
the plane 3x + 2y + Qz = 7 to the point P is 2. 

3. Solve Exercise 2 if the distance is 4. 

6. Angle between two planes. The angle 6 between two planes 
is defined to be the angle between the rays from normal to the 
two planes. It is not uniquely defined if either plane passes 
through the origin. If the equations are \x + p,y -f~ vz = p and 
Xix + my + v\z = pi, then cos = XXi + MMi + vv\. When the 
equations are given in general form, they must first be converted 
into normal form. 

EXERCISES 

1. Compute the cosine of the angle between the following pairs of 
planes: 

(a) 2x + 2y - z = 3,3x - 2y + 6z = 14 
(6) x + y + z = 1, 2x - y - 2z = 6 

(c) x - y + z = 2, 2x + y = 4 

(d) a? + T/ = 3, 3z + 4y = 2 

2. Show that the following pairs of planes are perpendicular. 
(a) 2x + 2y - z = 3, x - 2y - 2z = 1 

(6) * + y + z = 1, x - z = 6 

(c) 3x - 27/ + 4^ = 1, 2x + y - z = 1 

7. The line of intersection of two planes. Two nonparallel 
planes intersect in a line. If the planes have as equations 
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+ biy + ciz + di = 
b z y + c 2 z + d* = 0, 

then their line of intersection is orthogonal to both of their 
normals and the direction cosines of the line satisfy the equations 



+ bin + civ = 
a 2 X + &2M + c 2 j> = 0. 



However, the simplest method of determining the line is to find 
two points on it. The procedure is illustrated below. 

ILLUSTRATIVE EXAMPLE 

Give a set of equations in symmetric form of the line L of intersection 
of the planes 

2x + y - 22 = 1, 5x + 4y - 6z = ~2 
and determine direction cosines of L. 

Solution 

We solve for x and y by writing 8x + ty 82 = 4 ; therefore 3x 2z 
= 6 and 

2s + 6 22 - 9 

x = ^ > t/ = 2^ 2z + 1 = g -- 

Put 2 = 0, and obtain P t = (2, -3, 0) on L. 

Put 2 = 3, and obtain P 2 = (4, 1, 3) on L. Then L has as its equation 

x - 2 _ y + 3 _ 
2 " 2 ~ 3 

and the corresponding direction cosines are 

223 



Vl7' VT7' 

EXERCISE 

Give a set of equations in symmetric form of the line intersection of 
the following two planes, and determine its direction cosines, 
(a) 3z - 2y + z = 5, x + y + z = 1 
(ft) s - 2y + z = 6, re - z = 2 

(c) Zx + 4y + 2z = 1, x - y + z = 6 

(d) 6s + 2y - 3s = -5, x - y + 2 = 6 

(e) 4x - j/ + 2 = 2, 2x + y - 3^ = 4 
(/) * - V - 32 = 2, 2x + 2y - 2 = -3 
(g) 4x + 3y - 2z = 1 , x + 2y - 3z = 4 
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8. Angle between a line and a plane. Let the plane ax + by 
+ cz + d = and the line 

(14) x- xi _ y - yi = s - z\ 

ai bi ci 

be given. We require the angle <f> between the line and the 
plane. This angle is defined to be the complement of the angle 
6 between the directed normal to the plane and the line and is not 
unique unless the line is also directed. When both are directed, 
we will have 



where we have now made X, M> v and Xi, pi, v\ unique. Then 
cos = sin and 

(16) sin <t> = XXi + MMI + vv\. 

EXERCISE 

Find the sine of the angle between the following lines and planes: 
(a) 3z - 2y + 4z = 1, 3x + 4 = 2y - 1 = 6z + 7 

(6) x + 2y - 2z = 2, = = 



(c) 3z + 2y - 62 = 4, = = 2-3 

9. Pencils of planes. Let the equations f(x, yS$) = a\x 
+ biy + CiZ + di = and g(x, y, z) = a 2 x + b*y + c& + d* = 
define two distinct planes. Then the equation 

(17) sf(x, y, 2) + tg(x, y, 2) = 

is called the equation of the pencil of planes determined by the 
two given planes. Indeed, formula (17) is an equation of a 
plane for all real numbers s and t which are such that the vector 
s(a,i, 61, Ci) + t(a 2 , 6 2 , 02) 5^ 0. When the two given planes are 
not parallel, they intersect in a line L; formula (17) represents a 
plane Q(s, <) for all real numbers s and t not both zero, and every 
such plane Q(s, ) contains the line L. For a point Po = (#o, 
j/o, So) is on L if and only if both / = f(x Q , y , 2 ) = and g$ = 
g(x Q , t/o, o) =0; whence s/ + ^o = and P is on Q(s, t). 

When the two given planes are parallel, there exists a nonzero 
real number k such that (a^ 62, 0%) = /b(ai, 61, Ci) but d^ ^ kdi } 
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since the planes defined by /(x, y, z) = and gr(x, y, z) = have 
been assumed to be distinct. In this case formula (17) becomes 

(18) (s + kt}(aix + 6ij/ + ciz) + sdi + td, = 0, 

and this equation represents a plane for all real numbers s and t 
such that s 7^ kt. All such planes are evidently parallel to 
the two given planes, and we call the pencil defined by the two 
given planes a parallel pencil. 

Conversely, let Q be any plane through the line of intersection 
L in case the given planes intersect and parallel to f(x, y, z) = 
in the case of a parallel pencil. Then Q is uniquely determined 
by the property just described and the assumption that Q con- 
tains a point Po == (#o, 2/o, 2o) not on both of the given planes. 
The corresponding values / = /Or , yo, z ) and g Q = Q(XQ, yo, z ) 
determine a member Q(<7o, /o) of the pencil, and the given plane 
Q actually coincides with Q(gr , /o). For s/ + tgo = 0, when 
s = #o and / = / , and P is on Q(<7o, /o). It then remains 
only to see that the given values of s and t actually define a plane. 
This is true in the case of intersecting planes, since in that case 
we merely require that s and t should not both be zero, and we 
have already assumed that /o and g Q are not both zero. In the 
case of parallel planes, the hypothesis that s + kt = implies 
that </o = fc/o. Then /o cannot be zero, and formula (18) for 
(x, y, z) replaced by (X Q , y , z ) yields sdi + id* = g<4i fod* 
= fo(kdi d 2 ) =0, d 2 = kdi contrary to hypothesis. Thus 
we have shown that the pencil of planes defined by two distinct 
planes contains all planes through their line of intersection when 
they intersect, and contains all planes parallel to them when they 
are parallel. 

When two- given planes intersect, their equations may be put 
into normal form and thus written as 



(19) \ix + fjny + viz - pi = 0, \ 2 x + M2?/ + v& - p 2 = 0. 

Then the pencil of planes defined in this case contains the two 
planes whose defining equations are 

(20) (Xi X 2 )x + (MI M2)y + (vi ~ vj)z (pi p 2 ) = 0, 

(21) (Xi + X 2 )z + (MI + M2)y + (vi + v*)z - (pi + p t ) 0. 

Each of these planes consists of points P = (x, y, z) whose dis- 
tances d\ = \ix + my + v& pi and d* = X 2 # + /z 2 y + V2Z 
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pz from the two given planes are equal in absolute value. 
The planes are then the bisector planes of the dihedral angles 
formed by the two given intersecting planes. 

ILLUSTRATIVE EXAMPLES 

I. Find an equation of the plane containing the line of intersection of 
the planes 3x - 2y + 4z - 6 = 0, 2x + 3y + z - 12 = as well as 
the point (-1,2,4). 

Solution 

We compute /(-I, 2, 4) = -3 - 4 + 16 - 6 = 3, and 0(-l, 2, 4) 
= -2 + 6 + 4 - 12 = -4. This yields the values s = 4, t = 3, and 
we compute I2x - 8y + 16z - 24 + 6x + 9y + 3z - 36 to obtain 
the solution I8x + y + 19z - 60 = 0. 

II. Find an equation of the plane containing the line of intersection 
of the planes of Example I and normal to a line whose direction numbers 
(a, 6, c) have the property a = 36. 

Solution 

The coefficients of the required plane are a = 3 + 2t, b = 2 + 3, 
c = 4 + , d = -6- I2t. Then a = 3 + 2t = 36 = -6 + 9, 7t = 9, 
t = T. Multiply by 7 and obtain 2lx - Uy + 28z - 42 = I8x + 27y 
+ 9z - 108 = 0, and the answer is 39z + 13?/ + 37z - 150 = 0. 

III. Find the bisectors of the dihedral angles formed by the planes 
2x + y - 2z - 4, 2x - 3y + 62 = -2. 

Solution 
The required equations are 

2x + y -2z - 4 2x - 3y + 62 + 2 

- -- = '~ - > 



that is, 14z + 7y - Uz - 28 = (6x - 9y + 18 + 6). Hence, the 
solutions are 80: + 1 677 - 32z - 34 = and 20z - 2y + 4z - 22 = 0. 

EXERCISES 

1. Find an equation of the plane containing the line of intersection of 
the two given planes and the given point in the following cases: 
(a) 2x - 2y + z = 3, x + 2y + 2z = 5, P = (-1, 2, -2) 
(6) 3z + 2y + 6z = 1, 2x - &y + Zz = -4, P = (2, -1, -3) 

(c) a; - 2y - 2z = 3, 3x - 2y - 62 = 4, P = (-1, 0, 1) 

(d) 3x + 67; + 2z = -3, 2x + 2y - 2 = -2, P = (1, 1, 1) 

(e) ^ + y + z = -4, 2x - y - * = -3, P = (-3, 2, -1) 

(f) x + z= -2,y -2z = 4, P = (1, 2, -3) 
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(g) x + y + z + 1 = 0, 2* - 3y - 3z - 3 = 0, P = (3, -1, -2) 
(h) x = y + z, y - x + z, P = (4, -1, 2) 
(i) 2z + 4y - 82; = 7, 4z + 2t/ - 3z = 9, P = (-1, -1, -2) 
(j) 3z - 4y + 6z = 1, 5z - 2y + z = 8, P = (0, 0, 0) 

2. Find an equation ax + by + cz + d = of the pencils defined by 
Exercise 1 in which 2a = 36, and find one in which 4d = 3c. 

3. Find the bisectors of the dihedral angles formed by the pairs of 
planes of parts (a), (6), (c), (d), (e), (/), and (h) of Exercise 1. 




FIG. 10. 

10. Parametric equations of a plane. Three points Pi 

= Oi, 2/1, zi), P 2 = (x 2 , 2/2, 2 2 ), and P 3 = (x 3 , 2/3, z 3 ) are not 
collinear if the vectors 



P 2 ' = P 2 - P! = x 2 - 



j/ 2 ~ 



2 2 - zi), P/ = P 3 - PI 
== (z 3 ~ Xi, 2/3 2/i, 23 - 21) 

do not have proportional coordinates. Then PI, P 2 , P 3 determine 
a plane S. Translate the origin to 0' = Pi, and thus write 
P 7 = (x xi, y 2/1, z Zi) for the vector of the coordinates 
of any points P = (x, y, z). If s and are any real numbers, the 
line joining 0' to P 3 ' passes through <P 3 7 , Hence, sP 2 ' and /Ps 7 
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are in the plane S and their sum sP 2 ' + tP$ is the end of the 
diagonal of a parallelogram whose vertices are in /S; sP 2 ' + tP* 
is the vector (x Xi, y 3/1, z Zi) of a point of S. 

Conversely, if P is a point of S, we draw a line L 2 in S through 
P parallel to PiP 2 and a line L 3 in S through P parallel to PiPs. 
Then L 3 cuts PiP 2 in a point whose transformed vector is sP* 
and L 2 cuts PiP 3 in a point whose transformed vector is tfY. 
For a typical case, see Fig. 10. But then P f = sP 2 ' + P 3 ', and 
we have 

x xi = s(x 2 xi) + t(x* xi) 

(22) y - 2/1 = 5(2/2 - 2/i) + (2/3 - y\) 

z - Zi = sfa zi) + t(z* - 21). 

We have shown that the equations of formula (22) give the 
coordinates x, y, z of points P on the plane S through PI, P 2 , PS 
for all real numbers 5 and i and that every point on the plane S 
has coordinates which satisfy these equations for some real num- 
bers s and t. We are then justified in calling the equations of 
formula (21) a set of parametric equations of the plane S. 

EXERCISE 

Give a set of parametric equations for the plane through PI, P 2 , PS for 
each set of points listed in p]xercise 1 of Sec. 2. 



CHAPTER 3 
SURFACES AND CURVES 

1. Equations of a surface. A real single-valued function 
f(x, y, z) of three real variables x, y, z is a correspondence wherein 
every triple of real numbers xi, y\, z\ uniquely determines a real 
number /(#i, y\, Zi) called the value of the function at the point 
where x = Xi, y = y\, z = Zi. The functions commonly con- 
sidered in analytic geometry are usually defined by algebraic or 
trigonometric formulas. 

If f(x, y y z) is a real single-valued function, a point P 
(xi, yi, Zi), defined by real coordinates x\ 9 yi, Zi, is said to be a 
solution of the equation 

(1) f(x, y, z) = 

if f(xi, yij Zi) is the real number zero. An equation (such as 
x 2 + y 2 + z 2 = 0) may have only a single point as a solution or 
no real points as solutions. However, the solutions usually make 
up a geometric locus (i.e. f set) of points called a surface. 

The equation f(x, y, z) = is said to be an equation of a sur- 
face >S if S is the set of all solutions of the equation. A surface, 
such as a cone, a cylinder, or a sphere, may be defined geo- 
metrically, and we shall give geometric definitions of certain 
special types of surfaces. If we then give an equation f(x, y, z) 
= and state that it is an equation of S, the proof of the result 
requires that we verify the following two properties: 

1. Every point P on S is a solution of f (x, y, z) = 0. 

2. Every solution P of f (x, y, z) = is a point on S. 

We shall frequently define a surface S by prescribing an equa- 
tion f(x, y, z) = and stating that S is the surface of all points 
P which are solutions of the equation. As we have seen, the 
solutions might consist of a single point or no real points and the 
corresponding surface will then be called a point surface or an 
imaginary surface. The general geometric definition of those 
point sets which are true surfaces is beyond the scope of the 
present text. 

37 
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If a is any nonzero real number and /(#, y, z) = is an equa- 
tion of a surface, then af(x, y, z) = is also an equation of the 
same surface. Indeed, let g(x, y y z) be any function such as 
x 2 + 1 which has the property that g(xi, yi, Zi) ? for all real 
numbers x\, y\, z\. Then if /(z, y, z) = is an equation of S, so 
is 0(3, y, z)f(x, y, z) = 0. For any product 00i, y l9 Zi)f(xi, yi 9 Zi) 
is zero if and only if f(xi, y\, Zi) is zero. 

It will be convenient for us to speak of "the surf ace f(x, y y z) 
= 0" when we mean the surface defined by the equation f(x, y, z) 
= 0. Thus we may speak of "the plane ax + by + cz + d = 0," 
and we have already used this terminology in Chap. 1 and a 
similar terminology for lines in Chap. 2. 

2. Space curves. The intersection of two surfaces Si and 82 is 
the locus of all points on both Si and 82. This locus may be a 
surface; for example, the equation xy = is an equation of a 
surface S which consists of two planes, and the equation xz = 
is the equation of a surface which also consists of two planes. 
Their intersection is the plane x = 0. 

In most of the cases commonly considered in analytic geometry 
the intersection C of two surfaces Si and 82 is a space curve. 
This is a geometric manifold with one degree of freedom and the 
precise definition is beyond the scope of the present text. We 
shall assume that the equations 

f(x, y, z) = 0, g(x, y, z) = 

are equations of two surfaces and shall call the pair of equations 
a set of equations of the space curve which is their intersection. 
The intersection may then turn out to be empty, a single point, 
a genuine space curve, or a surface, and a revision of the geo- 
metric term space curve, by means of which the intersection has 
been described, may be necessary. Note that the curve of inter- 
section of two planes is a straight line and that a straight line is 
an instance of a space curve. 

As in Sec. 1, we shall speak of "the curve /(z, y, z} = 0, 
g(x, y, z) = 0" when we wish to refer to the curve of intersection 
of the surfaces defined by the equations f(x, y, z) =0 and 
g(x, y, z) = 0. 

3. Plane sections. The curve of intersection of a plane and a 
surface is called a plane section of the surface. In general, it is 
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an ordinary plane curve, but it may consist of a single point or 
of no points. 

The determination of a plane section may be carried out most 
simply if an equation of the surface can be obtained relative to a 
coordinate system in which the given plane is one of the coordi- 
nate planes. For example, if the given plane is the x, y plane, and 
a corresponding equation of the surface is/(#, y, z) = 0, then the 
plane section is the curve f(x, y, 0) =0 in the x, y plane. 

When a surface is given by means of an equation /(#, y, z) = 
and a plane by an equation ax + by + cz + d = 0, it is always 
possible to perform a transformation of coordinates so as to 
obtain an x r , y', z f coordinate system in which the given plane 
becomes the plane z f = and the surface is g'ven relative to the 
transformed coordinate system by an equation <t>(x', y', z'} 0. 
The plane section may then be analyzed by a study of the equa- 
tion 0(x', T/', 0) = 0. The equations that transform f(x, y, z) 
into 4>(x', y', z'} will be developed in Chap. 7. 

ILLUSTRATIVE EXAMPLE 

Discuss the sections of the surface 

** + ?/! , I 2 _ i 

16 "*" 4 "*" 9 

made by planes parallel to the x 1 y plane. 

Solution 

The required planes are the planes z = fc, where & is a real number. 
Then the sections are curves whose equations are 

* 2 9 ~ fc2 



. 
16 4 ~ 9 

If k* > 9, the corresponding plane does not intersect the surface. 
When k = 3, the intersection is the single point (0, 0, 3) and the plane 
is said to be tangent to the surface. Similarly, if k 3, the plane 
z = 3 is tangent to the surface at the point (0, 0, 3). If k 2 < 9, 
the corresponding sections are ellipses having center at (0, 0, &), semi- 
major axes ^ \/9 & 2 , semiminor axes f \/9 A; 2 , and a major axis, 
which is a segment of the line z = fc, y == 0, parallel to the x axis. 

EXERCISE 

Discuss the sections of the following surfaces made by planes 
parallel to the x, y plane: 
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() 1 + I + = 1 (*) ^ 2 " 2/ 2 - * 2 = 1 

(fr) ? + 1 " 15 = w 4 * 2 + 22/2 + 4 * 2 = * 

(c) x 2 = 4^ ((/) 4^ = ?/ + 2y 






(d) fg + ~ = 1 (A) z 2 + 2x + 4y* - z 2 + 2* = 

4. Algebraic surfaces. A surface is called an algebraic surface 
if it is defined by an equation f(x, y, z) = 0, where f(x, y, z) is a 
polynomial in x, y, z with real coefficients. The curve of inter- 
section of two algebraic surfaces is called an algebraic curve. 

Every polynomial /(z, y, z) is a sum of terms of the form 

ax'y'z*, 

where the coefficient a is a real number and r, s, t are natural 
numbers, i.e., nonnegative integers. The degree of such a term 
is r + s + t, and the degree n of f(x, y, z) is the largest degree of 
all of the terms of f(x, y, z) that have nonzero coefficients. We 
call f(x, y, z) a constant polynomial if n = 0. If n = 1, the 
polynomial f(x, y, z) = ax + by + cz + d and the equation 
/(x, i/, 2) = is an equation of a surface that we have called a 
plane. If n = 2, we call /(x, y, 2) = a quadric surface, and if 
n = 3 we call /(x, y, 2) = a cubic surface. 

If /(x, y, 2) = g(x, y, z)h(x, y, z), where g(x, y, z) and h(x, y, z) 
are polynomials with real coefficients and neither is a con- 
stant, we call f(x, y, z} a reducible polynomial. Otherwise, we 
call f(x, y, z) an irreducible polynomial. In the former case 
f(x, y, z) =0 is an equation of the surface consisting of all solu- 
tions of either of the equations. 

g(x, y, z) = 0, h(x, y, z) = 0. 

If either of these surfaces is an imaginary locus, such as x 2 + y 2 
+ z 2 + 1 = 0, the factor is deleted and /(x, y, z) is replaced by 
the remaining factor. However, when both surfaces contain 
real points, we call the surface defined by /(x, y, z) = a reducible 
surface. Then every irreducible surface is defined by an equation 
f(x, y, z) = 0, where /(x, y, z) is an irreducible polynomial. 

6. Cones. A cone is a surface S containing a fixed point V 
called the vertex of S such that all points of any line joining F to a 
point of S are points of S. We shall restrict our attention to 
algebraic cones, i.e., cones which are algebraic surfaces and shall 
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assume that n > I. Thus we do not include planes in the set of 
surfaces called cones. 

A polynomial f(x, y, z) is said to be homogeneous if all terms of 
f(x, y, z) have the same degree n. A homogeneous polynomial of 
degree n is called an nic form. We shall study quadratic forms, 
i.e., forms of degree two, later. If f(x, y, z) is homogeneous of 
degree n, then f(tx, ty, tz) = t n f(x, y, z). Also /(O, 0, 0) = if 
n > 0. 

Theorem 1. // f (x, y, z) has the property that g(x', y', z') 
= f (x 7 + XQ, y' + yo, z' + ZQ) is homogeneous of degree n > 1 in 
x', y', z' then f (x, y, z) = is an equation of a cone with vertex at 
(x , y , z ). 

For a translation of axes, which moves the origin to V = 
(XQ, ?/o, zo), replaces /(z, y, z) by g(x' , y', z'}. Since g(Q, 0, 0) = 0, 
the point V is on the surface S defined by f(x, y, z) = 0. If 
Q = (xi, yi 7 , zi) is a point on S, then g(xi, yi, Zi) and 
g(txi', tyi', fe/) = t n g(xi', ?/i 7 , */) = for every t. But t(xi', 
VL, Zi) = (x 1 ', y 7 , z 7 ) is the vector form of a set of parametric 
equations of the line joining V to Q, and we have proved that all 
points of this line are on S, S is a cone. 

The cones, which consist of all points on the lines joining the 
vertex to the points of a fixed plane curve, are of particular 
interest in geometry. Let us then derive the equations of such 
surfaces in the case where the given curve is a curve in the x, y 
plane and its equation is <t>(x, y) 0. Suppose that PO = 
(XQ, ?/o, ZQ) is the vertex. Then a point P = (x, y, z) is on the 
cone if and only if P P = t(P\ PO), where Pi = (x\, y\, 0) is 
a point on the curve. Consequently 

(2) X - XQ = t(Xi - XQ), y - 7/o = %l ~ 2/o), 

Z ZQ tZQ. 

We assume that P is not a point of the x, y plane and thus 
Z Q 7* 0. Then we may eliminate t and obtain 

ZQ(X XQ) = (Xi - XQ)(ZQ - z) , Z Q (y y ) = (j/1 ~ J/o) 

(ZQ ^ Z). 

We solve for Xi and y\ to obtain 

/ox ^ , ZQ(X - XQ) __ z Q (y - y ) 

(3) xi = x H 9 7/1 = t/o -\ -; 

ZQ Z ZQ Z 
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and substitute these formulas in <t>(xi, y^ = 0. The resulting 
equation in x y y, z is an equation of the given cone. 

ILLUSTRATIVE EXAMPLES 

I. Find an equation of the elliptic cone with vertex at (2, 1, 3), and 
base the curve 4# 2 + y 2 = 1. 

Solution 
We use formula (3) and write 

3(x - 2) 3(y + 1) 

- 



Then (3 z)zi = 3# 2z, (3 0)1/1 = 3y + z, so that the required 
equation is 

4(3* - 22) 2 + (3ff + 2) 2 = (3 - 2) 2 . 

The solution should be left in this form since the equation is not essen- 
tially simplified when the indicated operations are carried out. 

II. Find an equation of the hyperbolic cone through the point (3, 0, 0) 
and the hyperbola whose equations are z = 0, (x 2 /9) (?/ 2 /4) = 1 . 

Solution 

We use the method of derivation of formula (3) to write x 3 = - 3f , 
y = ty,, z = tzi. Then -3</ = (* - 3)f/i, -32 = (a; - 8)2,, and 



(x - 3) 2 2! 2 (x - WyS _ W 97/' 

9 4 -(*-<*) - -9 4- 

Hence 4(x 3) 2 = 43 2 Oj/ 2 is an equation of the cone. 

EXERCISES 

1. Find an equation of the cone with vertex at the given point P and 
passing through the curve whose equations are 2 = 0, <j>(x, y) = in the 
following cases: 

(a) (0, 0, 3), z 2 + y* = 9 (c) (1, -1, 2), z 2 + y* = 1 

(6) (0, 0, -2), x* + ?/ = 4 (d) (-1,1, -3), x 2 - 2/ 2 = 1 

(?) (-1,2, -l),3* + 2y= 1 

(/) (-1, 2, 3), * 2 + 2s + t/ 2 - 4v + 3 = 

(g) (2, I, -4), z 2 - 4z + 3j/ 2 - 6y + 6 = 

(A) (-1, 2, 3). * 2 = 8y (j) (2, 1, 4), x* - 37/ 2 = 1 

(t) (0, -1, 2), 2* 2 - j/' = 1 (*) (0, 0, 1), x 3 + y 3 = 1 
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2. Interchange x and z in the equations of the plane curves of Exer- 
cise 1, and leave the points P unaltered. Find equations of the corre- 
sponding cones. The equations are then x = 0, $(z, y) 0. 

3. Find an equation of the cone with vertex at ( 1, 1, 2) and passing 
through the curve of intersection of the surface x 2 + 4t/ 2 = 1 and the 
plane z = 3. 

6. Cylinders. A cylinder is a surface consisting all of the 
points on all the lines which are parallel to a given line and which 
pass through a fixed plane curve in a plane not parallel to the 
given line. It is a simple matter to obtain an equation of a 
cylinder with respect to a coordinate system which is chosen so 
that the plane curve lies in a coordinate plane. 

Let us assume that the curve lies in the x, y plane so that its 
equations are 2 = 0, ^(x, y) = 0. The fixed line will have direc- 
tion numbers a, 6, c, and a point P = (x, y } z} will be on the 
cylinder S if and only if 

m x ~ Xl - y ~ yi - z , 

w 5 ~c' 

where (xi, y\, 0) is on the given curve. Then 

/c\ a & 

(5) xi = x - -z, y l = y - -z, 

c/ c 

and 

(6) 

is an equation of the surface. Note that our hypothesis that the 
given line is not parallel to the given plane implies that the 
number c is not zero. 

If a = b = 0, the given line is parallel to the z axis and the 
equation we have derived reduces to <l>(x, y) = 0. Thus an 
equation in two of the three variables x, y, z represents a cylinder 
generated by a line parallel to the coordinate axis corresponding 
to the missing variable. 

ILLUSTRATIVE EXAMPLES 

I. Find an equation of the elliptic cylinder determined by the curve 
z = 0, 4.r 2 + y 2 = 1 and a line with direction numbers (2, 1, 3). 

Solution 
We use formula (6) to write the answer 

4(* - fz) 2 + (y + z) 2 = 1. 
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II. Find an equation of the hyperbolic cylinder determined by the 
curve x 5, 9z 2 4;/ 2 = 1 and the line with direction numbers (1, 
-2, 3). 

Solution 

As in formula (4), the points on the cylinder satisfy 

a? 5 _ y - y\ z z\ 
. ~T~ " ~~=2 ~ 3 ' 

where (5, y\, 21) is on the curve. Then 

yi = y + 2(x - 5), z l = z - 3(z - 5) 

and thus the required equation is 

9(2 - 3x + 15) 2 - 4(?/ + 2x - 10) 2 - 1. 

III. Find an equation of the cylinder determined by the curve of 
intersection of the plane 2x + 3?/ z = 1 with the surface 3z 2 ?/ 2 
+ 2z 2 = 1 and the line with direction numbers ( 1, 2, 3). 

Solution 
The points P = (#, T/, 2) satisfy 

x xi _ y - ?/i _ - gi 

" 



-1 " 2 ~ 3 f 

where cci, t/i, zi are on the curve. Then zi = 2 3(#i x), T 
y 2(^i a:) and 2#i + 3?/i zi = 2^i + 3?/ 6(^1 ) 
3(;ci - a?) = 1. Hence, xi = 3o; + 3y - z 1, Zi = z + 3x 
+ 3t/ - 2 1) = -60: - Jty + 4^ + 3, 7/1 = y + 2x - 2(3z + 3y z 
- 1) = -4x - 5y + 2z + 2. Then our solution is 3(3z + 3y - z - 
I) 2 - (4x + 5y -2z - 2) 2 + 2(6* + 9?/ - 4z - 3) 2 = 1. 

EXERCISES 

1. Find an equation of a cylinder determined by the given plane 
curve and set of direction numbers in each of the following cases: 
(a) x* + y* = 9, z = 0, (-1,2, 3) 
(6) *-y-4,-2 f (2,-l,4) 

(c) 2y-*= 1,*- -3, (-!,-!, 1) 

(d) 3z' + 2 2 ^ = 1, y = 3, (1, 2, 0) 

(e) * 2 - By, z = 0, (0, 1, 2) + 

(/) W * = 4?-2,* = 0, (-1,3,2) 
(g) 2 2 = 4x^ + 8*, y = 0, (0, 4, 1) 
(A) * + *- 1,0 = 0, (-1,1, 1) 
(t) *-y= 1, = 0, (-2,1, -1) 
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(j) * 3 = s 2 ,y = 3, (-2,2,1) 

(A) </ 2 -z 2 - -1,* = 2, (-2,2,2) 
(Z) x + z 2 = 1, y = 0, (1, 1, 1) 

(m) x-y = *j/, = 0, (1,0,1) 

(n) sy = -2,* = 1, (-2,1,1) 

2. In the following cases, the plane curves are not in planes parallel 
to coordinate planes. Find an equation of each of the corresponding 
cylinders. 

(a) z 2 + 7/ 2 = !,* = *, (-1,1, 2) 

(6) * 2 + s,* + * = 1, * = 2?j, (1, 2, 3) 

(c) 3z 2 - 4?/ + z 2 = 1, x - y = z, (1, -1, 1) 

(d) z 2 = ?/ + z, z = z - 2y, (-2, 1, 1) 

7. Surfaces of revolution. If a plane curve (or a line) is 
revolved about -a line in the plane of the curve, the resulting sur- 
face is called a surface of revolution. We shall obtain the equa- 
tions of such surfaces relative to a coordinate system where the 
given plane is the x, y plane and the given line is parallel to the 
x axis. 

Let <t>(x, y) = 0, z = be equations of the curve and z 0, 
y k be equations of the line. Then a point P = (x, y, z) is on 
the surface of revolution S if and only if P is on a circle whose 
center is (x, /c, 0) and whose radius is the distance \yi k\ from 
(x, k, 0) to the point (x, yi, 0) on the given curve. Thus P = 
(x, y, z) is at the same distance from (x, fc, 0), and 



\yi - fc| = V(x - xY + (y - A;) 2 + z 2 = V(y - A;) 2 + z\ 

It follows that (T/I k) 2 = (y fc) 2 + z 2 and that the required 
equation of the surface S is obtained by replacing (yi k) 2 by 
(y k) 2 + 2 in <t>(x, y) = or in a suitably modified equation. 
For the modifications, see the examples below. 

ILLUSTRATIVE EXAMPLES 

I. Find an equation of the surface of revolution obtained by revolving 
the curve z = 0, x 2 = 4y about the x axis. 

Solution 

If (x, y, 0) is on the curve both (x, y, 0) and (x, y, 0) are on the sur- 
face, the surface is the surface of revolution of the double curve x 4 = 16?/ 2 . 
Hence, an equation of the required surface is x 4 = 16(?/ 2 + z 2 ). 

II. Find an equation of the surface of revolution obtained by revolv- 
ing the curve z = 0, a: 2 = 2y about the line z = 0, y = 3, 
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Solution 

The equations of curve are z = 0, x 2 6 = 2(y 3), and therefore 
(# 2 6) 2 = 4(y 3) 2 may be converted into the required equation 
(x* - 6) 2 = 4[(y - 3) 2 + z 2 ], that is, z 4 - 12x 2 + 36 = 4y* - 24y + 
4z 2 + 36. Hence, the answer required is x 4 12z 2 = 4?/ 2 24?/ + 4z 2 . 

EXERCISES 

1. Find the equations of the surfaces of revolution obtained by revolv- 
ing the following curves in the x, y plane about the x axis : 

(a) z 2 + 2?/ = 1 (/) x z - 4x + 2y = 1 

(6) z 2 - y* = 4 fo) z 4 + 2x*y + y 2 = 1 

(c) 2x 2 - 3t/ 2 = 2 (A) a; 2 + 3?/ 2 = 2z?/ 

(d) a: 2 + 2x + ?/ 2 = 1 (f ) ^ 3 + y 3 - 1 

(e) x 2 + xy + y 2 1 (j) x 2 7/ 2 = # 3 ?/ 3 

2. Solve Exercise 1 for revolution about the y axis. 

3. Solve Exercise 1 for revolution about the following lines in the x, y 
plane: 

(a) x = -1 (6) y = 2 (c) 3 = 1 

8. Symmetries of surfaces. Let $ be a surface and T be a 
plane. Then $ is said to be symmetric with respect to T if every 
line perpendicular to T and cutting S in a point P not on T also 
cuts S in a point Q such that P and Q are on opposite sides of T 
and the same distance from T. If I 7 is taken to be the x, ?/ plane 
and /(#, y, z) =0 is the equation of S, then /S will be symmetric 
with respect to T if and only if f(x, y, z) =0 for every point 
P = (x, y, z) such that f(x, y, z) = 0. This is evidently satisfied 
when f(x, y, z) is a polynomial in x, y y and 2 2 , rather than z. 

A surface S is said to be symmetric with respect to a line L if 
every line perpendicular to L which cuts S in a point P not on L 
also cuts S in a second point Q not on L such that P and Q are on 
opposite sides of L and the same distance from L. If L is taken 
to be the x axis and f(x, y, z) = is an equation of S, then S will 
by symmetric with respect to L if and only if f(x, y,z)=Q 
for every point P = (x, y, z) such that /(x, y, 2) = 0. This is 
satisfied when f(x, y, z} is a polynomial in x, y 2 , 2 2 , yz. 

A surface S is said to be symmetric with respect to a point if 

every line joining to a point P on $ also joins to a point Q on 

> 
/S such that QO = OP, If is the origin, Q = P and therefore 
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we must have f( x, y, z) = for every point P = (#, y, z) 
such that /(#, ?/, 2) = 0. In the exercises below only symmetries 
with respect to planes parallel to coordinate planes and lines 
parallel to coordinate axes are to be discussed. 

ILLUSTRATIVE EXAMPLES 

I. Discuss the symmetries of the surface x 2 4# + 3y 2 + 2z 2 I2z 

__ I 

Solution 

We translate the origin so as to simplify the equation. The equation 
is (x - 2) 2 + 3y 2 + 2(2 - 3) 2 = 23, and thus the translation x' = x 
- 2, y' = y,z' = z - 3 converts the equation to x' 2 + 3?/ 2 + 2z' 2 = 23. 
This surface has symmetry with respect to the new origin, the x', y', 
and z' planes, and the #', ?/, and z f axes. Thus the original surface is 
symmetric with the planes x = 2, y = 0, and z 3, the lines of inter- 
section of these planes ; and their point of intersection (2, 0, 3). 

II. Discuss the symmetries of the surface 2x 2 + 4x + 3yz* +1=0. 

Solution 

The equation is equivalent to 2(x + I) 2 + 3yz 3 1, and thus to 
2^ /2 + 3y'z' 3 = 1, where x' = x + 1, y' = y, z' z. This equation is 
unaltered when we replace x' by x f and y', z' by y', z'. It is then 
a surface symmetric with respect to the point ( 1,0,0), the plane 
x = 1, and the line y z 0. 

EXERCISES 

Discuss the symmetries of the following surfaces: 

(a) 3x 2 + y 2 - z = (g) x 4 + 2x 2 = y 2 z 

(6) x 2 - 2yz = 1 (A) x 2 - x* = y 2 + yz + z 2 

(c) x 3 + y/ 3 = z 3 (i) x 2 = y* - 2y + z(y - 1) 

(d) x 3 = yV ( j) * 3 = (y 2 - 2y)z + z* 

(e) x 2 - 2x = 7/2 (A;) x 3 + ?/ - 2 5 = 
(/) 3a; 3 - .T = 7/2 (0 x 2 + y 2 z 2 = a: 3 ?/ 

9. Intersections of a line and a surface. Every line L is 
defined by a set of equations 

(7) a: = x () + t%, y = y + ^, 2 == 2 + # , 

where (a: , 2/o, o) is a point on L, (, ty, f ) is a set of direction num- 
bers of L, and any real value of t defines a corresponding point 
(x, y, z) of L. If f(x, y,z) =0 is an equation of a surface S, the 
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points of intersection of L and S are those points of L for which 
(8) $(*) = f(x* + , 2/o + fcf, zo + = 0. 

If /(#, y, 2) =0 is the equation defining an algebraic surface 
of degree n, the corresponding function <(/) is a polynomial in / 
tvhose degree cannot exceed n. Then either $() is identically 
zero and all values of t satisfy <f>(t) = or <t>(t) = has at most 
n real roots. In the former case, all points of L are points of S, 
and the line L lies wholly on the surface S. We have proved the 
following result: 

Theorem 2. // a line is not wholly on an algebraic surface of 
degree n, it cuts the surface in at most n points. 

A line then cuts a quadric surface in no points, one point, or 
two points or is a line all of whose points are on the surface. 

ILLUSTRATIVE EXAMPLE 

Find the points of intersection of the line j oining ( 1 , 2, 3) , (2, 1,4), 
and the surface x 2 7/ 2 + 2z 2 + 1 = 0. 

Solution 

We write x = - 1 + 3t, y = 2 - 3, z = 3 + t and have (-1 + 3) 2 
- (2 - 3<) 2 + 2(3 + O 2 + 1 = 9 2 - 6f + 1 - (9* 2 - 12$ + 4) + 
2(t* + 6t + 9) + 1 = 2* 2 + ISt + 16 = 0, * 2 + 9t + 8 = (t + 8)(t + 1) 
= 0. Hence, t = 8, 1, and the points are ( 25,26, 5), ( 4, 
5, 2). 

EXERCISE 

Find the points of intersection of the lines joining PI to P 2 and the 
surface/ = f(x t y,z) =0 in the following cases: 

(a) P l = (0, 0, 0), P 2 = (2, -1, 3),/ = 2*' + y - ^ + 1 
(6) Px = (0, 0, 0), P 2 = (2, - 1, 3), / s* 2* + y - 2* 

(c) Pi = (1, 1, 2), P,= (3,0,5), / = 2(*-l) + 2(-l)- 
(2 - 2)' - 1 

(d) Pi = (-1, 2, 3), P 2 = (1, 1, 2), / = s + t/ 2 - &' + 21 

(e) P! = (-4, 2, 3), P 2 = (-3, 3, 4), / =* ^ + y + - 2 
(/) P 1= (-1,2,3),P 2 = (l,l > 2),/ ss 3*-8j/-4' + 41 
(^) P! = (-1, 1,2),P 2 = (3, -4, l),/ = ^ + 2^ + 3^ -7 

10. Lines on a cylinder. The lines on a cylinder, consisting of 
lines parallel to the z axis and passing through a curve in the x, y 
plane, may be determined easily. We first prove the following: 

Lemma 1. Let f(x, y) be a homogeneous polynomial and. a, be 
real numbers not both zero. Then f(a, /3) = if and only if 
0x ay is a factor of f(x, y). 
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For if f(x, y) = (fix - ay)g(x, y), then /(a, ft) = (fia - aft). 
g(a y ft) = 0. Conversely, let /(a, ft) = 0, and define = /8x 
ay, 77 = aa; + fty. Since a 2 + /3 2 7* 0, we see by direct com- 
putation that 



X 



Then /(x, y) - 0(, ij), where 0(, r?) = a ^ + aig"-^ + 
+ a n rj n is necessarily homogeneous in and rj. When x = a and 
y = ft, we have = fta aft = 0, 77 = a 2 + /? 2 , and /(a, /?) = 
0(0, a 2 + 2 ) = a n (a 2 + ft*)* = 0. Then a n = 0, 0(, 77) SE fr 
lK, f) = /O, ?/) = (^ - 2/)00, y), where ^(x, y) = $(fix - ay, 
ax + fty). This proves the lemma. 

We use the lemma in the proof of the following theorem : 
Theorem 3. Let S be a cylinder defined by an irreducible poly- 
nomial equation f (x, y) = 0, of degree n > 1. Then there is only 
one line of S through each point of S. 

For if P is any point of /S, we may translate the origin to P 
and therefore assume that /(O, 0) = 0. Then f(x, y) = f n (x, y) 
+ / n -i(s, y) + - + fi(x, y) where each/fc(, y) is homogeneous 
of degree k. Every line L through the origin P has equations 
x = ta, y = tft, z = ty, and 

0(0 = /(te, 0) = */n(, 0) + f^Vn-iCa, /?)++ </i(, 0). 

Then L is on S if and only if 0(0 == 0, that is, / n (a, /S) = / n -i(, 0) 
=== /i(o:, /3) = 0. If (a, )8) ^ (0, 0), Lemma 1 implies 
that ftx ay is a factor of every /&(#, y) and thus is a factor of 
f(x, y) contrary to our hypothesis that n > 1 and that f(x, y) is 
irreducible. Hence, a = ft = and the only line through P is 

the z axis. 

EXERCISE 

Show that a surface f(x, y) = determined by a homogeneous poly- 
nomial f(x, y) of degree n consists of r ^ n planes through the z axis or 
is a point surface. 

11. Tangent lines and planes. A line x = XQ + t%, y = y Q 

+ ty, z = Zo + tf will intersect a surface /(x, y, 2) = in a point 
-Pi = (%i, yi, 21) if and only if there exists a real number ti which 
is a root of the equation 0(0 = f(x Q + , T/O + ft?, 2o + <f) = 
such that #1 = #o + *i, 2/i = 2/o+ <n?, i = ^o + /if. We shall 
say that the line is tangent to the surface at the point (of tan- 
gency) Pi if t\ is a root of multiplicity m ^ 2 of 0(0 = 0. 
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Let us investigate the problem of determining whether or not 
a given line through a point (z , 2/0, zo) on a surface f(x, y, z) = 
is tangent to the surface at the point. The point corresponds to 
the value t = and so we are assuming that <(0) = 0. We limit 
our discussion to functions f(x, y, z) having partial derivatives of 
all orders and see that <t>(t) has a Taylor series expansion 



which is an ordinary polynomial when/(x, ?/, z) is a polynomial. 
Then the given line is a tangent line if and only if <'(0) = 0. 

The function <'(/<) may be computed most easily by partial 
differentiation. Define f x (x, y, z) to be the partial derivative of 
f(x, y, z) with respect to x, that is, the derivative of f(x, y, z) as 
a function of x alone, and define f fj (x, y, z),f z (x, y, z) similarly. 
Then the derivative <f>'(i) of 4>(t) with respect to t is shown in 
elementary calculus to be given by the formula 

(9) *'(*) -f'^t+fy^t+f'Jt-te + fyl + te. 

Since x = rc , y = 2/o, z = z at t = 0, and (, 77, f) = (rr 2 , 2/2, 22) 
(#0, 2/o, z ), we see that a line joining the point (x , 2/0, z ) to a 
point P% (0:2, 2/2, ^2) is tangent to the surface at (#o, 2/o, Zo) if 
and only if P 2 is a solution of the equation 



(10) f x (x , y , ZQ)(X XQ) + f y (xo, 7/ 

2/0, ZQ)(Z 3 ) = 0. 



If the coefficients / x , / w , / 2 are all fcero at the point PO = (^o, 2/o, 
So), we call PO a singular point of the surface. All other points 
are called ordinary points. Since formula (10) is an equation of a 
plane at an ordinary point PO we have proved the following result : 

Theorem 3. A line passing through an ordinary point (x , yo, z ) 
of a surface f (x, y, z) =0 is tangent to the surface if and only if it 
is a line in the plane of formula (10). We call this plane the tangent 
plane to the surface at the given point. 

If PO is an ordinary point of two surfaces Si and $2, we define 
the angle between Si and $2 at PO to be the angle 6 between the 
two corresponding tangent planes at PO. Then Si and S 2 are 
said to be orthogonal at P if 6 = 90, i.e., the tangent planes are 
perpendicular. 
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The line perpendicular to the tangent plane of a surface S at 
an ordinary point PO of S is called the normal to S at PO. Equations 
of this line are 

(11) X - XQ = y - l/o 2 - Zp 



ILLUSTRATIVE EXAMPLES 

I. Find an equation of the plane tangent to the surface 3x 3 2y 2 + 
4z + 1 = at the point (-1, 1, 1). 

Solution 

We compute f x = 9#o 2 - 9, / = 4y = 4,/ z = 4 and the required 
equation is 9(x + 1) 4(y 1) + 4(z - 1) = Oz - 4?/ + 4z + 9 + 
4 4 = 0. 4ns. 9z + 4t/ - 4z + 9 = 0. 

II. Find equations of the normal line to the surface 3x 3 2y* + 4z + 
1 = Oat the point (-1,1, 1). 

Solution 
By our theory and the result of Example) I the answer is 

aHM y - 1 z - 1 



III. The point ( 1, 1, 1) is a point of intersection of the surface 
3z 3 - 2y 2 + 4z + 1 = and the surface 2x 2 - 4xy + z 3 = 7. Find 
the cosine of the angle 8 between the two surfaces at this point. 

Solution 

The plane tangent to the second surface is (4xo tyo)(x XQ) 
4#o(2/ l/o) + 3z 2 ( 20) = and therefore a set of direction numbers 
of the normal line is ( 8, 4, 3). A set of direction numbers of the 
first line is (9, 4, 4) and therefore 

| -72 - 16 + 12| _ 76 

COS v 



8)*+ (-4) 2 + 3 2 V& + (-4) 2 + (4) 2 

EXERCISES 

1. A homogeneous algebraic equation f(x, y, z) =0 of degree n ^ 2 
defines a cone with vertex at the origin. Show that the vertex is a 
singular point of the cone. 

2. Find an equation of the tangent plane to each of the following sur- 
faces at the given points: 
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(a) x 3 - 2y* + 5z 2 = 18, (0, - 1, 2) 

(6) 3z 2 + 47/ + 2 2 = 0, (0, -1,2) 

(c) x*-x*y + z* = 3, (1, -1, -1) 

(d) 2xy-z* = -I, (!,-!,- 1) 

(e) a;?/ - 82 = 3, (2, 3, 1) 
(/) xyz = 6, (2, 3, 1) 

(g) x*z + V = 2, (-1,0,2) 

(A) zt/ 2 + yz 2 + **' = 2, (-1,0,2) 

3. The problems in Exercise 2 may be grouped in pairs such that the 
given point P is a point on a pair of planes. Find the cosine of the 
angle between the planes for each pair. 

12. Tangents to quadrics. A quadric surface is defined by a 
polynomial equation 

(12) f(x, y, z) = ax* + by 2 + cz* + 2(dxy + exz + gyz) 

+ 2(hx + py + qz) + s = Q, 

where the coefficients a,b, . . . ,s are real numbers and a, 6, c, 
rf, e, g are not all zero. Then the tangent plane to this surface at 

(X Q , 7/o, 2 ) 1 s ( ax <> + dy<> + ez <> + ty (X XQ) + (6?/o + dx Q + QZ Q 

2/o) + (CZ Q + ex Q + gy Q + q)(z ZQ) = ax Q x + by Q y + 
d(xy Q + Xoy} + e(xz Q + x z) + g(yz Q + z y) + h(x + XQ) 
+ 2/o) + q(z + ZQ) + s [ax 2 + by<? + cz 2 + 2(dx Q y Q + 
+ gyozo) + 2(hx Q + py Q + qz ) + s]. Since (X Q , T/ O , 2 ) is 
on the quadric, this equation reduces to 



(13) axx Q + 67/1/0 + CZZQ + d(xy + x Q y) + e(xz 

o + yoz) + h(x + XQ) + p(y + T/O) + q(z + 2 ) + s = 0. 



This formula arises from the original formula (12) of a quadric 
surface by the replacement of square terms such as a: 2 by corre- 
sponding products such as XX Q , terms 2xy by XT/O + 2/#o, and terms 
2x by x + XQ. It may then be easily remembered. 

ILLUSTRATIVE EXAMPLE 

Find an equation of the tangent plane to the quadric surface 3x* + 4y* 
-2z 2 + 6xy - 3xz + 6z - 3z - 7 = at the point ( 1, 3, 2). 



Solution 

tie equation is 3x + \2y - 

- f (z + 2) - 7 33 6z + 9y - 4z - 13 = 0. 



The equation is -3z + \2y - 4z + 3(3x - y) - |(2x - z) + 3(z - 1) 



See. 12] SURFACES AND CURVES 53 

EXERCISE 

Write an equation of the tangent plane to each of the following 
quadrics at the given points: 

(a) x 2 - 2y* + W = 18, (0, -1, 2) 
(6) 3z* + 4y* + 2z = 0, (0, 1, -2) 

(c) x*-xy + z* = Z,(l, -I, -1) 

(d) 2xy - 3x2 + 52 = 2, (1, -1, 2) 

(e) 2x* -y* + z* + 2xy + xz + x = Q, (0, -1, 1) 



CHAPTER 4 
SPHERES 

1. Equations of spheres. Let r be a positive real number and 
PO = (#o, 2/o, Zo) be a fixed point. Then we define a sphere S 
of radius r and center at P to be the locus of all points P whose 
distance from P is r. It follows that a point P is on S if and 
only if \/(x z ) 2 + (y yo} 2 + (z z ) 2 = r, that is, if and 
only if P is a solution of the equation 

(1) (x - x.r + (y - 2/o) 2 + (z - z ) 2 = r 2 . 

We have then succeeded in deriving an equation of an arbitrary 
sphere. Evidently a sphere is a quadric surface. 

Theorem 1. Let a, b, c, d, e be real numbers such that a 5^ 
and put 

=Xri a ' a ' p a 2 

Then the equation 

(3) f (x, y, z) = a(x 2 + y 2 + z 2 ) + 2(bx + cy +'dz) + e = 

has no real points as solutions ifp<0 and is called an equation of 
an imaginary sphere. If p = 0, the only solution of f (x, y, z) = 
is the point PO = (XQ, yo, Zo) and the equation is said to be an equa- 
tion of a point sphere. The only remaining case is the case p > 
and f (x, y, z) =0 is an equation of the sphere whose center is P 
and whose radius is VP- 

For a point is a solution of f(x y y, z} = if and only if it is a 
solution of (\/a)f(x, y, z} = 0, that is, of 

2 i 7 2 2 . ?6 , 2c 2d ^ -e 
a a a a 

This equation is equivalent, in view of the definitions of formula 
(2), to 

(4) (x - z ) 2 + (y - 2/ ) 2 + (z - z ) 2 = p. 

There are no real solutions (x, y, z} of this equation if p < 0, 

54 



Sec. 2] SPHERES 55 

its only real solution is (X Q , y , ZQ) if p = 0, and it reduces to 
formula (1) with r 2 = p if p > 0. 

EXERCISES 

1. Write an equation of a sphere with center at Po and radius r in the 
following cases: 

(a) Po = (0, 0, 0), r = 5 (d) PO = (1, 6, 2), r = 6 

(6) Po = (-1, 2, 0), r = 3 (e) P = (1, 2, 2), r = 3 

(c) Po = (-1, -1, 3), r = 7 (/) Po = (-2, 6, 3), r = 7 

2. Determine the center and radius of the following spheres: 
(a) x 1 + 7/ 2 + z 2 + 6z + 4y + Sz + 25 = 

(6) z 2 + y 2 + z* + 2x - Qy + 4* + 10 = 

(c) 2(x 2 + 7/ + 2 2 ) + 2x - 3y + 2^ - 1 

(d) x 2 + 7/ 2 + 2 2 - 3x + 5y - 7z = 
(6) 3(a; 2 + ?/ + 2 2 ) + x + y + z - 6 

3. Show that the lines passing through the center of a sphere are 
normal lines to the sphere. 

2. Spheres satisfying given conditions. The equation of 
formula (3) is a linear homogeneous function of the five param- 
eters (variables not regarded as point coordinates) a, 6, c, d, e. 
It may then be made to satisfy four conditions such as passing 
through four points. 

ILLUSTRATIVE EXAMPLES 

I. Find an equation of the sphere with center at ( 1, 2, 3) and pass- 
ing through the point (1, 1, 2). 

Solution 

The restriction on the center is a set of three conditions. We use 
the formula (1) and write (x + I) 2 + (y - 2) 2 + (z - 3) 2 = r 2 . Then 

(1 + I) 2 + (1 - 2) 2 + (2 - 3) 2 = 4 + 9 + 1 = 14 = r 2 . 

Ans. (x + I) 2 + (y - 2) 2 + (z - 3) 2 = 14. 

II. Find an equation of the sphere through the points (1, 1,0), 
(2,1,1), (3, -1,4), (-1, -1,2). 

Solution 

We substitute the coordinates of the four given points to obtain the 
equations 

2a + 2(6 - c) + e = 

6a + 2(26 + c + d)+e = 

26a + 2(36 - c + 4d) + e = 

6a + 2(-6 -c + 2d) +e - 
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Then e = 2(a + b e), and we subtract the first of our equations 
from the remaining equations and divide by 2 to obtain 

2a+ b + 2c+ d = 

12a + 26 + 4d = 

2a - 26 + 2d = 

This gives 2c = (2a + 6 + d), and we need to solve 

6a + 6 - -2d 
a - 6 = - d 

The solution is 7a = 3d, and we avoid fractions by selecting d = 14, 
a = 6, 6 = a + d = -8, 2c = -(12 - 8 - 14) = 10, c = 5, e = 
-2(6-8-5) = 14. Thus an equation is 6(o; 2 + 7/ 2 + z 2 ) + 2(-8z 
+ 5y - 142) + 14 = 0. 

Ans. 3(z 2 + y 2 + z*) - 80; + 5y - 14z + 7 = 0. 

EXERCISES 

1. Find an equation of a sphere with center at P and passing through 
Pi in the following cases: 

(a) P = (0,0,0),P! = (-1,2,3) 

(b) P = (-1,1,0),P! = (1, -1,2) 

(c) Po= (1, -l f 3),P!= (2,3, -4) 

(d) Po = (-1, -3,2),^ = (-2, 1, -4) 

(e) P = (0, -1,0),P!= (4,0, -3) 
(/) P = (1,2,3),P 1 = (-1,0,2) 

2. Find an equation of a sphere passing through the points Po, Pi of 
Exercise 1 and the points P%, PS as follows: 

(a) P 2 = (1, -2, 1), P 3 = (1, 1, 2) 
(6) P 2 = (1,2, -1),P,= (0,0,1) 

(c) P 2 = (2, 1, 1), P 3 = (1, 0, -3) 

(d) P 2 = (1, -1,2),P 3 = (2,3, -4) 
(e)P 2 = (-2, -1, 1),P,= (-1,1,1) 
(/) P 2 = (3, 2, 1), Pi = (2, 2, -1) 

3. Linear families of spheres. Consider the two spheres Si 
and Sz denned by the equations 

(5) f(x, y, z) s ai(z 2 + y 2 + 2 2 ) + 2(6^ + Cl y + d^z) + ei = 0, 

(6) g(x, y, z) = a 2 (x 2 + y z + z 2 ) + 2(6 2 a; + c 2 y + d^) + e 2 = 0. 
Then the equation 

(7) s/(x, y, 0) + ty(z, y, ) = 

ia an equation of the form of formula (3) in which 

(8) a = sai + ta z , b = sbi + tbz, c = sci + tc 2 , 

d = sdi + td%, e = sei -\- tei. 
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It follows that if 

(9) a = soi + to 2 ^ 

the equation of formula (7) is an equation of a (point, imaginary 
or ordinary) sphere. We shall therefore call the family of sur- 
faces defined by formula (7), as s and t range over all real numbers, 
a family of spheres generated by the two given spheres. 

If Po = (XQ, 7/0, z ) 1 s an Y point on both f(x, y, z) = and 
g(x, y, z) = 0, then P is also on the surface sf + tg = 0. Thus 
all spheres of the family generated by the two given spheres con- 
tain all of points of intersection of these two spheres. Conversely, 
if S is a sphere through the curve C of intersection of two spheres 
Si and /S 2 , and S contains a point PO not on C, then S is the sphere 
S(s, t) of the family defined by formula (7) when we take 5 = 
0(^0, ?/o, Zo) and t = f(x Q , yo, z ). For there is at most one 
sphere through the curve C and a point not on C. The solution 
may, of course, be an imaginary sphere or a point sphere. 

The values of s and t for which a = sai + ta% yield the 
member of the family of spheres whose equation is 

(10) 2o 2 (6iz + ciy + d&) - 2ai(b*x + c z y + d z z) 

+ (0261 Oie 2 ) = 0. 

This is the equation of a plane except when 02(61, Ci, d\) = 
01(62, 2, ^2). The plane is called the radical plane of the two 
spheres. It contains the curve C of intersection of the two 
spheres and C is then the curve of intersection of the radical 
plane with either sphere. It follows that if two spheres intersect 
in at least two distinct points they intersect in a circle in their 
radical plane. 

EXERCISES 

1. Show that the line joining the centers of two intersecting spheres 
cuts their radical plane in the center of the circle of intersection. 

2. Find the center and radius of the circle of intersection of the sphere 
x 2 + y 2 + z 2 = 4 and the sphere (x - I) 2 + y* + z 2 = 9. 

3. Find an equation of the sphere through the circle of intersection 
of the spheres of Exercise 2 and the following points determining the 
center and radius of the solution sphere in all cases where it is not a 
point or an imaginary sphere. 

(a) (-1,2,3) (e) (-1,1, -3) 

(b) (4, 0, 0) (/) (1, 1, 1) 

(c) (-3,0,0) (g) (-1.-1.0) 

(d) (1, -1,2) (A) (2, -1,3) 



58 SOLID ANALYTIC GEOMETRY [Chap.5 

4. Angles between spheres. The tangent planes to the 
spheres of formulas (5) and (6) at a point of intersection (X Q , yo,z ) 
are given by 



(11) ai(xx + yy Q + zz ) + bi(x + XQ) 

+ di(z + z ) + ei = 0, 

(12) a^(xx Q + yyo + zz ) + b^(x + XQ) + c z (y + yo) 

+ d 2 (z + z ) + <? 2 = 0. 

Then the corresponding sets of direction numbers are (a\x Q + 61, 
aiy Q + ci, diz + di) and (0,2X0 + 6 2 , a 2 yo + c 2 , a 2 z + d 2 ). The 
square of the length of the first of these vectors is 



+ z 2 ) + 2 

= bi 2 + cS + dS - 

since (XQ, y Q , z ) is a point of fi(x, y, z) = 0. Similarly the square 
of the length of the second vector is &2 2 + C2 2 + d^ 
The inner product of the two vectors is 



g = (aizo + bi)(a 2 x Q + 62) 

+ 2/o 2 + z 2 ) 



6162 + CiC 2 + did 2 = ia 2 [ai(z 2 + 2/o 2 + z 2 ) 

+ i/o 2 + z 2 ) 



Then g = (a 2 ei + i^) + 6162 + ^ic 2 + eW 2 and the cosine 
of the angle between two planes at any point of their circle of 
intersection is given by the formula 

(13) 

_ |2(6i6 2 + CiC 2 + did*) 

COS \j - ..... ' 



Since this formula does not involve the coordinates X Q , y , z , we 
have proved the following result : 

Theorem 2. Two spheres intersect at the same angle at all points 
of their circle of intersection. 



CHAPTER 5 
QUADRIC SURFACES 

1. Ellipsoids. In this chapter we shall discuss the geometric 
properties of the quadric surfaces given by a special set of equa- 
tions. This discussion will yield the properties of all quadrics, 
since we shall prove in Chap. 7 that any equation of a quadric can 
be carried into one of the equations we shall discuss, by a suitable 
choice of coordinate axes. 

If a, b, c are any positive real numbers, the equation 

~2 ? .2 2 2 

m ?L -I. *L -u L 1 

(L) a 2 + 6 2 + c 2 

is an equation of a surface called an ellipsoid. Such a surface; 
is symmetrical with respect to the origin, thg coordinate planes,, 
and the coordinate axes, since it is a function of x 2 , ?/ 2 , and z 2 . 

The six points (a, 0, 0), (0, b, 0), (0, 0, c) are called the 
vertices of the ellipsoid. They are the points of intersection of 
the coordinate axes and the surface. The three line segments 
joining the pairs of vertices on each coordinate axis are called 
the axes of the ellipsoid. The axes intersect at the origin of 
coordinates, and we call this point the center of the ellipsoid. The 
line segments joining the center to the vertices have lengths 
a, b, c and these lengths are called the semiaxes of the ellipsoid. If 
we order these three positive real numbers 2a, 26, 2c, the largest 
of the corresponding segments is called the major axis of the 
ellipsoid, the next largest the mean axis, and the smallest the 
minor axis. 

We shall show in Chap. 7 that every plane section of a quadric 
is a conic. In particular we shall show that the plane sections 
of an ellipsoid are ellipses. The plane sections by the planes 
z = k are ellipses 

* y _ _ * _ c-* 

(2) + i - 1 -- ~ 



provided that c 2 > k 2 . If c 2 = k 2 , the corresponding planes 
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z = c are tangent to the ellipsoid, and there are no real points 
of intersection if k > c. The student should make a similar 
analysis of the plane sections by planes x = k and y = k and 
should examine Fig. 11 carefully. 

z 



(o,o,c) 




FIG. 11. 

If two of the numbers a, 6, c are equal, the ellipsoid is a surface 
of revolution about that axis which corresponds to the remaining 
letter; for example, the ellipsoid 



(3) 



= 1 



is an ellipsoid obtained by revolving the ellipse 

rr 2 z 2 
(4) 73 + z - t - 1, y = 



about the 2 axis. If a > c, this ellipsoid of revolution is the 
result of a revolution of an ellipse about its minor axis and is 
called an oblate spheroid. If a < c, the ellipsoid is the result of a 
revolution about its major axis and is called a prolate spheroid. 
If a = c, the ellipsoid of revolution is a sphere. 
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If a T* 0, /8 9* 0, y 7* 0, d, e, g, h are real numbers, the equation 
(5) f(x, y, 2) = ax 2 + 0y* + yz* + (dx + ey + gz) + h = 
is equivalent to 



where 

d 2 e 2 n 2 

m ,_- + |. + + . 

Suppose that a, 0, 7 all have the same sign. Then if p = 0, the 
only real point satisfying f(x, y, z) = is the point (x , 2/0, 2o) 
defined by 

(8) x = - ^, t/o =- ^- 20 = - , 

and the equation is called an equation of a point ellipsoid. If p 
and a. have opposite signs, there are no real solutions of f(x, y, z) 
= and this is an equation of what is called an imaginary 
ellipsoid. If a and p have the same signs, the equation /(x, y, z) 
= is clearly an equation of the ellipsoid whose semiaxes are 



(9) 



/P /P 

\' \0' 



whose axes are on the lines of intersection of the planes x = X Q , 

y = 2/0, z = 2 , and whose center is the point (0*0, 2/o, 2 ). 

EXERCISES 

1. Give the coordinate of the center, the lengths of the semiaxes, and 
the equations of the lines on which the axes lie for the following ellipsoids: 

(a) 9x 2 + 367/ 2 + 4z 2 - 18x + 144y +117 = 
(6) x 2 + 2i/ 2 + 3z 2 -2a? + 8y-fts + ll=0 

(c) 3rc 2 + 4i/ 2 + 52 2 + 6z - 16j/ + 10s + 23 = 

(d) 3a; 2 + 3?^ 2 + 2z 2 - 6 - 12y - 12 + 29 = 

2. Give an equation of an ellipsoid with axes parallel to the coordinate 
axes, major axis 2a parallel to the y axis, minor axis 2c parallel to the 
z axis, and center at P in the following cases: 

(a) P = (-l,2,3),a = 4,6 = 3,c= 1 
(6) P (1, -2, 0), a = 5, 6 = V% c = 2 

(c) P = (-1, -3, 1), a = 2, 6 = \/3, c = \/2 

(d) P = (4, -1, 2), a = 3, 6 = 2, c = V~2 
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3. Give the equations of Exercise 2 if the major axis is parallel to the 
z axis and the minor axis to the x axis. 

2. Quadric cones. Every set of positive real numbers a, 6, c 
defines an equation 



which is homogeneous in x, y, z and thus is an equation of a 
quadric cone with vertex at the origin. It is evidently a surface 
symmetric with respect to the origin, the coordinate axes, and the 
coordinate planes. 

The nature of such a surface is indicated by its sections with 
planes z = k ^ 0. These sections are the curves 

(ID * = *, x i + l = 

and thus are ellipses with semiaxes 

ah bh 

- ) - ; 

C C 

where h = \k\. The student should analyze the sections made 
by the planes x = k and y = k. 

Every plane through the origin and not parallel to the x axis 
is a plane x = ay + )8z, where a and ft are real. The intersection 
of such a plane with the cone consists of points satisfying the 
equation x = ay + /3z and the equation 

fo )2 4- v* z * - n 
> - + - - 0. 



The coefficient of y 2 is a positive real number, and we may desig- 
nate it by 

(14} T 2 = 4- - 

u * ; T a 2 ^ 6 2 

The coefficient of 2/2 may then be designated by 275 and that of 
z 2 by 5 2 p. Then formula (13) becomes 

(15) (yy + Sz) 2 - pz 2 = 0. 

If p < 0, the only points on the given plane and cone are the 
points z = yy + 82 = 0, x = ay + ftz, and thus yy = 0, y = 0, 
x = 0. Hence, when p < 0, the only point of intersection of the 
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plane and the cone is the origin. When p ^ 0, the intersection 
consists of the line 

(16) x = ay + pz, yy + (d + Vp)z = 
and the line 

(17) x = ay + fr, yy + (8 - VP)Z = 0. 

These lines coincide if p 0. 

Every plane through the origin is either the plane z = which 
cuts the cone of formula (10) in the origin or is a plane which is 
surely not parallel either to the x axis or to the y axis. Since x 
and y have symmetric roles in formula (10) if a plane is parallel 
to the x axis and therefore not parallel to the y axis, we may 
repeat the argument above with the roles of x and y interchanged. 
We have therefore proved the following theorem: 

Theorem 1. The intersection of the quadric cone of formula (10) 
with a plane through its vertex consists either of the vertex or of two 
lines that may be coincident. 

3. Hyperboloids. If the numbers a, /3, 7 of formula (5) are not 
zero but do not all have the same sign, we may multiply by 1 
if necessary and hence assume that 

(18) a > 0, j8 > 0, 7 < 0. 

We may then convert formula (5) into formula (W and see that 
if p = formula (6) becomes the equation 

na\ ( x "" *)* -j- ^ ~~ y) 2 (* - *) 2 n 

(19) - -f + - p -- *~~ = 

of a quadric cone with 

1 a l l 

- -,' ft = F , -y = -, 

and vertex (x , yo, Zo) given by formula (8) . Assume then that 
P 7* 0. 

If p > 0, formula (6) becomes the equation 



where 

(21) a* =?, 1* = ?, c = - ^- 

a 7 
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This is an equation of a surface called a hyperboloid of one sheet. 
The point XQ, yo, ZQ is called the center of this surface, and the 
surface is symmetric with respect to the translated origin, the 
translated coordinate axes, and the translated coordinate planes 




x 



FIG. 12. Hyperboloid of one sheet. 

after a translation of axes that carries the origin to the center. 
Let us assume then that the origin is the center and therefore 
study the equation 

(22) % + jg - fj = 1. 



(23) 



The sections of this hyperboloid by the plane z = k are ellipses 

2 ,,,2 1.2 _|_ C 2 



x | ?/ 2 _ 

To I TV 



z = k 



with center at the origin for all real values of k. The semiaxes of 

these ellipses are 

n h 

Uc / V"^ : n t/ 



(24) 



and increase as \k\ increases. 

The sections of the hyperboloid of formula (22) by planes 
y = k are hyperbolas 
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for all values of k except k = 6, and the semiaxes of these 
hyperbolas are 

c 



(26) 



- H 



When k = 6, the corresponding sections are pairs of straight 
lines 



(27) 



z = - 
a 



y = k. 



A similar analysis of the plane sections by planes x = fc should be 
carried out by the student. 

When p < 0, formula (6) becomes 



(28) 

where 

(29) 



(x- 



I 



(y - y,Y _ (z - 



a 2 = - 



This surface is called a hyperboloid of two sheets. The point 




c* a* 
FIG. 13. Hyperboloid of two sheets. 
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(#o, yo, Zo) is called the center of the hyperboloid, and we may 
translate the axes to the center and thus simplify the equation. 
Let us assume then that the origin is the center and therefore 
study the equation 

? 2 r 2 i/ 2 

< 30 > --- ' 



We shall call the numbers a, b, c of formulas (19), (20), and (28) 
the semiaxes of the three types of surfaces we are considering. 

The sections of a hyperboloid of two sheets by the planes 
z = k are given by 



t/ 2 _ /c 2 - c 2 
----- ~~~ 



and therefore are ellipses with semiaxes 
(32) 




provided that \k\ > c. If, however, c < k < c, the plane 
z = k does not intersect the hyperboloid and thus there are no 
points on the hyperboloid of formula (30) between the plane 
z = c and the plane z = c. The points (0, 0, c) and (0, 0, c) 
are the respective intersections of the planes z = c and z = c 
with the hyperboloid and are called its vertices. It should be 
noted that a hyperboloid of two sheets actually consists of two 
separated surfaces such, that z ^ c on one of the parts and 
z ^ c on the other. 

The sections of a hyperboloid of two sheets by planes x = k 
are all hyperbolas 



with semiaxes 



+ a 2 b Vk* + a\ 



a a 

Similarly the sections by planes y = k are all hyperbolas. 

We have now completely classified all surfaces given by equa- 
tions of formula (5) for apy ^ 0. 

EXERCISES 

1. Give a necessary and sufficient condition that the surfaces of 
formulas (19), (20), and (28) shall be surfaces of revolution. 
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2. Classify the following surfaces and give the center, the axes of 
symmetry, and the semiaxes in each case : 
(a) x* + 2ij* - z* + 6z - y + 11 = 
(6) x* - 2?/ 2 - 3z 2 + 4z - 6y + 12z + 11 = 

(c) 2x* - 3?y 2 + 40 2 + 6x - OT/ + 120-0 

(d) 2x* - 27/ 2 - 3* 2 + 16?/ + 12^ = 44 

(e) 4* 2 + 7/ 2 - 9z 2 + 8x + 2y - 18z + 11 = 
(/) 4z 2 - 4i/ 2 - 9^ 2 - &c + 8y - 862 = 

4. Lines on a hyperboloid. We shall begin our discussion of 
lines on a hyperboloid by proving the following: 

Theorem 2. A hyperboloid of two sheets contains no lines. 

For the line x = XQ + t, y = T/ O + to?, 2 = z + # will lie on 
the hyperboloid 



if and only if 

(*o + ^) 2 _ (2/c + ^?) 2 , (go 



for all values of t. This equation is then an identity in t and the 
coefficients of 1, t, and i 2 in the two members must be equal, i.e., 

*o 2 _ yo 2 , Zo 2 , a:o^ _ /o? , 2of f 2 _ 7 2 , f 2 

-- 



Then 



and 

2g o 2 + 



1 |_ i_ 

o i^ 9 i /r \> Vf - 

6 J c 2 (be) 

However, this expression is a sum of three real squares and can 
vanish only when the squares vanish separately. Thus rj = f 
= 0, and the last equation of formula (36) implies that 2 = 0, 
= 0. Then our line degenerates to x = #o, y = 2/o, 2 = 2 , that 
is, no line lies wholly on the surface. 

We next pass on to the case of a hyperboloid of one sheet. 
The lines 

(39) - + - = 
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defined for all real numbers X, consist of points wholly on the sur- 
face defined by formula (22). For if P = (XQ, 2/o, Zo) * s a point 
on a line of formula (39), then 



and Po is on the hyperboloid. Conversely, if P is on the surface, 
then 



If the denominator of 
(42) X = 



(ZQ/C) 



is not zero, the point Po lies on the unique line of the family of 
formula (39) determined by this value of X. If both numerator 
and denominator are zero in formula (42) P is on the line uniquely 
determined by X = xo/a. Finally, if the numerator of formula 
(42) is not zero and the denominator is zero, then P is on the line 

^ b ~~ ' a c ~~ 

This line is regarded as being that member of the family of lines 
of formula (39) defined for the infinite value of X. 

The family of lines defined by formula (39) is called a regulus 
of the corresponding surface, and we have proved that through 
every point of the surface there passes one and only one line of 
the regulus. The family of lines defined by 

is a second regulus of the surface of formula (22), and it should 
be evident that it can be similarly proved that one and only one 
line of this second regulus passes through a point on the surface. 
The two reguli of a hyperboloid have no lines in common; for 
let us suppose that a particular value of X and a particular value 
of M give the same line. The point y = 6, x = aX, z = c\ is on 
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the line of formula (39) and when substituted in formula (44) 
gives /i = 0. Then the point y = 6, x = z = is on this second 
line and substitution in formula (39) yields A = 0. However, 
the two resulting lines 

(45) l=-- c > y = b; l = - c > y = b 

have only the point (0, 6, 0) in common and cannot coincide. 

If P is a point on a hyperboloid of one sheet, there are exactly 
two lines on the surface that pass through P. For the line 
through XQ, 7/0, Zo an d having direction numbers , rj, f will lie on 
the surface 8 of formula (22) if and only if 

(so + <) 2 , (yo + fr) 2 _ (so + #) 2 

a 2 "*" 6 2 c 2 

is identically zero in t. This requires that (XQ, y Q , z ) shall be a 
point of iS and that 



s I */___ __^ S r\ _-^ - _L ^^* * n 

tf "*" 62 ^2 - u > ^T f -Jj 72 - 

Then (|, 17, f) is a point not the origin and on the plane 



_ __ 

O I 19 9 

a 2 6 2 c 2 



and the cone 



By Theorem 1 this plane through the origin cuts the cone in 
exactly two lines through the origin. Let the coordinate of two 
nonzero vectors on these two lines be (1, ijj, f i) and ( 2 , 172, ^2), 
i-espectively, so that every point on the first line is a vector 
(> I, f) = p(i> ^i> f i) an( i every point on the second line is a 
vector (, r;, f) = #(2, ^72, ^2). Since all nonzero scalar multiples 
of a vector (, ry, f ) define a set of direction numbers of the same 
line as (,*?,"), we have shown that there are not more than 
two lines through a point on the surface of a hyperboloid. We 
have also proved the existence of two lines lying in two reguli and 
the proof is complete. 

EXERCISES 

1. Substitute the value y = 6 in formulas (39) and (44), and thus 
determine a point on a corresponding line of the reguli. Determine a 
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second point with y = 6, when X 7* and when ju 7* 0. However, 
when X = 0, show that (a, 6, c) is a point of the line and that, when 
jut = 0, (a, 6, c) is a point of the line. Use these values and the infinite 
case to determine a set of direction numbers for any line of the two reguli. 
2. By symmetry we see that a hyperboloid of formula (22) has the 
regulus 



and the regulus 

y z /i i x \ (y r z \ 1 * 
_ = ^ 1 1 -^ i, ~i_.j__i == i 

6 c ' \ aj ' \b cj a 

By substituting the points of Exercise 1 show that the first of these two 
reguli is precisely the regulus of formula (44) and that the second is the 
regulus of formula (39). 

3. Show that every plane through a line of the regulus of formula (39) 
also contains a line of the regulus of formula (44). HINT: The planes 
through a line of formula (39) are all of the form 



4. Show that every plane of Exercise 3 is tangent to the surface at the 
point of intersection of the two generators. 

5. Write out the equations of the reguli of the following hyperboloids, 
and determine which pair of lines of the reguli pass through the corre- 
sponding point. 

r 2 7/ 2 z 2 

(a) T + ? -le- 1 '* 2 '- 1 '** 

W ? + 2/ 2 -^ = i, (3,+, -i) 

(c) ^^- 2 -*' = !, (3 V2, 0, 1) 

6. Show that there are no lines on the surface of an ellipsoid by using 
the method of the proof of Theorem 2. 

6. Paraboloids. By consideration of symmetry we see that a 
discussion of formula (5) in the case where two of the numbers 
a, /3, 7 are not zero and the remaining number is zero need only 
be carried out for the case where a/3 ^ 0, 7 = 0. We may also 
assume that a > 0, since we may multiply formula (5) by 1 if 
necessary. Let us also assume in this section that g ^ 0. We 
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may then convert formula (5) into 

(46) a(x - a; ) 2 + j8(y - y ) 2 + 2g(z - z ) = 0, 
where 

(47) xo = > Vo = > ZQ = 



71 



h . d* 



We translate the axes to (#0, 2/o, z) and divide the equation of 
formula (46) by \2g\ to convert it to the form 



(48) 

if a/3 > 0, and to 

(49) 

u, 

if a/3 < 0, where e = 1 or 1 . 



x^ , y 2 _ 

a 2 "I" 6 2 ~ *> 



x 2 __ y 2 __ 




FIG. 14. Elliptic paraboloid. 

Equation (48) is an equation of a surface called an elliptic 
paraboloid. Its plane sections, by planes z = fc, are ellipses 



(50) 
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if cfc > 0. The plane z = cuts this parabola at the origin and 
this point is called the vertex of the paraboloid. The plane sec- 
tions x = k and y = k of the paraboloid are parabolas. 



(51) 
and 

(52) 



Js iv j 

? - - 55' - * 



= v X = k. 




FIG. 15. Hyperbolic paraboloid. 

Equation (49) is an equation of a surface called a hyperbolic 
paraboloid. Its sections, by planes z = k, are all hyperbolas 



(53) 

except for the section 






- - - ek 

70 t*^ 



by the plane 2 = 0. This section consists of the line 

z = o 



= 

a 6 



and the line 



a 



= 0. 
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The hyperbolic paraboloid of formula (49) has two reguli. 
They are given by the equations 



and 

(66) !_J.^ 

The student should show, as an exercise, that the lines of the first 
regulus are parallel to the plane bx = ay and that those of the 
second regulus are parallel to the plane bx = ay. These are 
then two distinct lines through each point of the surface. The 
hyperboloid of one sheet and the hyperbolic paraboloid are some- 
times called ruled surfaces. 

It is easy to show that there are no lines on an elliptic parabo- 
loid; for a line is defined by direction numbers (f, 77, f) j (0, 0, 0) 
and a point (z , 2/0, ZQ) such that 



f % JL ^ ^- 
I"9~i"r9/"~ s = 

\a 2 6 2 / 



_ \** / U ' V S/ I \f7 

and thus 

Ji + 52 = > 2 

Then = T; = f = 0, a contradiction. 

We have already found two distinct lines through each point 
(#o, 2/o, Zo) f the surface (49), and they are the line of formula 
(54) determined by 

(x Q /a) + (yo/b) 1 



and the line of formula (55) determined by 
(x/a) - (y/b) 1 



M = 



Z Q (xo/a) + (2/0/6) 



To prove that these are the only two lines, we examine the 
equation 



74 SOLID ANALYTIC GEOMETRY [Chap. 5 

and thus see that 



tf fro Wo _ n 

55' -jr - &r - <f ~ - 



Then 



(o_ 
> 
fl 



- 

eab 



i.e., all sets of direction numbers , 17, f are scalar multiples of the 
two sets 

b x Q 3/0 \ /- b x Q y Q \ 

,-; r -- -I; II, -- ; - 5 -\ -- r I- 

a ecr ab/ \ a eer eab/ 

This proves that there are exactly two lines on the surface through 
each point of a hyperbolic paraboloid. 

EXERCISES 

1. Determine a set of direction numbers for the lines of formulas (54) 
and (55) by using the value z = and the value y = kb. 

2. Write out the equations of the reguli of the following paraboloids 
and of the particular lines of the reguli through the corresponding 
points: 

() - = 22, (4, 3, f) 
(6) ^-^= -&, (2,1,0) 

< c > f ~ 1 = z > < 6 > 6 > ~ 5 > 

6. Cylinders. The equation of formula (5) reduces to 



if a/3 j and g = 0. If a > 0, ft > 0, this equation is an equa- 
tion of what is called an elliptic cylinder. Its sections by planes 
z = k are evidently ellipses, and we have already proved that 
there is one line on the surface through each point of it. 

If a > 0, /3 > 0, and p = 0, the surface defined by formula (56) 
consists only of the point 

d d 

X= ~2~a' y== ~2(3' 

and if aft > 0, p < 0, there are no real points on the surface. 
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If a > and j8 < 0, formula (56) is an equation of a hyper- 
bolic cylinder when p = and of two planes through the z axis 
when p = 0. This completes the case of surfaces defined by 
formula (56). 

There remains the case where two of the numbers a, 0, 7 of 
formula (5) are zero. We may then assume that a 1, since 
we can always divide the equation by a 7* 0. Then formula (5) 
becomes 

(57) z 2 + 2(dx + ey + gz) + h = 0. 

In Chap. 7, we shall show how to select a coordinate system 
such that a surface defined by formula (57) has the form of (57) 
with g = 0. Then the equation reduces to 

(58) (x + d)* = -2ey + d* - h. 

This is an equation of a parabolic cylinder when e j 0; it is an 
equation of two distinct parallel planes when e 0, d 2 ? /i; it 
is an equation x = d of two coincident planes when 6 = 0, 
d 2 = ft; and is a surface with no real points when e = 0, d < h. 

7. Classification of quadric surfaces. A quadric surface is 
defined by an equation f(x, y, z) = 0, where f(x, y, z) is a poly- 
nomial of degree two. We shall prove in Chap. 7 that it is 
always possible to select a coordinate system such that the sur- 
face is defined by an equation of formula (5) where the numbers 
a, )8, 7 are independent of the choice of the coordinate system. 
Let us call these numbers a set of characteristic roots of the sur- 
face. They are unique apart from a proportionality factor t, 
which may be introduced by multiplying /(x, y, z} = by t ^ 0. 

Quadric surfaces may be classified according to the properties 
of their sets of characteristic roots and the number of lines of the 
surface through each point of the surface. The classification 
follows: 

1. CHARACTERISTIC ROOTS ALL OF THE SAME SIGN. Such 
surfaces are either ellipsoids, point ellipsoids, or imaginary ellip- 
soids. There are no lines on such surfaces. 

2. THREE NONZERO CHARACTERISTIC ROOTS NOT ALL OF THE 
SAME SIGN. 

a. The quadric cone. This is a surface containing a point (its 
vertex) through which pass infinitely many lines of the surface. 

6. The hyperboloid of one sheet. There are precisely two lines 
of the surface through every point of it. 
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c. The hyperboloid of two sheets. There are no lines on this 
surface. 

3. A ZERO ROOT AND TWO NONZERO ROOTS HAVING THE SAME 
SIGN. 

a. The elliptic paraboloid. There are no lines on this surface. 

b. The elliptic cylinder. There is exactly one line of the surface 
through each point of it. This surface may reduce to a point or 
it may be imaginary. 

4. A ZERO ROOT AND TWO NONZERO ROOTS HAVING OPPOSITE 
SIGNS. 

a. The hyperbolic paraboloid. There are two distinct lines of 
the surface S through each point of S. 

b. The hyperbolic cylinder. One and only one line of S passes 
through each point of S. 

c. Two distinct planes through a line. There are infinitely 
many lines of S through each point of S. 

5. THERE is ONLY ONE NONZERO CHARACTERISTIC ROOT. 

a. The parabolic cylinder. One and only one line of S passes 
through each point of S. 

b. Two distinct parallel planes. There are infinitely many 
lines on S through each point of S. 

c. Two coincident planes. This is the only place in the classi- 
fication where lines and sets of characteristic roots are not 
adequate to separate cases. It is evident that the line criterion 
separates (c) and (a) into distinct types. The separation of (fe) 
and (c) requires a count of planes rather than of lines. 



CHAPTER 6 
THEORY OF MATRICES 

1. Matrices. A rectangular array is called a matrix. An m 
by n matrix is an array 



l 



(1) 






The horizontal lines in this array are called its rows. They are 
n-dimensional vectors and the ith row is (an, . . . , a n ). Every 
(row) vector is thus a one by n matrix. 

The vertical lines in A are called its columns. The jth column 
of A is 



(2) 



and A has n columns. Each column may be thought of as being 
an m-dimensional column (i.e. y vertically written) vector, and 
is an m by one matrix. 

The scalar in the ith row and jih column of A has been desig- 
nated, by implication, as a/, where the first subscript always will 
indicate the label of the row and the second subscript the label 
of the column in which ay appears. It will be convenient to use 
the notation 



(3) 



A = 



(i = 1, 



,m;j = 1, . . . , n), 



instead of the more cumbersome notation of formula (1). 

A matrix A is called a square matrix if it has as many rows as 
columns, that is, A is an n by n matrix. A square matrix having 
n rows and columns is called an n-rowed square matrix or a square 
matrix of order n. 

77 
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ORAL EXERCISES 

1. Read off the elements an, a 2 4, #51, &26, a 3 6, a w> a ^ m the following 
five by six matrix: 



2 


-1 


3 


1 


4 


5^ 





-2 


-4 


6 


7 


8 


1 





-2 


4 


3 


-5 


-1 


2 


1 


-2 


1 


1 



0020-1 
2. Read off the third row and the fourth column of this matrix. 

2. Addition and scalar multiplication. The sum A + B of 

two rectangular matrices A and B is defined only when A and 
B have the same size. If 

(4) A = (a</), B = (ft*) 

(i = 1, . . . ,m;'j = 1, . . . ,n) 

then the sum of the m by n matrix A and the m by n matrix B 
is the m by n matrix 

(5) C = A + B = (*,-), ca = an + 6 t7 

(i = 1, . . . ,m;j = 1, . . . , n). 

In words, matrices are added by adding corresponding elements. 

Lemmas 1 to 5 of Chap. 1 are special cases of the corresponding 
properties of matrix addition. The student should formulate and 
verify the matrix properties. Note that the zero m by n matrix 
is the matrix whose elements are all zero. We shall use 
the symbol for such a matrix no matter what its size is, and the 
size will always be given by the context. The matrix A is the 
matrix whose elements are the negatives of the elements of A ; 
B A is obtained by subtracting the elements of A from those 
of B. 

If a is a number and A is an m by n matrix, the scalar product 
is the m by n matrix obtained by multiplying every element of 
A by a. Thus if A = (a t /), the element in the iih row and jih 
column of aA is aa^. Clearly 1A = A, ( 1)A = A, oA = 
if a = 0. We also have the properties 

(6) o(6A) = (ab)A, (a + V)A = aA + bA, 

a(A + B) = aA + aB 

for all numbers a, b and all mby n matrices A, B. The verifica- 
tion of these properties is very simple and will not be given here. 



A = 


r 
2 

\4 


| 
1 

2 


2 


1 


4 ) 

O/ 






/ 2 


1 


3> 






/ 4 


1 


2 
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ORAL EXERCISES 

1. Form the matrix sums A + B in the following cases: 



(a) 



-2 1 
-1 

2. Give A and B for the matrices A and B of Oral Exercise 1. 

3. Give 2A and 35 for the matrices A and B of Oral Exercise 1. 

4. Give 2A + B for the matrices A and B of Oral Exercise 1. 

3. Matrix multiplication. The product AB of two rectangular 
matrices A and B is defined only if the number of columns of A 
is equal to the number of rows of B. 

Let us then assume that A is an m by n matrix and that B is 
an n by t matrix. Then the ith row of A is an n-dimensional 
vector 

(7) A t - = (an, aw, . . . , a in ). 

The kth column of B is also an n-dimensional (column) vector, 
which we may write horizontally as 



(8) ' B<*> = (6u, 6 2 ,, . . . , b nk ). 
The inner product 

(9) g ik = A; B<*> = a a &u + + a in b nk 

of the ith row of A and the fcth column of B is defined for i 1, 
. . . , m and fc = 1, ...,/, and therefore the m by i matrix 

<? = (gik) (i = 1, . . . , m; fc = 1, . . . , /) 
is also defined. We call (7 the product of A and 5 and write 

G = AB. 

It is a matrix whose element in the ith row and fcth column is the 
inner product of the ith row of A and the fcth column of B. 

Lemma 1. Matrix multiplication is associative, that is, (AB)C 
= A(BC) for all m 6y n matrices A, n by t matrices B, anrf t by s 
matrices C. 
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For we may write A = (a t /), B = (&,-*), C = (ck P ), where i = 1, 
. . . , m;j = 1, . . . , n; k = 1, . . . , t; p = 1, . . . , s. The 



n 

general element of A B is y a^b,* and that of (AB)C is d tp = 

y-i 

t n t 

/ (/ ciijbjkj Ck P . The general element of BC is ^ 



/t 

that of A(jB(7) is0 t > = N a t /f N ftytCfcpJ. The two finite double 



sums c?i P and gfi p are sums of exactly the same products a^ 

are equal for all values of i = 1, . . . ,m,andofp = 1, . . . , s. 

This proves that A (5(7) = (AB)C. 

Lemma 2. Matrix multiplication is distributive with addition, 
i.e., 

(10) A(B + C) = AB + AC, (B + C)D = BD + CD 

for all m by n matrices A, n by t matrices B and C, and i by s 
matrices D. 

For let A = (a,,-), B = (6,-*), C = (c#), where i = 1, . . . , m; 
j = 1, . . . , n; k = 1, . . . , . Then the element in the zth 

n 

row and fcth column of A(B + C) is } a a(bjk + Cjk) and the 



matrix equality A(B + (7) = AJ3 + AC is equivalent to the 
formula 

n n n 

(11) 2^ a#(&/* + c/ft) = ^ a^-b/fc + 2^ OtfC/jt. 

y - 1 y = i j - 1 

This formula is evidently correct. The proof that (B + <?)) = 
BD + CD is carried out similarly. 

EXERCISES 

1. Form the matrix product AB in the following cases: 



(\ A I -*~ * * -O 

a) 




/ 2 1 3 4\ 
= ( -1 2 0) 
\ 3 -1 2 -3/ 



J5 - 2 - 1 

- 1 2 - 
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(c) 



(-1 2 1 0\ 
2 -1 1 3), 
-1 1 -1 I/ 



B 




2. Form the matrix product BA in those cases of Exercise 1 where it 
exists. 

3. Form the matrix products A(BC) and (AB)C if A and B are given 
by part (6) of Exercise 1 and 

1 - 

2 ( 
-1 2 

3 - 



C = 



4. Form A(B + C) and AB, AC, AB + AC if A and B are given by 
part (c) of Exercise 1 and 

1 -3 

-2 -1 

2 -1 -1 
-1 1 2> 

4. Transposition. If the rows and columns of an m by n 
matrix A are interchanged, the result is an n by m matrix A* 
(read A transpose) called the transpose of A. The operation of 
transposition thus consists of writing the column vectors of A as 
row vectors of A*. If an is the element in the ith row and jth 
column of A, then an is the element in the jth row and ith column 
of A*. 

If P = (xi, . . . , Xn) and Q = (T/ I; . . . , y n ), the matrix 
product PQ* of the one by n matrix P and the n by one matrix 
Q* is given by 



2/1 



x n y n . 



Hence, the inner product P Q of two n-dimensional vectors P 
and Q coincides with the matrix product PQ*. 

Theorem 1. The transpose of a product of matrices is the 
product of the transposes of the factors in reverse order, that is, 



(12) 



(AB)* B*A*. 
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For let A = (ay) and B = (6^), where i = 1, . . . , w; 
j = 1, . . . , n\ k = 1, . . . , t. Then the element in the kth 

n 

row and iih column of (AB)* is dk = /, a tfb/*- The element in 

;=i 

the fcth row and jih column of J3* is &#, the element in the jth 
row and iih column of A* is a t; -, and therefore the element in the 

n 

kih row and iih column of B*A* is y &/&# = c. Hence, 

(AB)* = #*4*. 

ORAL EXERCISE 

Show that if A is an m by n matrix then C = ^4^1* = C*. 

5. Special matrices. The elements a of any matrix A = (an) 
are called the diagonal elements of A and the line of these ele- 
ments is called the diagonal of A. The elements a t -/ with j > i 
(column label greater than row label) are said t<j lie above the 
diagonal of A and those with j < i are said to lie below the diagonal 
of A. 

A matrix is called a triangular matrix if A is a square matrix 
such that either all of the elements above the diagonal in A are 
zero or all of the elements below the diagonal in A are zero. 
Examples are 



'2 


.0 


-1 

2 




4 
-1 
6 



3 \ 
5 I 

-;/ 


i 

2 
-1 

4 
5 



6 
5 
1 

J 






2 
6 





3 

7 


1 






-8 



If A = (an) is a square matrix such that a# = for every 
3) we call A a diagonal matrix. The notation 

(IS) A = diag {ai, a 2 , . . . , a n } 

will be used for an n-rowed diagonal matrix whose iih diagonal 

element a is the number a*. 

A diagonal matrix whose diagonal elements are all equal is 
called a scalar matrix. The scalar matrix whose diagonal ele- 
ments are all unity is called the identity matrix and will usually 
be designated simply by /. Then we may verify by direct 
multiplication 
(14) I m A = AI n 
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for any m by n matrix A where we have used subscripts on / to 
indicate the sizes of the identity matrices involved. Observe 
that a scalar matrix whose diagonal elements are all a is the 
scalar product al and that the matrix product (a/) A is equal to 
the scalar product a A for any matrix A. 

6. Products in terms of rows and columns. Let A = (a,,) 
and B = (6 y ^-), where i = 1, . . . , m\ j = 1, . . . , n; k = 1, 
. . . , t. Designate the jth row of B by 

Bj = (hi, 6/2, . . . , 6,-n) (j = 1, . . . , n) 
and form 



Then the element in the kih column of the vector C is a t -i&u 
+ a t -2&2* + * + dinbnk. This is the element in the ith row 
and fcth column of ^45, and we have proved the result about rows 
in the following theorem. The result about columns is proved 
similarly. 

Theorem 2. The ith row of AB is that linear combination of the 
rows of B whose coefficients make up the ith row of A. The kth 
column of AB is that linear combination of the columns of A whose 
coefficients make up the kth column of B. 

We apply Theorem 2 in the case where A is a diagonal matrix 
so that each row of A has only one nonzero element and this 
element is in the ith column of A. This yields the following 
result : 

Theorem 3. Let A = diag {ai, . . . , a n } and B be an n by 
n matrix. Then the ith row of AB is ai times the ith row of B. // 
D = diag {di, . . . , d n j the jth column of BD is d/ times the }th 
column of B. 

The result given in formula (14) is clearly a special case of 
Theorem 3. 

EXERCISES 

1. Show that if A is a square matrix such that AD = DA, where D is a 
diagonal matrix having distinct diagonal elements, then A is a diagonal 
matrix. 

2. Use the property of Exercise 1 to prove that if A is an n-rowed 
square matrix commutative with all n-rowed square matrices (f.e, 
having the property that AB = BA for every B) then A is a scalar 
matrix. HINT: Form the products AC } - and C,A where C/ has 1 in the 
first row and jth column and zeros elsewhere. 
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7. Partitioning of a matrix. The elements which appear in r 
of the rows and 5 of the columns of an m by n matrix A form an 
r by s matrix which is called a submatrix of A. The r rows need 
not be adjacent. When all r rows and all 5 columns of a sub- 
matrix of A are adjacent, then we shall refer to the submatrix as 
being a block of elements of A. 

Every matrix A may be partitioned into four blocks, and we 
may write 



(15) A = , , At 

Here A = (a*,-) for i 1, . . . , m and j = 1, . . . , n and the 
partitioning is completely determined when we write 

(16) A! = (a) (i = 1, . . . ,r;j = 1, . . . , s). 
Then we are assuming that 

A 2 = (an) (i = 1, . . . , r; j = s + 1, . . . , n), 

A 4 = (a iy ) (i = r + 1, . . . , m; j = s + 1, . . . , n). 

Let A be partitioned as in formula (15). This partitioning is 
determined by the fact that A i is an r by s block in the first r 
rows and first 5 columns of A. Then we shall say that an n by 
t matrix B is partitioned similarly to A if 

(17) B = (f 3 B 

where BI has s rows. The number q of columns in BI is com- 
pletely arbitrary and so is the number of rows in AI. We then 
multiply A by B and have 



A n 



= ( AlBl 
VUBi 



Formula (18) states that if two matrices are partitioned simi- 
larly we may multiply them as if they were two-rowed square 
matrices whose elements are, of course, not numbers but blocks. 
The formula is very easy to derive. We first note that formula 
(18) states that the block of AB which makes up its first r rows 
and first q columns is A\B\ + A%B$. The first r rows of A make 
up a block of elements that may be designated by (AiA 2 ) and 
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the first q columns of B may be designated by 



n 

The element in the iih row and kth column of A B is y a,-/&# = 

y = i 

a n 

y dijbik + 2, Uijbjk- If i = 1, . . . , r and fc = 1, . . . , 5 the 
y=i y = 3+i 

A. 

sum 2 Uijbjk is the element in the ith row and kth column of 



AuBi and } &#&/* is the element in the ith row and fcth column 

._Lj 

of A 2 J5 3 . This proves that 

(A i Az) ( B j = AiBi + A^Bz. 

The remaining relations in formula (18) are proved similarly. 

It should be observed that formula (18) is the result of focusing 
attention on certain rows of A and certain columns of B and of 

n 

breaking up the general sum ^ a;/&/& into the sum of two partial 

s n 

sums ) dijbjk and ) aabft. The formula is of particular use 

./-/ . Lj 

in cases where a block is a zero matrix or an identity matrix. 
For example, let A and B be n-rowed square matrices, and write 

(19) A = 

In this formula we are using the same symbol / for an r-rowed 
identity matrix and for an (n r)-rowed identity matrix, the 
matrix C is an r by n r matrix, and the matrix is an n r 
by r zero matrix. Then 



It follows that in this case the last n r rows of AB coincide 
with the corresponding rows of B. However, each of the first r 
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rows of AB is the sum of the corresponding row of B and a 
linear combination of the last r rows of B. 

EXERCISES 

1. Show that if 



then A 2 and A 4 are zero matrices. 
2. Show that if 



(A, o VB, o \ = 

\A, Aj\0 B,) 



then AiBi and A\B> 2 are identity matrices. This will be shown later to 
imply that A\Bi = B^Ai. Use this property to prove that A 2 = 0. 

8. Determinants. Every n-rowed square matrix A = (a#) 
has a determinant that is a certain function of the elements of A . 
We designate the determinant of A by |A|, by |a#|, or by 



(21) 



#12 ' ' * a in 



We shall refer to the determinant of an n-rowed square matrix 
as an n-rowed determinant or a determinant of order n. Rec- 
tangular matrices with m 7^ n do not have determinants, but 
every r-rowed square submatrix of any matrix A has a determi- 
nant that is called an r-rowed minor of A. 

We shall define \A\ by an induction on n. If n = 1, the 
matrix A = (an), and we define |A| = an. Assume then that 
we have a definition of all (n l)-rowed determinants. 

Each element a,-,- of an n-rowed square matrix A defines a row 
and a column of A. Delete this row and column of A, and obtain 
a submatrix AH of A. This submatrix AH has n 1 rows and 
columns and thus has a determinant |A#| that is an (n 1)- 
rowed minor of A defined for every i and j. Designate by &>,- 
(not by bis) the number 

(22) (-l)*\Av\ 

and call this signed minor the cof actor of an in A. Then &#.is the 
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element in the jih row and iih column of a matrix 

(23) adj A = (6yO (j, i = 1, . . . , n), 

which is called the adjoint of A. 

Let us now form A (adj A) = C = (CM), where 

(24) CM = ttil&lfc + ' * ' + CLi n b nk 

is obtained by multiplying the elements of the iih row of A by 
the cofactors of the elements of the fcth row of A and adding the 
resulting products. We similarly form (adj A) A = D = (djw), 
where 

is obtained by multiplying the elements in the iih column of A 
by the cofactors of the elements in the kih column and adding 
the resulting products. In particular CH is the sum of the prod- 
ucts of elements in the iih row of A by their cofactors and d kk is 
the sum of the products of the elements in the kih column of A 
by their cofactors. 

It can be shown that the n numbers c and the n numbers d kk 
are all equal. We define the common value of these 2n numbers to 
be the determinant of A. We have thus given not only a definition 
of a determinant but what are called the expansions of it accord- 
ing to any row and any column. 

EXERCISES 



1. Expand the following determinants according to the first row: 



(a) 



(d) 



1 -1 2 




2 1 -2 




2-1 1 1 


3 1-2 


(6) 


-1 1 


(c) 


1101 


1 2 3 




-2 1 




220-1 




1-3-2 1 


2 -1 -2 1 




3-1 2-1 




-1211 




-1101 




1 -3 -1 -2 


(e) 


1232 




21-23 




2-3-1 1 





2. Expand each of the determinants of Exercise 1 according to the 
second row. 

3. Compute the adjoint of each of the following matrices: 

/-I 2 1\ / 1 -1 1\ / -1 2\ 

(a) ( 2 -1 31 (6) ( 2 2-1) (c) (2-1 3) 

\ 1 2 -I/ \-l 1 3/ \-4 2 I/ 
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9. Properties of determinants. We shall assume the following 
properties of determinants. 

Lemma 3. A square matrix and its transpose have the same 
determinant. 

Lemma 4. Let B be the result of interchanging two rows (col- 
umns} of a square matrix A. Then |B| = |A|. 

Lemma 6. Let B be the result of multiplying a row (column) of 
A by a number a. Then |B| = a|A|. 

Note that this result implies that \aA\ = a n \A\. 

Lemma 6. Let B be the result of adding a scalar multiple of a 
row (column) of A to another row (column). Then |B| = |A|. 

Lemma 7. Let a row (column) of A be a scalar multiple of 
another row (column). Then |A| = 0. 

If we multiply the elements of the iih row (column) of A by 
the cof actors of the elements in its kih row (column), the result 
is the expansion according to the kih row of the determinant of a 
matrix B obtained by replacing the kih row (column) of A by 
its iih row (column). Then two rows (columns) of B are equal 
and 1 75 1 = 0. But \B\ is the number CM of formula (24) in the 
row case and is the number dki of formula (25) in the column 
case. Hence, Cik = Cki = for i ^ k and c = da = \A\. This 
yields the following result: 

Theorem 4. Let A be an u-rowed square matrix. Then 

A(arf? A) = (a#A)A = |A|I 

is the scalar product of the n-rowed identity matrix I by the deter- 
minant of A. 

Any square matrix A can be converted into a diagonal matrix 
by a finite sequence of transformations of the types given in 
Lemmas 4 and 6. Indeed if A j 0, we can carry any nonzero 
element ai of A into its first row and column. We can then sub- 
tract multiples of the first row from the remaining rows and 
multiples of the first column from the remaining columns and 
convert our matrix into 





where A\ has n 1 rows and columns. The same procedure 
may be applied to A \ by applying transformations to the last 
n 1 rows and columns of B], and we ultimately convert A into 
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a diagonal matrix D. If A has been converted into D by trans- 
formations involving an even number of row and column inter- 
changes, then | A | = \D\. Otherwise, we change the sign of a 
row before obtaining the final diagonal matrix D and will have 
|A| = \D\. Then |A| may be computed by the use of 

Lemma 8. The determinant of a triangular matrix T is the 
product of the diagonal elements of T. 

It is clearly sufficient to consider the case where the elements 
above the diagonal in T are all zero. Then 


Si T 

where T\ is a triangular matrix of n 1 rows. Expand T 
according to its first row, and see that |2P| = ailT 7 ^. Apply the 
same procedure to T\, and ultimately obtain the result of the 
theorem. 

We should observe that if the transformations used to convert 
A to a diagonal matrix D such that \A\ = \D\ are applied to any 
square matrix B of any size the resulting matrix C will have the 
property that \B\ = \C\. For C results from B by a finite number 
of transformations of the kind in Lemmas 4, 5, and 6 and there 
are only either an even number of those in Lemma 4 and none of 
those in Lemma 5 or there are an odd number of those in Lemma 
4 and one of those in Lemma 5 with a = 1. We use this 
result to prove 

Lemma 9. Let 



(20) C = 



A 
G B 



AH 
B 



where A a^id B are square matrices. Then 

(27) |C| = |D| = |A| |B|. 

To prove this result, we apply a sequence of transformations 
to C and to D which convert -A to a diagonal matrix A Q and B to a 
diagonal matrix B Q such that \A\ = \A Q \ and \B\ |JBo|. These 
transformations convert C and D, respectively, to 

A 



(A H \ 
\0 Bo/ 



Then |C | = |C|, |D | = |D|. But the new matrices C and D 
are triangular matrices and their diagonal elements are the 
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diagonal elements of A and B Q . Then |Co| = |Ao| \B Q \ = |/) | 
and |C| = |D| = \A\ \B\ as desired. 

We use Lemma 9 in the proof of the very important result that 
we state as follows: 

Theorem 5. The determinant of a product of n-rowed square 
matrices is the product of the determinants of its factors, i.e., 
|AB| = |A| |B|. 

To prove this result, we consider the 2n-rowed square matrix 



Form the product 
(29) D 



o / -i B 



as in formula (19). Only the first n rows of C have been altered 
and each of these rows is a corresponding row of C plus a linear 
combination of the last n rows of C. By Lemma 0, we have 
\D\ = \C\. It requires n interchanges of rows and n changes of 
signs of rows to replace D by 



(30) G = ln AB 

Then \0\ = (~l) 2n |/)| = |C|. By Lemma 9 we have |0| = 
|/| |AJS| = |AB| and by this same lemma |C| = |A| \B\. 
Hence, |AJ3| = |A|-|B| and our theorem is proved. 

EXERCISES 

1. Use the properties of determinants to reduce each determinant of 
Exercise 1 of Sect. 8 to triangular form and compute its value by the 
use of Lemma 8. 

2. What are the determinants of the following matrices? 







8 
5 



10. The inverse of a matrix. A matrix A is said to be non- 
singular if A is a square matrix and |A| 5^ 0. A square matrix 
A whose determinant is zero is called a singular matrix. 




2 


1 


3 


8 


9 


-3\ 





2 


1 


6 


1 


2 





4 


3 


9 


8 


7 











4 

















-1 


5 





p 








9 


6 


2 



Sec. 10] THEORY OF MATRICES 91 

If A is a nonsingular n-rowed matrix, we define a matrix 

(31) A" 1 = 4-, adj A 

\ A \ 

which we call the inverse of A. By Theorem 4 

(32) AA- 1 = A~ 1 A = I 

is the n-rowed identity matrix. A singular matrix does not have 
an inverse. For if A A" 1 = 7, then |AA->| = |A| \A~ l \ = 1 
and | A | ^ 0. Note that 

(33) |A-'| = |A|-'. 

If A is nonsingular, the matrix equations 

(34) AX = B y YA = C 

have unique solutions X and Y for any given matrices B and C. 
Indeed if AX = B, then A-'(AX) = A~ 1 B = (A-'A^ = /Z 
= X. Similarly if FA = C, then CA- 1 = (FA)A- 1 = Y(AA~ l ) 
= Y. We have proved that the equations of formula (34) have 
the unique solutions 

(35) X = A- 1 ^ F = CA' 1 . 

In particular the equations AX = /and FA = / have the unique 
solution X = F = A" 1 . It follows that AX = I if and only if 
XA = /, that is, X = A" 1 . 

The following results are almost trivial. 

Theorem 6. The inverse of a product is the product of the 
inverses of the factors in reverse order. 

For if C = AiA 2 A r , then CAr l A^r l - Ar 1 = 1, 
that is, C~ l A r ~- l A r -iT l Ai" 1 . In the simple case of two 
factors, we are stating that ABB~ l A~ l = / and therefore 



Theorem 7. The inverse of A* is the transpose of A"" 1 . 

For A A" 1 = /,/* = /= (A -1 )*A* and therefore (A*)" 1 = 

EXERCISE 

Compute A~ l by using formula (31) in each of the following cases: 

/2 -1 1\ /I 2 0\ /I 1 2\ 

(a) (O 4 l), (6) (4 7 ), (<) (l 2 3). 
\0 3 I/ \6 -3 I/ \3 -2 2/ 
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11. Linear systems of equations. A system of ra linear equa- 
tions 

dnXl + ai2#2 4" ' * * + di n X n = 2/1 

j I I 

/rtr\ d%\X\ ~t" d%%X2 "p * * " ~p d%nXn == V2 

(36) 



r &m2#2 i * ' ' i d mn x n y m 

in n unknowns x\, . . . , x n is called a linear system. If A is the 
m by n matrix, A = (a,-/) of the coefficients, and P = (xi, . . . ', 
Xn), Q = (2/1, , J/w)i we may write any linear system as an 
equivalent vector matrix equation 



(37) 
where 



= Q*, 



g* = 



Any system of n equations in n unknowns with \A\ 7^ may 
be solved by determinants or by elimination. Whatever pro- 
cedure is used, the solution must be 

P* - A- 1 ^*. 

It follows that to find the inverse of a given nonsingular matrix 
A we form the equations 



/ oox 
(OQ) 



^nXn = 2/2 



dnlXi + ' ' + d nn Xn = 2/n 

and solve these equations for Xi, . . . , x n in terms of #1, . . . , 
2/ n? thereby obtaining solutions of the form 

Xi = 6ll2/l + ' ' * + &ln2/n 
(39) 2 21 1 2*1 n 

Then, if J5 = (&</), we see that B = A"" 1 is the required inverse 
matrix. 
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ILLUSTRATIVE EXAMPLE 

Compute the inverse of 



CO 3 2 3v 
12121 
-123 5J' 
-1 3 4 7/ 



Solution 

We require the solution of the system of equations 

3x 2 + 2x 3 + 3x 4 = y l 

Xl + 2X 2 + X 3 + 2^4 = 2/2 

xi + 2x 2 + 3x 3 + 5^4 = 2/3 
-Xi + 3x 2 + 4x 3 + 7x 4 = 2/4. 

We first eliminate x\ and obtain 

4o- 2 + 4x 3 + 7x 4 = 2/2 + 2/3 
5x 2 + 5x 3 + 9x 4 = 7/2^-f- 2/4- 

Then 36x 4 35x 4 = x 4 = 4(?/ 2 + 2/4) 5(?/ 2 + 2/3) = 2/2 %3 + 
4?y 4 . We also compute x 2 + x ;{ + 2x 4 y\ 7/3 and have 3(2/4 2/3) 

2/i = 3x 2 + 3x 3 + 6x 4 (3x 2 + 2x 3 + 3x 4 ) = x 3 3x 4 , so that x 3 = 
3ij 4 - 82/3 - 2/i + 32/2 + 157/3 - 12?/ 4 = -2/1 + 3?/ 2 + 12t/ 3 - 9y 4 . We 
now obtain x 2 = 2/4 - 2/3 - ( 2/i + 3?/ 2 + 127/ 3 9?/ 4 ) - 2( y 2 ^2/3 
+ 4?/ 4 ) = 2/1 2/2 ~ 3?/ 3 + 27/ 4 , xi = 2(7/1 2/2 37/ 3 + 27/ 4 ) 

(-2/1 + 37/ 2 + 127/a - 97/ 4 ) 2(~7/ 2 57/3 + 47/ 4 ) = -2/1 + 2/2 + 47/3 

37/ 4 . It follows that 

\ 




The result should be checked by computing A A' 1 = /. 

EXERCISES 

1. Compute A" 1 by the method above for the exercise of Sec. 11. 

2. Compute A~ l for each of the following matrices: 










-1 


2 


-2 


0> 





-1 


1 


1 


3 


1 


-4 


3 





2 


-1 


_.. 1 > 


-1 


1 





1> 


1 








1 


4 


3 


-3 


3 





-1 


2 


1' 
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12. Homogeneous systems. A linear system 

OllXi + ' ' + dinXn = 
t n = 



+ * * ' + a mn X n = 

is called a homogeneous system. It requires that 

(41) AP* = 

and therefore that P = (x\, . . . , x n ) be a vector such that 

(42) 4<-P = (* = 1, . . . ,n), 

where A is the ith row of F*. Thus we see that a homogeneous 
system of equations proposes the problem of finding an n-dimen- 
sional vector P which is orthogonal to m given n-dimensional 
vectors. The zero vector is a trivial solution of the problem but 
we are usually interested in finding a nontrivial solution, i.e., a 
nonzero vector P. 

Theorem 8. A homogeneous system of m equations in n 
unknowns has a nontrivial solution ifm < n. 

A single equation a\x\ + * + a m x n = in at least two 
unknowns has a nontrivial solution. For we put #3 = = 
x n = and see that if ai = a 2 = the vector (1, 1, 0, 0, . . . , 0) 
is a nontrivial solution. But if one of ai and a% is not zero the 
vector ( a 2 , ai, 0, . . . , 0) is a nontrivial solution. We now 
make an induction on m and assume that the theorem is true for 
m 1 equations in more than m 1 unknowns. ^/ 

The theorem is surely trivial iFalTcoefficients are zero. Hence, 
we may assume that at least one coefficient is not zero. There 
is surely no loss of generality if we assume that this coefficient 
is an. This amounts to a relabeling of equations and variables 
if necessary. We may then subtract multiples air^a^ of the 
first equation from all the other equations and obtain our equiva- 
lent system. The first equation may be written as 

(43) x\ = air l (aizx<i + + ai n x w ), 

and the remaining equations are free of #1. We then have m 1 
equations in n 1 unknowns. These equations are homogene- 
ous and have a nontrivial solution (#2, . , x n ). If. x\ is 
defined by formula (43), the vector (xi, . . . , #) is a nontrivial 
solution of our original system. 
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Theorem 9. A homogeneous system AP* = of n equations in 
n unknowns has a nontrivial solution P ^ if and only if \ A| = 0. 

For if \A\ 7 0, the equation AP* = implies that A~ l (AP*) 
= P* = 0, P = 0. Let us then suppose that \A\ = 0. If 
n = 1, the system becomes a\x\ = 0, where \A\ = a\ = and 
XL j& defines a nontrivial solution. Assume then that the 
theorem is true for the case of n 1 equations in n 1 unknowns. 
As in the proof of Theorem 8, we may assume that a\\ ^ and 
pass to an equivalent system with matrix 

^ _ /an BI 



where B\ is an (n l)-dimensional vector, the zero represents a 
column of n 1 zeros, and 2 is (n l)-rowed square matrix. 
By Lemma 6 \A\ = |B| = 0, and by Lemma 9 |J5| = an|J3 2 |. 
Since an 9^ 0, the number |jB 2 | = 0. We may then solve the 
homogeneous system of n 1 equations in n 1 unknowns 
with #2 as matrix and obtain a nontrivial solution (#2, . . . , x n ). 
Determine x\ by formula (43), and obtain a nontrivial solution 
(xi, . . . , Xn) of the original system. 

13. The characteristic equation. If A = (a#) is any n-rowed 
square matrix, we may subtract x from the diagonal elements of 
A and thus form the matrix 



(A - xl) = 




Compute the determinant of this matrix. It is a polynomial 
in x whose leading coefficient is (~l) n . It follows that the 
polynomial 

f(x) = t-\Y\A - xl\ = & + aix~ l + + On 

has leading coefficient unity. We call f(x) the characteristic 
determinant of A. It is actually equal to \xl A\. 

The principal minors of A are the determinants of those square 
submatrices of A whose diagonal elements are all diagonal ele- 
ments of A. It can be shown that the coefficient a of f(x) is 
actually equal to ( 1)* times the sum of all i-rowed principal 
minors of A. We shall not prove this result, but it is important 
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to observe that 

ai = an + + a nn 

is the sum of the diagonal elements of A and that 

/(O) =an== |_ A | = (_1)|A|. 

The equation /(#) = is called the characteristic equation of A 
and its roots are called the characteristic roots of A. 

EXERCISE 

Compute the characteristic determinants and characteristic roots of 
the following matrices: 



CO 1 0\ 
001) 
2 1 -2/ 



(6) 




1 

1 



4 -6 

3 

33) (/) | - 

3 



"U 

\ 

3 ) 

O/ 




14. Similar matrices. Two n-rowed square matrices A and B 
are said to be similar if there exists a nonsingular matrix L such 

that 

B = L- 



Theorem 10. Similar matrices have the same characteristic 
determinants and consequently the same characteristic roots. 

For \B - xl\ = IL-^L - *7| = IL-'CA - */)| = IL- 1 ! |A 
xl\ \L\ = \A xl\. Then if A and B are similar, their 
characteristic determinants are the same, and we have proved 
the theorem. 

We also can prove the following: 

Theorem 11. Let A be similar to a diagonal matrix D. Then 
the diagonal elements of D are the characteristic roots of A. 

For if D = diag {di, . . . , d n ], the matrix D xl = diag 
{di x y . . . , d n ~ x] and, by Lemma 8, \D xl\ = (di x) 
' (d n #). But |Z) #7 1 = \A a;7| and di, . . . , d n 
must be the characteristic roots of A. 

Theorem 11 states that, while it is actually not true that all 
square matrices are similar to diagonal matrices, if a matrix A 
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is similar to a diagonal matrix D the matrix D is uniquely deter- 
mined up to an arbitrary permutation of its diagonal elements. 
Suppose then that A is given and that we determine the charac- 
teristic roots of A and so prescribe Z). We propose to try to find 
a nonsingular matrix L such that 

AL = LD. 

The jth column of AL is found by our row by column rule of 
matrix multiplication to be the product AQ (J \ where Q 0) is the 
jth column of L. The jth column of LD is the product Q (j) d 7 - of 
the jth column of L by the jth diagonal element c/ 7 of D. Then 
AL = LD if and only if AQ& = d/Q'>. But this latter equation 
is equivalent to the equation 

(A - d,./)<p - 0. 

We are thus led to the determination of each column of L as 
a vector that is a solution of a linear homogeneous system with 
matrix A dl, where d is a corresponding characteristic root. 
Since \A dl\ = for each d, a corresponding nonzero vector 
Q can always be determined. A set of solutions Q (1) , . . . , Q (n) , 
which are the columns of a nonsingular matrix, can be found only 
when A is similar to a diagonal matrix. 

EXERCISE 

Find a nonsingular matrix L such that L~ 1 AL is a diagonal matrix in 
each of the following cases where L exists: 




(6) 



(d) 



/ 1 Ov 

001 
\-l -3 3/ 

r l \ 
i i ii 

\0 1 O/ 



15. Real symmetric matrices. A matrix is called a symmetric 
matrix if it is equal to its transpose. Then a symmetric matrix is 
necessarily a square matrix A = (a t v) such that a,-/ = a# for every 
i and j. 

Theorem 12. The characteristic roots of a real symmetric 
matrix are all real. 
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For let a + bi be a characteristic root of A where a and b are 
real and i 2 = 1. Then A (a + bi)I is singular and so is the 
product 



B = [A - (a + bi}I][A - (a - 

= A 2 - 2aA + (a 2 + 6 2 )7 = (A - a/) 2 + 6 2 7. 

The matrix B is a singular matrix with real elements and, by 
Theorem 9, there exists a real nonzero vector P such that PB = 0. 
Then 



PBP* = JP(A - a7) 2 P* + 6 2 PP* = QQ* + 6 2 PP* = 0, 

where Q = P(A - a/) and Q* = (A - a7)*P* = (A - a7)P? 
since A is symmetric. Since Q is a real vector, QQ* = Q Q ^ 0. 
But PP* = P P > and QQ* + 6 2 PP* = is impossible if 6 ^ 0. 
Hence, b = 0, and we have proved that all characteristic roots 
a + bi of A are real. 

16. Orthogonal matrices. A matrix L is called an orthogonal 
matrix if L is an n-rowed square matrix such that LL* I is the 
ri-rowed identity matrix. Then 

LL* = L*L - /, L" 1 - L*. 

A set of ra distinct vectors Pi, . . . , P m is said to consist of 
pairwise orthogonal vectors if P P 7 - = for i ^ j, that is, any 
two distinct vectors in the set are orthogonal. Let PI, . . . , P n 
be the rows of an orthogonal matrix L. Then the element in the 
iih row and jth column of LL* is P< P,-. Since LL* = /, we see 
that 

Pi-Pi = 1, Pi Pf = (i ^ j] i y j = 1, . . . , ri). 

Thus the rows of an orthogonal matrix are pairwise orthogonal unit 
vectors. Conversely if the rows of L are n pairwise orthogonal 
n-dimensional unit vectors, L is an orthogonal matrix. 

The columns of L are the rows of L*. Since L*L = /, we see 
that the columns of an orthogonal matrix are also n pairwise 
orthogonal unit vectors. Conversely, if the columns of a square 
matrix L are pairwise orthogonal unit vectors, L is an orthogonal 
matrix. ( J^UA o^ ) 

Theorem 13. Let L and M be u-rowed orthogonal matrices. 
Then LM is an orthogonal matrix. 

For LL* = 7, MM* = /, (LM)(LM)* = LMM*L* = L/L* 
- LL* = 7. 
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Theorem 14. Let L and M be orthogonal matrices. Then 



/ 

\ 



'L 

U) 

is an orthogonal matrix. 

For proof we merely compute 

/L o VL o y /L o VL* o \ /LL* o \ 

\0 M/\0 ^7 \ Af/\0 M*/ \0 MM*/ 
Note that in Theorem|3 the matrices L and M must have the 

same size but that this is not necessary in Theorem 14. 

Theorem 15. Let PI, . . . , P m be m pairwise orthogonal 

\\-dimensional nonzero vectors. Then m ^ n and there exists an 

orthogonal matrix whose first m rows (columns) are scalar multiples 

nf P, P 

OJ 1 1, . . . , m. 

For the equations Pi P = 0, . . . , P m P = form a 
homoj eneous system of m linear equations in the coordinates 
Xi, . . . , x n of P. If m < n, there exists a solution P ? of 
this system. It follows that if m < n we can find a vector 
Pm+i T^ such that PI, . . . , P m +i are pairwise orthogonal. If 
m + 1 < n, we can extend the set of vectors again and thus 
ultimately obtain vectors P w +i, . . . , P n such that Pi, . . . , P n 
are n pairwise orthogonal nonzero vectors. Define the scalar 
multiples 



Then C/i, . . . , U n are n pairwise orthogonal unit vectors and 
are the rows of an orthogonal matrix L as well as the columns of 
the orthogonal matrix L*. If m ^ n, we form L as above and 
Pi PI = for i = 1, . . . , n, and j > n implies that C7 t - P/ = 
for i 1, . . . , n and j > n. But then LP/* = 0. However, 
L is nonsingular and necessarily P/* = 0, contrary to hypothesis. 
It follows that m ^ n. 

EXERCISE 

Find an orthogonal matrix whose first two rows are scalar multiples 
of Pi and P 2 in the following cases: 
(a) P a = (1,2,-2),P 2 = (2,1,2) 
(6) P 1= = (1,0,1),P 2 = (1,0, -1) 

(c) P! = (1,2,3),P 2 = (1, -5,3) 

(d) P! = (1,1,1),P 2 = (0,1, -1) 

(e) Pi = (6, 3, 2), P, (1, 0, -3) 
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17* Orthogonal reduction of a symmetric matrix. A real 
symmetric matrix A is said to be orthogonally equivalent to a real 
symmetric matrix B if there exists an orthogonal matrix L with 
real elements such that 

B = L*AL. 

Then B = L~ 1 AL is similar to A and has the same characteristic 
roots as A. 

Theorem 16. Every real symmetric matrix A is orthogonally 
equivalent to the diagonal matrix D = diag {di, . . . , d n }, where 
di, . . . , d n are the characteristic roots of A arranged in any 
prescribed order. 

For the equation (A di/)P = has a nonzero solution P. 
Then P may be taken to be a unit vector and, by Theorem 15, 
we may determine an orthogonal matrix L whose first column 
is P = C/i. Designate the jth column of L by Uj, and see that 



Then Ui*AUi = bn is zero for i ^ 1, and 6n = d\. It follows 
that the element bn in the tth row and first column of B = L *AL 
is zero if i 5^ 1 and bn = d\. Since #* = L *A*L = L ML 
= B is a symmetric matrix, we have 



where A i is an (n l)-ro\ved real symmetric matrix. Also 
- x 



- xl\ = 



Ai - xl 



= |4 - a:/ 1, 



and d2, . . . , d n are the characteristic roots of A\. By the 
proof above applied to A\ there exists an (n l)-rowed orthog- 
onal matrix L\ such that 



where A 2 is an (n 2) -rowed real symmetric matrix whose 
characteristic roots must be d 3 , . . . , d n . Then 

N = Lo 
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is the product of two orthogonal matrices and is an orthogonal 
matrix having the property that 

-(; "Hi -)-(' <) 

\o L I / \o L! / \p Az y 

A repetition of this process yields the result desired. 

The proof above is not intended to provide a procedure for 
finding an orthogonal matrix L such that L*AL = diag {d\ 
. . . , d n }. The procedure that is indicated by the proof is not 
recommended, because it is certainly too clumsy for use. The 
best procedure to employ is that of Sec. 14 where we keep in 
mind at every stage that the columns of the matrix L which we 
are determining are pairwise orthogonal nonzero vectors (see 
Sec. 4 of Chap. 7 for examples of this procedure). 

18. Uniqueness of characteristic unit vectors. If d, is a 
characteristic root of a symmetric matrix A and P/ is a solution 
of the equation (A d//)P/* = 0, we call P, a characteristic 
vector of A corresponding to the root d,-. We now prove the 
following : 

Theorem 17. Let dj be a simple root of the characteristic equa- 
tion of A. Then the corresponding characteristic unit vector is 
unique apart from sign. 

For we know that there exists an orthogonal matrix L such 
thatLML = diag {di, . . . ,d n }. Then AP,-* = dy/P/* is equiv- 
alent to L*ALL*P/* = d ; L*P ; * and therefore to 



(diag{d l9 . . . , 
where Q/* = L*P/*, P/* = LQ/*. But then 

(diag {di - dj, d 2 d/, . . . , d d/})Q/* = 0. 

This means that if Q/ = (Xi, . . . , X n ), then X*c,- = 0, where 
Ci di dj. Then X t - = 0, for i 7* j, Q/* is a vector with X/ in 
the jth row and zeros elsewhere, and P 7 -* = X,-J7y* where Uj* is 
the jth column of L. If P/ is a unit vector, it follows that 
Pi Py = Xy 2 C7,- Uj = Xy 2 = 1, X/ = + 1, Py = C7 y as desired. 

We note that if dj is a double root and thus is equal to d& for 
some value of fc, then Q/ is a vector with X/ in the jth row, X* 
in the fcth row and zeros elsewhere, and P, = X;7/ + X^C/*. 
Thus PJ is not unique. Indeed we may show similarly that 
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Pk = fjijUj + UkUk. Then 

p. . p. = x, 2 + X* 2 = 1, Qj - Qj - M/ 2 + p. k * = 1 
since Uj - U = (7* // = 0. If also P 7 P^ = 0, we have 

X//iy ~f" X^jLtfc == 0. 

It follows that the matrix 

% 



\ 



is an orthogonal matrix and therefore P/ and PA are expressible 
in terms of Uj and [7& by means of a two-rowed orthogonal matrix 
in the case where dj is a double root. We shall not consider 
further cases of this study but pass on to its application to three- 
dimensional geometry. 



CHAPTER 7 
ROTATIONS OF AXES AND APPLICATIONS 

1. Orthogonal transformations. If four points in space are 
not coplanar, they determine a tetrahedron, i.e., a pyramid with 
a triangular base. In particular, the origin and the three unit 
points 

U= (1,0,0), F= (0,1,0), W= (0,0,1), 
on the positive rays of the coordinate axes, determine a tetra- 
hedron the vertex of which is the origin and which is called the 
tetrahedron of reference. Conversely, if a tetrahedron of reference 
is given, the coordinate system is completely determined. For 
the tetrahedron determines the origin, the unit distance, and the 
positive rays on the coordinate axes. 

Let a coordinate system called the initial system be given so 
that every point P has initial coordinates x, y, z and P is the linear 
combination 

(1) P - (x, y, z) = xU + yV + zW , 

with coefficients the coordinates of P. We propose the study of 
all other rectangular Cartesian coordinate systems with the same 
origin as the initial system. Every such transformed system 
is determined by three new unit vectors /', F', W and every 
point P has transformed coordinates x', y', z f . Then 

(2) P = x'U' + y'V + z'W 

where the coefficients x' ', y', z 1 are the transformed coordinates of 
P. Let the initial coordinates of U' 9 V, W be given by 

(3) U' = (Xi, MI, n), V = (X 2 , M2, 2 ), W 9 = (X 8 , M3, *s). 

Then P = x'(\i, MI, "i) + 2/'(X 2 , M2, ^2) + 2 7 (X 3 , MS, v 3 ) = (x, y, z] 
and the relations between the initial and transformed coordinates 
of all points are given completely by the set of equations 
x = 

(4) y = 

Z = 

103 
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These equations may be written in matrix form as 





(5) 



Then the columns of L are the pairwise orthogonal unit vectors 
I}', V, Wj and L is an orthogonal matrix. 

A system of equations of the form given in formula (5) with L 
an orthogonal matrix is called an orthogonal transformation of 
coordinates. We have proved that any two rectangular coordi- 
nate systems with the same origin are related by an orthogonal 
transformation. Conversely, an orthogonal transformation may 
be interpreted as relating two coordinate systems in which the 
columns of the matrix L give the initial coordinates of the unit 
vectors on the transformed coordinate axes. 

If L is orthogonal, Lr l = L. Then the solved form of formula 
(5) is 

(6) 

These equations may be written in full as 

x' = \\x + MI?/ + v\z 

(7) y' = X 2 # + niy + viz 

z f = \ 3 x + n*y + v*z. 

In our applications of this theory to quadric surfaces we will 
determine vectors of integers which are scalar multiples of 
U', V, W '. The scalar multipliers are square roots of rational 
numbers and appear as common denominators of each equation 
in formula (7); this is not true of formula (4). Then it will be 
more convenient for us to express our final answers in the form 
given by formula (7). It should be noted that 

U = (1, 0, 0) -> (Xi, X,, X 8 ) 
7 - (0, 1, 0) -> (MI, M, MS) 

w = (o, o, i) - K ? 2 , vd 

where the set of coordinates of each of our points given after the 
arrow is the set of x', y', z 1 coordinates. These coordinates are 
obtained by substitution of the set of coordinates before the arrow 
for x, y, z in formula (7). 
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We have now shown that the columns of L are the x, y, z 
coordinates of the unit vectors on the x', y', z' axes and that the 
rows of L are the x' * y', z f coordinates of the unit vectors on the 
x, y, z axes. 

EXERCISES 

1. Convert the matrices 



( 
1 2 2 1, B = 



/2 -2 1\ 

= 1 2 2, 

\2 1 -2/ 




into orthogonal matrices that are scalar multiples aA, @B of the given 
matrices. 

2. Let L in (5) be the orthogonal matrix a A of Exercise 1. Give 
the x', y', z' coordinates of the points whose x, y, z coordinates are 

(a) (-1,2,2) (d) (-1,1,0) 

(b) (0,1,1) (e) (-1,1, -2) 

(c) (1, 1, 1) (/) (6, -3, 2) 

3. Let the coordinates given in Exercise 2 be the x', y', z f coordinates 
of points. What are their x, y, z coordinates? 

4. Give an x', ?/', z' equation of each of the following planes using the 
transformation of Exercise 2 : 

(a) 2x + y + 2z = 1 (c) 2x - y = 1 

(b) 2x - 2y + z = 3 (d) x + y + z = 

6. Give the x', y', z 1 equations in parametric form of the lines joining 
the following pairs of points where the transformation of Exercise 2 is 
used and the coordinates given are #, y, z coordinates: 

(a) (-1, 2, 2), (2, 1, 3) (d) (3, 4, 1), (-3, -4, -1) 

(6) (1, 2, -2), (0, 0, 0) (e) (1, 0, -1), (0, 1, 0) 

(c) (2, 1, 2), (-2, 2, 1) (/) (2, 3, -4), (1, 2, 3) 

6. Let L be the orthogonal matrix (}B of Exercise 1. Give the follow- 
ing equations: 

(a) The x', y', z f equation of the plane x = 0. 
(6) The x'j y'j z r equation of the plane y = 7. 

(c) The or', y' y z' equation of the plane 2 = 1. 

(d) The x, y, z equation of the plane y' = 0. 

(e) The x, y, z equation of the plane z' = 7. 
(/) The x, y, z equation of the plane x' = 1. 
(g) The x, ?/, 2 equations of the line x' = y' == 1. 
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2. Products of orthogonal transformations. The result of 
applying two successive orthogonal transformations of coordi- 
nates is an orthogonal transformation of coordinates called their 
product. Thus the product of the orthogonal transformation 





(8) 

with matrix L and the orthogonal transformation 

(x'\ fx' 

y'} = M[y' 
z'/ V 

with matrix M is the product orthogonal transformation 

i x \ 

(10) [y U A' 

W 

whose matrix N is obtained by substitution in formula (8) of the 
values of x', y f , z' in terms of x n ', y", z" as given by formula (9). 
But then 




= L 




= (LM) 



1 X "\ 

(y") 

\z"/ 



and we have proved that LM = A r , that is, the matrix of a product 
of two orthogonal transformations is the product of the matrices of 
the transformations. 

EXERCISE 

Give the equations of the orthogonal transformation with matrix 
(aA)(0B) of Exercise 1 of Sec. 1. 

3. Reflections and rotations. If the direction on a coordinate 
axis is changed, the resulting transformation is an orthogonal 
transformation of coordinates defined by one of the matrices 

-1 0\ /I 0\ /I 0> 
1 , (0 -1 01, 10 1 
I/ \0 I/ \0 -Ij 

Such a transformation is called a reflection of axes. 
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Suppose now that the vertex of the tetrahedron of reference 
is held fixed in space and that the tetrahedron is rotated about 
this vertex to some new position. The resulting orthogonal 
transformation is called a rotation of axes. It should be evident 
that every product of a finite number of rotations of axes is a 
rotation of axes. 

A rotation of axes about the z axis may be conceived of as an 
ordinary rotation of axes of plane analytic geometry. Its equa- 
tions are 

x = ux' vy' 

(11) y = vx' + uy' 

z = z' 

where u = cos 0, v = sin 0, and the angle of rotation 6 is measured 
in a counterclockwise direction from the unit point U to the unit 
point U'. The matrix form of this rotation is given by formula 
(5) where 

/u -v ON 

(12) L = [ v u 

\0 

and the determinant of L is 1. Similarly the rotations about the 
x and y axes are space rotations whose matrices 



(13) 



are orthogonal matrices having determinant 1. We shall call 
any rotation of axes about a coordinate axis a planar rotation of 
axes and shall discuss such rotations further in Chap. 8. 

The product of two reflections of axes is a rotation of axes; 
for example, 



o o\ A o oN 








I -1 



is the matrix of a rotation of axes. Indeed the product of two 
reflections of distinct axes may be seen geometrically to be the 
planar rotation about the remaining axis through 180. In the 
example above, the result of replacing y by y and z by z is 
the same as the rotation about the x axis through 180. 
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A reflection of axes is not a rotation of axes. For a rotation of 
axes that carries two of the points C7, V, W into a corresponding 
pair of the set U', V, W' is a rigid motion of the tetrahedron 
and must carry the remaining point of C7, V, W into the remaining 
point of U'j V, W. We shall use this result in the derivation of 
the following basic theorem. 

Theorem 1. An orthogonal transformation is a rotation of axes 
if and only if the determinant of its matrix is 1. Every rotation 
of axes can be expressed as a product of three planar rotations and 
every orthogonal transformation not a rotation of axes can be 
expressed as the product of a rotation of axes and a reflection of axes. 

For consider an orthogonal transformation in which C7, V, W 
are the unit vectors on the positive rays of the x, y, z axes and 
f/', V, W are the unit vectors on the positive rays of the 
x', y', z 1 axes. Apply a planar rotation about the z axis which 
carries the x axis into the line of intersection of the x', y f plane 
with the x, y plane. This rotation determines an xi, ?/i, z\ coor- 
dinate system in which the z\ axis coincides with the z axis. 
Rotate about the x\ axis, which is in the x' , y' plane, so that the 
z\ axis is carried into the x'y' plane. The result is an # 2 , 2/2, 2 
coordinate system in which the # 2 , 2/2 plane coincides with the 
x' j y' plane. If the angle from the corresponding unit vector t/ 2 
on the #2 axis to the unit vector V% on the ?/ 2 axis is not measured 
in the same direction as the angle from U' to V, we increase the 
planar rotation just described by 180 and replace F 2 by F 2 
and thus restore the angular orientation of the x'y' plane. Thus 
we may assume that the angle from t/ 2 to F 2 is measured in the 
same direction as that from V to V. A rotation about the z 2 
axis that carries Ut into U' also carries V% into V. 

We have now proved the existence of three planar rotations 
with corresponding matrices Li, L 2 , L 3 such that the product of 
these rotations is a rotation of axes 

M /zA 

(14) [y]= L,L,Ls y" = N 

W W 

with matrix N = LiL 2 L 3 and U" = U', V" = V. Then W" is 
a point on a line through perpendicular to the x', y' plane and 
W' = eTF", where 

(15) e = 1. 
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The transformation denned by 



(16) 




A"\ 

(y"\ 
\z"/ 



is an orthogonal transformation such that x' = x", y' y", 
z' = e 2 ". Then 

/I 0\ 

(17) R = L*AT = ( 1 01 = N*L = N~ 1 L, 

\0 / 
and 

/I 0\ 

(18) L = NR = LiLtL-AO 1 I. 

\0 / 

Since \N\ = 1, we see that |L| = |fl|, 

(19) = |L|. 

The matrix R is the matrix of a rotation of axes only when 
= 1, 72 = 7. Since TV is the matrix of a rotation of axes and 
R = L*N, we see that K is the matrix of a rotation of axes if and' 
only if L is the matrix of a rotation of axes. Then L is the 
matrix of a rotation of axes if and only if e = |L| = 1. When 
|L| = 1, we see that L = ATR is the product of the matrix of a 
rotation of axes and the matrix R of a reflection of axes. This 
proves the theorem. 

EXERCISES 

Determine whether or not the following matrices are the matrices of 
rotations of axes: 




\ I/ 2 3 6 \ 

(6)7 6 2-3 

/ \3 -6 2/ 





(d) 



1 1 


-1 


1 \ 

V3 



2 


V3 
1 

V2 
-1 

V6 


V3 
1 

V2 
1 
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4. Orthogonal reduction of a real quadratic form. A reaj 
quadratic form in x, y, z is a polynomial 

(20) /(z, y, z) = ax 2 + by* + cz 1 + 2dxy + 2exz + 2gyz, 
where a, 6, c, d, e, g are real numbers. Then 

(21) /(*, y, z) = PAP* 

where P = (x, y, z) and A is the real symmetric matrix 



(22) 



A rotation of axes is a linear transformation 

(23) P* = LQ* 

where L is an orthogonal matrix of determinant 1 and Q 
(x',y',z'}. Then 

(24) P = QL* 

and a rotation of axes replaces f(x, y, z) by 

(25) 0(z', y', z') = aV 2 + &Y 2 + cV 2 + Wx'y' 

+ 2e'x'z' + Ig'y'z' 
where 

(26) g(x', y' y z'} = f(x, y, z) = QL*^LQ* = QfiQ*. 

Hence 

/a' d' 6 r \ 

(27) 5 = ( d' b' g'\ = L*AL. 

Theorem 2. Every real quadratic form PAP* may be reduced to 
a real diagonal form 

(28) ax 72 + /3y' 2 + TZ' 2 

by a rotation of axes where a, /3, 7 are the characteristic roots of A. 
For, by Theorem 16 of Chap. 6, there exists an orthogonal 
matrix L such that 



L *AL<> 



(a 0\ 
/3 0). 
y/ 
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If |Lo| = 1, we take L = L and the orthogonal transformation 
with matrix L is a rotation of axes. If |L | = 1, we may 
change the sign of one column of Lo and replace L by an ortho- 
gonal matrix L such that |L| = |L | = 1. Evidently L*AL 
= LQ*AL Q . Then the rotation of axes with matrix L replaces 



f(x, y, z) by 



y', 



= ax'* 



ILLUSTRATIVE EXAMPLES 



I. Reduce the quadratic form 4x 2 + 3y 2 z 2 12xy + 4x2 8yz to 
diagonal form dix' 2 + d 2 y /2 + d 3 z' 2 by a rotation of axes such that 
di ^ d z ^ da. Give the diagonal form and the equations of rotation. 

Solution 
The matrix of the given form is 



and 



(4 1 4V 

\ 2 -4 -I/ 



4 - a; -6 2 

~6 3 - x -4 
2 -4 _i _ x 

4 - x 6 2 

2(1 - x) 9 + x 

2 4 -1 - x 



Then/(x) = 2[8(1 - a;) - 2(a^ + 9)] - (1 + x)[(4 - x)(a; + 9) - 
12(1 - x)] = -20(x + 1) - (x + 1) - x 2 + (4 - 9 + 12)x + 36-12 
- (x + l)(x 2 - 7x - 44) = (x + l)(x + 4)(x - 11). This yields the 
diagonal form llx /2 y' 2 4z /2 . 

We now solve the system of equations 



-7x - 

-6x 

2x - 



2z 



= 

= 
= 



with matrix A - (11)7. Then 6x - 12?/ - 862; = 0, -20?/ - 40z = 0, 
y == -2, 2x = 4?/ + 12z = -82; + 12z = 42, x = 2^;. The remain- 
ing equation merely yields the check 142 + 122 + 22 = 0. We have 
shown that (2, 2, 1) is a characteristic vector of A (ll)/. 

We next determine a characteristic vector of A + I and thus solve 

5x &y + 22 = 

-6x + 4y - 42 = 

2x - 4y = 0. 
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Then x = 2y, -12y + 4y = 4z = -8?/, z = -2i/. This yields the 
vector (2, 1, 2). The final vector is a characteristic vector of A + 47 
but may be obtained as a vector orthogonal to the vectors already 
determined. Then we have 2x 2y + z = 0, 2x + y 2z = 0, 
3i/ = 3z, y z, 2x = 2y 2 = 2, and we have shown that the columns 
of the matrix 

/ 2 2 l\ 



\ 1 -2 27 

are the required characteristic vectors. We compute 



2 2 

-2 1 

1 -2 




-6 -3 
-3 -6 



= 36 - 9 = 27 



and see that 




is the matrix of a rotation of axes. The equations of this rotation in 
solved form are (x', y 1 ', z') = (x, y, z)L, 



' = 2x - 2y 
/' = 2x + y 
z' - x + 2y 



z 

2z 
2z. 



II. Reduce the quadratic form x 2 y 2 7z 2 + 16x7/ + &cz + Syz 
to diagonal form by a rotation of axes. 



Solution 



The matrix of the form is 



/-I 8 4\ 
(8-1 4 
\ 4 4 -77 



and the characteristic determinant is 

-l-o; 8 4 

8 -1 - x 4 

4 4 -7 - on 





-1 - x 


4 
7 




= 


8 -9 -x 


4 






4 2x + 18 


~~"~ X 






-1 - x 


4 


= (x + 9) 


8 - 


1 4 




4 


2 -7 - z 
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Then f(x) - -\A - xl\ = -(x + 9)[4(16 + 4) - (1 + x). (x + 7 
_ 8)] = (x + 9)(z 2 - 1 - 80) = (x + 9)(z + 9)(z - 9). The diago- 
nal form is 9(z' 2 - 7/' 2 - z' 2 ). 
We first solve the system 

-10z + % + 4z = 
8z - 1% + 4z = 
4z + 4y - 162: - 

with matrix A 9/ and obtain 18i/ = 862, i/ = 2z, z = 42; y = 2z. 
This yields a characteristic vector (2, 2, 1). However, the system with 
matrix A + 97 consists of the three equations 

Sx + 8y + 4z = 
8* + 8?/ + 42 = 
4x + 4y + 2z = 0. 

Moreover, the condition that this vector be orthogonal to (2, 2, 1) is 
2x + 2?/ + z 0. Hence, any choice such as ( 1, 1, 0) is a character- 
istic vector. The remaining characteristic vector satisfies 2x + 2y + z 
= and x y = 0, so that x = y, z = 4y and (1, 1, 4) is a charac- 
teristic vector, 

i2 -1 1 



3 V2 VT8 

2 1__ _l_ 

3 \/2 \/T8 

1 _~4^ 

^ 



L = 



is the matrix of an orthogonal transformation carrying the given 
quadratic form into the diagonal form above. But 



2 -1 1 
2 1 I 
1 -4 






~~~ t 





2 


1 


1 


1 





-4 



- 2(-8 - 1) = -US 



and thus L is not the matrix of a rotation of axes. We therefore change 
the sign of the second column and obtain the following equations of 
rotation: 

3x' = 2x + 2y + z 
\/2?/' = x y 
3 \/2z r = x + y - 4^. 

EXERCISES 

1. Show that (2, 1, 2) is a characteristic vector of the matrix of 
Illustrative Example II corresponding to the root 9, and obtain the 
corresponding equations of rotation. 
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2. Reduce each of the following quadratic forms to diagonal form 
d lX '* + dty' 2 + d&' 2 with di = c? 2 = ds by a rotation of axes. Give 
the diagonal form and the equations of rotation solved for x' , ?/', z' as 
your answer. 

(a) 2x 2 + y 2 - 4.xy - 4yz 

Ans. x'* - 2y /2 + 4z' 2 ; 3tf = 2x + y - 2z; 3y' = x + 2y + 2z; 
3z' = 2x - 2y + z. 

(6) 3.r 2 - 3i/ 2 - 5z 2 - 2.c?/ - 6xz - 6yz 

(c) 3x 2 + y 2 + z 2 - 2xy + 2xz - 2yz 

Ans. 4*' 2 + z' 2 ; \/6x' = 2x - y + z; W = y + z; \/3z' = -x 

- y + z. 

(d) -2.r 2 + 4?/ 2 + 6s 2 + 20;?/ + 6x2 + GT/Z 
(c) 4x 2 + y 2 - 8z 2 + 4xj/ - Axz + 8^/z 

Ans. 5x' 2 + 2// /2 - 10z /2 ; V^' = 2x + y; \^y f = -x + 2y + z; 
\730z' = x - 2?y + 5z. 

(/) 3.c 2 + 3z 2 + 4xy + $xz + 4yz 

6. Quadric surfaces. Suppose that p, g, r, s are real numbers 
and that fo(x, y, z) is a real quadratic form. Then the equation 

(29) f(x, y, z) s / (a:, y, z) + 2(px + qy + rz) + s = 

is the general form of an equation defining a quadric surface. 
By Theorem 2 there exists a rotation of axes which replaces this 
equation by 



(30) <Kz', y', z') = ax' 2 + /ty' 2 + 7z' 2 + 2(px> + ay' + rz') 

+ 5 = 

for real numbers a, /?, 7, p, cr, r, 5 such that a, /3, 7 are not all zero. 
The vector (a, ft, 7) of the characteristic roots of the matrix of 
/o(z, 2A 2) is unique up to a nonzero real factor t that enters when 
we multiply /(x, y,z) =0 by J. Thus we may always take 

(31) a > 0, j8 0, a |8 7. 

In Chap. 5 we studied all the quadric surfaces defined by the 
equation of formula (30) except that where = 7 = and 
p 7* 0, r 7* 0. Then 



0(*', y', *') a a' + + 2(*t/' + r*0 + 5- = 
and we may perform a second rotation of axes defined by 
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(32) 



x" = x' 



This planar rotation of axes replaces 4>(x', y', z') = by 
(33) F(x", y", z") = a (x" + |Y 



-.'- 0. 



a 
The translation of axes whose equations are 




(*A\ r '" x " 4. *L, ? /" = y" 4- H , Z '" = 2" 

(M) xx-f^; y y -r r-sr-\ ' z z 

replaces F(x"j y", z") == by 



y'" = 0. 

This is an equation of a parabolic cylinder with vertex at the 
origin. Then we have proved that every quadric surface is one 
of the surfaces discussed in Chap. 5. 

Let us close this section with a summary of the allowable 
operations on a quadratic equation f(x, y, z) =0 defining a quad- 
ric surface, and of their effect on the quadratic form / (x, t/, z). 

The first operation is that of a translation of axes. It should 
be evident that only the linear and constant terms of f(x, y, z) 
are affected by this operation and that /o(#, y, z) is unaltered. 

The second operation is that of a rotation of axes. This 
replaces fo(x, y, z) by an orthogonally equivalent quadratic form 
having the same characteristic roots. If A is the matrix of 
/ (x, y, z) and L*AL is the matrix of the equivalent form, then 
P(A - dl) = if and only if PLL*(A - dI)L = 0, (PL) 
(L*AL dl) = 0. Thus the vector P is a characteristic vector 
of A corresponding to a root d if and only if PL is a characteristic 
vector of L*AL corresponding to the same root. 

The final operation is that of multiplying /(#, y, z) by a nonzero 
constant k. This multiplies /o(#, j/, z) by k and the matrix A of 
fo(x, y, z) by k. The corresponding characteristic roots are 
multiplied by k and indeed \iL*AL = diag{a, 0, y} then L*(kA)L 
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= diag{ka, fc/S, ky}. Also P(A - dl] = if and only if P(kA 
&d/) = and therefore the characteristic vectors correspond- 
ing to the given quadratic form are unaltered by this change. 

EXERCISE 

Reduce the following equations to simplified form by a rotation of 
axes. Give the center or vertex of the quadric, the type of quadric, the 
semiaxes, and the equations of rotation. 

(a) 2x 2 + 2?/ 2 - z 2 + Sxy - 4x2 - tyz = 2 

(6) 2x 2 + 7/ 2 + 2z 2 + 2xy - 2yz = 1 

(c) 4x 2 + 6y 2 + 4z 2 - 4x2 + 1 = 

(d) 2x 2 + 7/ 2 - 4x?j - 4yz + 12x + 6y + 62 = 1 

(e) x 2 + |/ 2 + 2 2 -- 4yz 4x2 4x?/ = 7 

(/) 5x 2 + 5?/ + 3z 2 - 2x7/ + 2x2 + 2yz + 2x - y = 

(^ 2/ 2 + * 2 - xy + xz + yz - 2x + 2y - 2z + 1 = 

(h) x 2 + 4?/ 2 + 9z 2 - 4xy + 6x2 - 12yz + 4x - 8y + 12z + 4 = 

(f) 16x 2 + 9?/ 2 + 4^ 2 + 24x7j - 16x2 - 12yz + 7x + 2y - 12z = 
( j) 2x 2 + 3y + 42 + 4 = 

(t) 2^ 2 + 5x + 12j/ + 120 + 18 = 

(0 a; 2 + y 2 + 4^ 2 - 2xy - 4x0 + 4^ + 6x + 12y + ISz = 

(m) 2x 2 + 2?/ 2 4^ 2 5x?/ - 2x2 2yz 2x 2y + z = 

(n) 3x 2 + i/ 2 + z 2 + 47/s + 12x + 2y - 20 + 9 = 

(o) 4x 2 + 4^/ 2 + 90 2 + Sxy + 12x0 + 12yz + lOx + y + 42 + 1 =0 

(p) 2x 2 + 4?/2 + 62 + 2y - 4x + 5 = 

(g) 3x 2 + 3?/ 2 + 2 2 - 2xy + 6x - 2y - 22 + 3 = 

6. Plane sections of quadrics. If f(x, y, z) = is an equation 
of a quadric surface and ax + by + cz + d = is an equation of 
a plane, we may rotate axes so that the x* ', t/ 7 plane is parallel to 
the given plane. We may then translate axes so that the given 
plane becomes the plane z 1 = 0. It follows that the plane sec- 
tions of any quadric are the sections by the plane z = of the 
general quadric of formula (29) for a properly selected coordinate 
system. 

As a consequence every plane section of a quadric is a conic 

(35) 2 = 0, ax 2 + by* + 2dxy + 2hx + 2py + r = 0, 

as studied in plane analytic geometry. Such a conic may be 
carried by a rotation of axes leaving z unaltered into a conic 

(36) z = z' = 0, ax' 2 + /3?/ 2 + 2yx' + 2dy' + 6 = 0. 
Apply this rotation to the original quadric and thus carry the 
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quadratic form / (x, y, 2) = ax 2 + by 2 + cz* + 2dxy + 2exz + 
2gyz into 

(37) (*', y', z') s az' 2 + /V 2 + yz'* + 2\x'z' + 2/iyV. 



If the given quadric is an ellipsoid, then <o(#', y', z') > for all 
real values of x', y', z 1 not all zero and this must be true of 
4>(x', y', 0). Hence, formula (37) is our equation of an ellipse 
and we have proved the following: 

Theorem 3. The plane sections of a quadric are conic sections. 
In particular, all plane sections of an ellipsoid are ellipses. 

7. Points of symmetry. We shall close this chapter by apply- 
ing notations of axes to obtain a discussion of the symmetries of 
riondegenerate quadrics, i.e., quadrics that are not planes, pairs 
of distinct planes, point loci, or imaginary loci. It is recom- 
mended that the details of the proofs be omitted in the classroom and 
that only the results be presented. 

We first observe that the nondegenerate quadrics defined by 
equations of the form 

ax* + J3y 2 + 7Z 2 + 5 = 

are symmetric with respect to the origin. We shall call such 
quadrics central quadrics and shall say that the origin is a center. 
If a/37 7* 0, the corresponding surfaces are ellipsoids, hyperbo- 
loids, or quadric cones. Then the origin is the only point of 
symmetry. Indeed, we shall prove the following result: 

Theorem 4. A nondegenerate quadric S is symmetric with 
respect to a point P if and only if S is a central quadric with P as a 
center. If S is an ellipsoid, hyperboloid, or quadric cone, the 
point P is the only point of symmetry. However, if S is a cylinder, 
the line through P parallel to a generating line of S is a line of 
centers oj S. 

For we select a coordinate system with a given point of sym- 
metry P as origin. Then if f(x, y, z) = is an equation of S 
relative to this coordinate system, it must be true that f(x, y, z) 
s= /( x, y, z). Hence f(x, y, z) has no linear terms. A rota- 
tion of axes about P reduces f(x, y, z) to the form ax 2 + /3y 2 
+ yz 2 +5 = and S is a central quadric, P is a center of S. 
Let (, ir], f) be a second point of symmetry, and carry out a 
translation of axes x = x 1 + , y = y r + rj, z = z f + f , which 
moves the origin to (, 77, f). Then f(x, y, z) a(x f + ) 2 + 
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/%' + u) 2 + y(z' + D 2 = *(*', y', z'), *(-*', y',*') = *(,,), 
and a = /fy = 7f = 0. If apy ^ 0, then (, 17, f) = (0, 0, 0) 
and P is the only point of symmetry. Since $ is nondegenerate, 
only one of the numbers a, |8, 7 can be zero, and we can assume 
the coordinate system chosen so that a/3 ^ 0, 7 = 0. Then 
= rj = and the surface S is symmetric with respect to all 
points (0, 0, f), that is, all points on the z axis. But the z axis is 
the line through P parallel to the generating line x = y = of 
the cylinder ax 2 + /3?/ 2 +6 = 0. This proves the theorem. 
8. Planes of symmetry. Define a quadric surface S by a 
quadratic equation /(x, y, z) = 0, and define a plane // by an 
equation \x + py + vz = p, where p ^ and X 2 + ju 2 + *> 2 = 1 
Then a line L normal to // and through a point (x , yo, Zo) of // 
is a line x = + Xi, 27 = 2/o + M^ 2 = 2 + yf. We find the 
points of intersection of L and >S by finding the roots of the 
equation 



(38) <t>(t) s /(a? + H y + ^, z + vt] = 0. 

The distance from // to any point of L is \(X Q + \t) + p(y Q + p,t) 
+ v(z, + vf) - p = (X 2 + M 2 + ^ 2 )^ + (Xx + W Q + vz, - p) = t 
and thus the solutions of <() =0 define not only the points of 
intersection of the normal lines with S but also their distances 
from //. 

A plane // is a plane of symmetry of a quadric surface S if 
every normal to // either does not cut S or cuts $ in two points 
on opposite sides of // and at the same distance from //. While 
this definition formally includes a plane having the property that 
all normals either do not cut the surface or lie wholly on it, we 
shall exclude such planes. Thus the planes z = k are not 
regarded as planes of symmetry of the cylinders /(x, y) = 0. A 
plane of symmetry then defines a set of normal lines and corre- 
sponding pairs of points of intersection Pi and P 2 whose corre- 
sponding distances Zi and 2 from // must have the property 
t\ + h = 0. Then necessarily the coefficient of t in 4>(t) is zero 
for all (XQ, 2/o, Zo) in 77. 

If f X) f y , f z are the partial derivatives of f(x, y, z) with respect 
to x, y, z, the coefficient of t in formula (38) is 



<iz , 2/o, Zo s rco, 2/o, *o ttv&o, j/o, Z v,xo, y , ZQ). 
We thus seek to determine all values of p, X, /*, v such that 
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o, 2/0, ZG) = for (x Qj 7/0, 2 ) on the plane H defined by 
\x + vy + vz p = 0. 
The central quadric 

(39) ax* + fty* + T z 2 + 5 = 0, 

defined for real a, /J, 7, 5 such that a/?7 y 0, has the corresponding 
function 

(40) 0i(zo, 2/o, Zo) =5 aX^o + /fyyo + yvZo. 

If X ^ 0, then (pX" 1 , 0, 0) is in // and therefore a\x Q ap = 0. 
Then p = 0. Similarly p = if /z =^ or v ^ 0. Since X 2 + /z 2 + 
*/ 2 = 1, we have proved that H is a plane through the center of S. 

If X = M = 0, then v 9^ and 77 is the plane 2 = 0. If 
X = and M ^ 0, then T/ O = jjr l vzo and (7 P)PZQ = for all 
values of 0:0 and 2 . Then (7 ff)v 0. If /S is a sphere, any 
plane through the center of S is a plane of symmetry. If S is 
not a sphere, we may take 7 =^ /? and see that *> = 0; H is the 
plane y = 0. 

There remains the case X -^ 0, #o = X"" 1 ^^ + ^o)(|8 a) 
M2/o + (7 a) 720 = for all T/ O and z . Then (ft )/* = 
(7 a)y = 0. When S is not a surface of revolution, ft 7^ a, 
7 T^ a so that /z = p = and # is the plane x = 0. If S is a 
surface of revolution defined by a = ^ 7, we have y = 0; H 
is a plane X# + /zy = through the axis of revolution x = y = 0. 
But all such planes are planes of symmetry. We have proved the 
following: 

Theorem 6. A plane is a plane of symmetry of a sphere if and 
only if it is a plane through the center of the sphere. If S is a central 
quadric of revolution, its planes of symmetry are the planes through 
the axis of revolution and the plane through the center perpendicular 
to the axis of revolution. The only planes of symmetry of a central 
quadric ax 2 + /3y 2 + 7Z 2 +5 = 0, which is not a surface of 
revolution, are the planes x = 0, y = 0, and z = 0. 

The equation 

(41) ax 2 + fty* = 252 

defines a paraboloid for all nonzero real numbers a, ft, 5. Then 

(42) 0i(z , 2/o, Zo) = aXzo + 0M2/o S". 

If 7 T* 0, then (0, 0, pv~ 1 } is on H and X (0, 0, pz/- 1 ) = 6v = 0. 
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But 5 7* 0, v 7* 0, which is a contradiction. Hence, /* = and 
Xa: + M2/ = P is an equation of //, <i = aXz + /?M2/o. If X -^ 0, 
then (X- 1 ?, 0, 0) is on # and <t>i(\~ l p, 0, 0) = ap = 0, p = 0. 
If X = 0, then ^ and (0, /r- 1 ??, 0) is on H, ^(O, M" 1 ?, 0) 
= )8p = 0, and p = 0. We have thus shown that every plane 
of symmetry of S is a plane X# + /xy = 0. If a = 0, so that S 
is a surface of revolution, all such planes are planes of symmetry. 
But if a 7* j8, we have \X Q = M2/o, <i = (/S <*)M2/o. When 
X ^ 0, <i = (/3 a)jLty = for all values of j/ and /* = 0, H is 
the plane z = 0. Otherwise, X = and H is the plane T/ = 0. 
We have proved the following: 

Theorem 6. A paraboloid ax 2 + /3y 2 = 25z has the planes 
x = 0, y = as planes of symmetry. These are its only planes of 
symmetry unless the paraboloid is a surface of revolution. In this 
latter case, the planes of symmetry are all the planes through the 
axis of revolution. 

The only remaining nondegenerate quadrics are the cylinders. 
The first case to be studied is that of a cylinder 

(43) ax 2 + /ty 2 = 5 

where a/38 7* 0. We have already seen that a plane z = k is 
not a plane of symmetry and therefore all planes of symmetry 
are planes \x + ny + vz = p, where X ^ or /z ^ 0. Now 
4>i(xo, j/o, zo) = oi\x Q + foy Q and (X~ ! p, 0, 0) or (0, pr l p, 0) is 
on H, <t>i(x Q , 7/0, Zo) = otp or Pp at these points and therefore 
p = 0. If X = 0, then M ^ 0, 1/0= M" 1 ^, aX# /^z = 
for all values of #o, ZQ, X = v = and // is the plane y = 0. If 
X 7* 0, then X = \~ l (nyo + PZ Q ), <f>i(x Q , y Q , z ) = (/3 a)/x^ 
avz Q = for all 2/0 and z , (/? ~~ <*)M = ^ = 0. When = a, 
the surface $ is a circular cylinder and all planes through its axis 
x = y = are planes of symmetry. In cases 7* a, then 
ju = v = and // is the plane a; = 0. We have proved the fol- 
lowing result: 

Theorem 7. The planes of symmetry of a circular cylinder are 
the planes through its axis. If a noncomposite cylinder ax 2 + 
/3y 2 = 5 is not a circular cylinder, the only planes of symmetry are 
the planes x = and y = 0. 

The only remaining case is that of the parabolic cylinder 
x 2 = Sy, where 5^0. As before, every plane of symmetry is a 
plane \x + py + vz = p, where X ^ or /z 5^ 0. Now <i>i(xo, 
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yo, Z Q ) X#o 5/x. If /z 7^ 0, then (0, ^p, 0) is on H, Sp, = 
contrary to hypothesis. Hence, M = 0, X ^ 0. Then (X-*p, 0, 
0) is on H, <t>\ = Xzo = p = 0. It follows that # is a plane 
Xx = fjiz, ( \~ I VZQ, y Q , ZQ) is on H for all y , z , <i = 0, and H is 
the plane x = 0. We have proved the following: 

Theorem 8. The only plane of symmetry of a parabolic cylinder 
x 2 = 5y is the plane x = 0. 

9. Lines of symmetry. If L is a line of symmetry of a quadric 
surface S, we may select a coordinate system such that L is the 
z axis. Let an equation of S relative to this coordinate system 
be f(x, y, z) = 0. Then f(x, y, z) = f(-x, y, z) = /(z, -y, 2), 
so that /(x, y, 2) involves only even powers of x and y. Then 



(44) f(x, y, z) = ax* + /3?/ 2 + yz* + 2Sz + c = 0. 

Suppose first that S is a central quadric with three nonzero 
characteristic roots. Then aj3y ^ 0, and we may select the 
coordinate system so that 6=0. The given line is then a line 
of intersection of two distinct planes of symmetry, i.e., the planes 
x = and y = 0. When a, /3, 7 are distinct, the surface S has 
only three planes of symmetry and thus the only possible lines 
of symmetry are the three intersections of pairs of such planes. 
These lines are called the principal axes of the quadric. When 
iS is a sphere, all lines through the center are lines of symmetry. 
However, if S is a quadric of revolution, the intersection of two 
distinct planes through the axis of revolution is the axis of revolu- 
tion and the intersection of a plane through the axis of revolution 
and the unique plane perpendicular to the axis of revolution is a 
line through the center perpendicular to the axis of revolution. 
These lines are all lines of symmetry. We have proved the 
following : 

Theorem 9. The lines of symmetry of a sphere are all the lines 
through its center. If S is a noncylindrical central quadric of 
revolution, its lines of symmetry are the axis of revolution and all 
lines through the center of S and perpendicular to the axis of revolu- 
tion. The lines of symmetry of a noncylindrical central quadric, 
not a surface of revolution, are its three principal axes. 

We next assume that S is a paraboloid. Then necessarily 
7 = 0, a ? 0, j8 7* in formula (44), and the line L of symmetry 
is the intersection of two planes of symmetry of the surface. By 
Theorem 6, if S is not a surface of revolution, L is the unique line 
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of intersection of its two planes of symmetry. If S is a surface 
of revolution, all planes of symmetry intersect on the axis of 
revolution and L is this axis. 

Theorem 10. A paraboloid has a single line of symmetry. 

If S is a cylinder, then apy = 0. If 7 = ft = 0, the equation 
is ax 2 + 25z + e = 0, and the line L lies on the plane of sym- 
metry x = 0. But this is the case of a parabolic cylinder and 
there is only one such plane. Hence, all lines of symmetry are 
lines of the plane x = 0. If the line L is defined by equations 
x = 0, Ay + \iz p, then the rotation of axes defined by 

x' = x, y' = X?/ + nz, z 1 = py + Xz 
has as solved form 

x = x' y y = X?/ - MZ', z = M' + Xz' 
and replaces the given equation by 

ax' 2 + 2d(fjiy' + Xz') + c = 0. 

This rotation replaces the equation of L by x' = 0, y' = p. A 
translation of axes x" = x', y" = y r + p, z" = z' yields 

ax" 2 + 26[n(y" + p) + Xz"] + e = 0. 

This must be unaltered by the replacement of y" by y". 
Then JJL = 0, and the line is given by y = k p. However, 
every line x = 0, y = k is a line of symmetry of ax 2 + 25z + e 
= 0, since x = y = is such a line, and the translation of axes 
x' x, y f = y + p does not alter the given equation. Inter- 
changing the roles of y and z, we have the following results: 

Theorem 11. The lines of symmetry of the parabolic cylinder 
x 2 = Sy are all the lines of intersection of the plane x = with 
planes z = k perpendicular to the lines on the cylinder. 

If 7 = a = and ft ^ 0, we also obtain a parabolic cylinder. 
Thus the only remaining cases are 7 = 6 = 0, /(x, y, z) = ax 2 
+ f3y 2 + 6, and 7^0, f(x, y, z) = ax 2 + yz 2 + 2dz + e, or 
f(x, y, z) = /3y 2 + 7Z 2 + 25z + e. In any case we may reduce 
the equation to ax 2 + fiy 2 + ^ = by a translation and rotation 
of axes and the line of symmetry becomes the line of intersection 
of two planes of symmetry. But, by Theorem 7, all planes of 
symmetry of a nonparabolic quadric cylinder intersect in a single 
line and we have proved the following: 

Theorem 12. A nonparabolic cylinder has a single line of 
symmetry. 



CHAPTER 8 
SPHERICAL COORDINATES 

1. Azimuth and elevation. It is sometimes convenient to 
locate points in space by three-dimensional vectors of coordinates 
called spherical coordinates. These coordinates have a position 
in space geometry like that of polar coordinates in plane geometry. 
They have the advantage over rectangular coordinates of being 
optically measureable from a fixed origin O of coordinates. 

The first spherical coordinate is the range r. This is the 
distance 



(1) r = \OP\ = Vx 2 + y* + z 2 

from the origin O to the arbitrary point P in space, and we have 
given a formula expressing r in terms of the rectangular coordi- 
nates of P. We note that r g and that r = only if P is at O, 
There are several physicogeometrical methods for the approxi- 
mate measurement of r and, in particular, the measurement of 
the distance of a physical object from an origin of measurement 
is possible by radar devices. 

The remaining two coordinates are angles a, e measured (in 
radians) as in Fig. 16. The set of numbers (r, a, e) is then a 
coordinate vector of three real numbers. If the angles are 
measured in degrees, they must be converted into radians. 
Actual measurement of a and e may be achieved by a very simple 
and accurate instrument called a surveyor's transit. 

The spherical coordinate angles a, e may be defined most 
easily relative to a given rectangular coordinate system. We 
assume that P = (x, y, z) is an arbitrary point in space and con- 
struct a half plane which contains P and has the z axis as its 
edge. In Fig. 10, this half plane intersects the x, y plane in a 
ray OA. Look down from the positive z direction on the x, y 
plane and so define a clockwise direction for angular measurement 
in the x, y plane. The angle a is then measured in a clockwise 
direction from the positive y axis to the ray OA. Evidently a. 

123 
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may be restricted to lie in the interval 
(2) ^ a < 2ir 

radians. The angle a is called the azimuth of P. 

z 



[Chap. 8 




FIG." 16. 



The angle e is measured from the ray OA to the ray OP, and 
we may assume that 



(3) 



2 = 



Here c > for points with z > and c ^ otherwise. We call 
c the elevation of P. Note that the x, y plane is the plane whose 
spherical coordinate equation is e = 0. 

It is sometimes convenient to use an angle f called the zenith 
of P rather than the elevation angle . This angle is measured 
from the positive z axis to the ray OP and f = 7r/2 e. Then 



^ 



TT. 
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Rectangular and spherical coordinates are related by the 
formulas 

(4) x = r sin a. cos , y = r cos a cos , z = r sin e. 

The last of these formulas comes from the fact that in Fig. 16 

> 
the distance BP = z and |OP| = r, so that z = r sin e. Then 

= r cos and y = CC = \OB\ cos a, and s = CB = 



sin a. This verifies formula (4). 

ORAL EXERCISES 

1. The following vectors are the rectangular coordinates of points. 
Give their spherical coordinates. 

(a) (0, 1, 0) (d) (0, 0, 5) ( ff ) (V2, V6, 2 x/2) 

(6) (1, 0, 0) (e) (1, 1, 1) (A) (V6, V6, -2) _ 

(c) (0, 0, 1) (/) (1, 1, 0) (i) (- V3, 1, 2 V3) 

2. The following vectors are the spherical coordinates (r, a, e) of 
corresponding points. Give their rectangular coordinates. 

(a) (0,a,e) ^( 2 'I'l) (ff) 0' T ~ l) 

<&) (l.O.e) (e) 3,. (A) 

(0 



3. What geometrical object is described by the equation a = 0? 
By the equation a(a TT) = 0? 

2. The angle between two vectors. If P and P are two points 
in space distinct from the origin 0, we have already given a 
formula for the cosine of the angle 6 = ZPOP in terms of the 
direction cosines of the rays OP and OP Q . We shall now derive 
a formula in terms of the spherical coordinates of P and P Q . 

If P = (r, a, e) and P = (TO, a , e ), the unit vector on the 
ray OP is U = (1, a, e) and that on the ray OP is f/ = (1, , e ). 
The rectangular coordinates of [/ are (sin a cos e, cos a cos e, 
sin e) and those of UQ are (sin ao cos , cos ao cos e , sin o), so 
that cos 6 = (sin a sin o + cos a cos ao) cos e cos o + sin e sin en. 
Then the required formula is 

(5) cos 6 = cos c cos e cos (a a ) + sin c sin CQ. 
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If is an angle so small that its sine is approximately equal to 
6 (measured in radians), we may replace formula (5) by a much 
simpler approximation formula. We shall use the notation 

a == b 

(read a is approximately equal to 6) throughout this chapter and 
the degree of accuracy of the approximation will of course depend 
on the actual approximations made and needs to be investigated 
in all specific cases. 

Formula (5) implies that 

1 cos 6 = cos e cos eo[l cos (a <*o)] + 1 
(6) cos cos sin e sin eo 

= COS COS [1 COS (a ao)] + 1 COS (e e ). 

If </> is any angle, the formula 1 cos = 2 sin 2 0/2 is a well- 
known formula of elementary trigonometry. We apply it to 
formula (6) to obtain 

. 9 6 . 9 fa a \ , ./ e \ 

sin 2 s = sin 2 1 = 1 cos c cos + sin 2 I -= 1, 

and therefore have the exact formula 



/*7\ / 2 ( a ~~ a O\ , o A CQ\ 

(7) sin ^ = <i/sm 2 1 ^ 1 cos e cos e + sin 2 I ^ 1 

for acute angles 6. When 6 is small, the angles a a and 
e e are small and cos e = cos e , sin ^(a ) i( ~ o), 
sin ^(e eo) = -2-(e CQ). Then formula (7) becomes 



(8) 

EXERCISES 

1. Compute in milliradians (thousands of radians) in the following 
cases by the use of the exact formula (7) and five-place tables of loga- 
rithms. Give your answer correct to the nearest tenth of a milliradian. 

(a) a = 1727', a Q = 1842', c = 15, e = 16 

(b) a = 2432 X , a = 24, e = 72, c = W4Q' 

(c) a = 19217 / , a = 191, e = -18, e = -19 

2. Compute in milliradians correct to the nearest tenth of a milli- 
radian by the use of formula (8), and compare the results with those of 
Exercise 1. 
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3. Parallax. The changes in spherical coordinates that are 
the results of a translation of axes are called parallax corrections. 
There are many physical situations where measurements of 
range, azimuth, and elevation are made by instruments located 
at a point O, but are required relative to a parallel coordinate 
system with origin at a nearby point 0'. In such cases, formulas 
for parallax corrections are needed. 

Let us suppose that the #, y, z coordinates of 0' are (a, b, c) 
and that the corresponding spherical coordinates are d, , e . 
Thus 



a 



(9) d = \/a* + 6 2 + c 2 , sin e = y tan o = T 

Suppose then that P (r, a, e) has rectangular coordinates 
x, y, z, and define 

/m\ ^ d a ax + by + cz 

(10) d = -> cos = T > 

r dr 

so that 9 is the angle between the ray 00' and the ray OP. Then 

(11) c? cos = (a sin a + 6 cos a) cos + c sin e, 

and d cos may be computed when a and e are given. The 
numbers a, 6, c may be regarded as fixed, and it is possible to 
construct a table of corresponding values of d cos in which the 
table is entered with a and e and d cos 6 is read from the table. 
We may, of course, apply formula (7) to compute d and then com- 
pute d cos 0. 

When many parallax corrections to the same point 0' must be 
made, it is worth while to construct a table as described above 
and to translate it into a contour graph. This graph will consist 
of a number of curves. Each curve will consist of points whose 
x coordinates are a radians and y coordinates are e radians, and 
d cos will be the same for all points on the same curve. If 
enough curves are drawn, it is possible to read off d cos directly 
from the graph by finding the curve on which a point with pre- 
scribed coordinates lies. Linear interpolation between curves 
is sufficiently accurate, of course, when many curves are drawn. 
It is usually possible to draw the curves by plotting a relatively 
small number of points on each curve. 

Let us suppose now that d cos can be quickly determined by 
a table or a contour graph whenever a and e are given. We 
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require the range r', the azimuth a', and the elevation e' of the 
point P relative to the coordinate system obtained by a transla- 
tion of axes carrying the origin to 0'. Then the parallax correc- 
tions are defined to be 

(12) A(r) = r' - r, A (a) = a' - a, A(e) = e' - e, 

and we propose to find formulas for the computation of A(r), 
A(), A(e). 

We use the translation formulas x f = x a, y' = y &, 
2' = z c, and compute r' 2 = a/ 2 + 2/ 2 + z' 2 = x 2 + y 2 + z 2 
- 2(az + by + cz) + a 2 + 6 2 + c 2 = r 2 - 2dr cos + d 2 by 
formula (10). Then 

(13) r' 2 = r 2 (l - 26 cos + 6 2 ) = r 2 [(l - 6 cos 0) 2 

+ (5 sin 0) 2 ]. 

The assumption that d is small compared with r implies that 
1 > 5 cos 0, 



(14) ^Kl- 

Expand this expression by the binomial theorem which states that 



V-*- I V) A I fl/iy | fy V/ | 

and which yields a convergent series if |6| < 1 and n = ?. Then 

/-, r-\ t /-, /.N f\ , 5 2 sin 2 

(15) / = r(l- 5 cos0)|_l + 2(1 ^ gcQsg)2 

^^ Q + . . .] 

8(1 - d cos 0) 4 ^ J 

and therefore dr = d yields 

(16) / = r - d cos + 0/ d2s ^ n2 ^ 
v y 2(r d cos 0) 

The parallax correction 

(17) A(r) = d cos 

is usually quite adequate and may be computed as was indicated 
above. The more accurate correction 

/io\ A/ \ ,1 * i d 2 d 2 cos 2 

(18) A(r) = -d cos 6 + ^r^T^g) 
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may also be computed easily from given values of d cos 0, but the 
difficulty in making exceedingly accurate measurements of r 
usually make this more refined correction meaningless. 

We next observe that r' sin e' = z r = z c = r sin e c. 
This gives 

/im . , r sin e c 

(19) sine' = - f , 

and 

/OA\ / (r r 7 ) sin e c 

(20) sin e' sin e = ' 

Then we use a well-known formula of trigonometry to replace 
formula (20) by 

(21) 2 sin ^ cos --^ = (r - r') sin e - c 

& & r 

If a, 6, c are small compared with r; the values of A(e) will be 
small, and we may use the approximations 



sin^^>, 



(22) 



cos s = cos I H I = cos 





and thus obtain 

(23) 



r cos 6 



This correction is proportional to r' and there is little loss of 
accuracy if the approximation 

(24) 1 d cos sine -c 

v ' v r cos e 

is used. The answer is the number of radians in A(e) and is 
converted to milliradians by multiplication by 1,000. 

As in the case of the computation of A(r), the computation of 
many values of A(r) for fixed a, 6, c may be achieved best by a 
graph of curves giving equal values of A(e) f or r = 1 and a and e 
as variables. The computation of A(r) is then a matter of reading 
a result from a graph and dividing the result by r. 

To derive a formula for A (a), \ve note that 

, . x' r sin a cos 6 a 

tan a = : = =-" 

y r cos a cos e 6 
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Then 

tan a 1 tan a = 

(r sin a cos e a) cos a (r cos a cos e b) sin a 
cos a(r cos a cos e b) 

This formula reduces to 

sin (a' a) a cos a + 6 sin a 



cos a cos a 1 cos a(r cos a. cos 6) 
It follows that 



/oc v / / \ / 

(25) Sin (a a) = COS a: ---------- 7- 

r cos a cos e 6 

We make the approximations sin (of a) = a' a = A (a), 
cos a/ == cos a, to obtain 

(26) A( 



- , 
r cos e b sec a 

But then 

, N a cos a + b sin a 



r cos 



+ (a cos a + b sin a) ( - 7 J 

\r cos e b sec a r cos e/ 

_._ a cos a + b sin a . / 6 sec a \ 
"" r cos e \r cos e 6 sec a) 

___ a cos a + 6 sin a / 6 \ 

r cos e \ r cos e cos a. &/ 

If b/r is small, we may use the approximation 
(27) A(a) = - 



r cos 



and may compute A (a) by a process similar to that used in the 
computation of A(c). Otherwise, it will be necessary to multiply 
the result by the factor 

(28) 



cos a cos 
where ft = 6/r. 
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EXERCISES 

1. Let a = 2 yards, 6 = 3 yards, and c 6 yards. Construct a 
table of values of d cos to tenths of a yard for a = 70, 71, 72, 73, 
74 and e = 15, 16, 17, 18, 19. 

2. Compute the values of A(r) by the use of formula (17) for the fol- 
lowing values of (r, a, e), where we are giving a and e in degrees and r in 
yards. 



(a) (1,000, 70, 18) 
(6) (972, 71, 17) 
(c) (500, 74, 15) 



(d) (300, 73, 15) 

(e) (100, 72, 19) 
(/) (75, 71, 18) 



3. Compute A(r) to tenths of a yard by the use of formula (18), and 
compare results. 

4. Compute A(c) in milliradians by the use of formula (24) and com- 
pare the result with that obtained by replacing r by r'. 

5. Compare A (a) by the use of formula (27) for parts (a) and (6) of 
Exercise 2. Compute the correction factor of formula (28), and com- 
pare results. 

4. Other spherical coordinates. The spherical coordinates 
a, e described in Sec. 1 are defined for every point P in terms of a 
half plane through P and the z axis. It should seem natural then 
to define sets of corresponding angles relative to the remaining 
coordinate axes. 




FIG. 17. 

We first construct a half plane through P and the y axis as in 
Fig. 17. Let this half plane cut the x, z plane in a ray OB and 
define the clockwise direction on the x, z plane as viewed from 
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the positive y direction. Then we measure the y azimuth angle 
/3 in this clockwise direction from the positive x axis to the ray 
OB and 

(29) ^ ft < 27r 

radians. It follows that ^ /3 g TT f or points having positive 
or zero elevation and TT < /3 < 2?r for points having negative 
elevation. 

An angle 6 may also be measured from the positive y axis to 
the ray OP. Then 

(30) o ^ e g TT, 

and the angle is a zenith angle rather than an elevation angle. 
We shall call the y zenith of P. It is called the angle off the nose 
in aerial navigation. 

The set of spherical coordinates just defined is determined by 
the relations 

(31) x r sin cos |8, y = r cos 0, z = r sin sin /3. 

We leave the verification of this formula to the reader. It may 
be combined with formula (4) to yield 

(32) cos 6 = cos a cos c. tan /3 = - -- 

; sin a 

These last expressions may be used to determine 6 when a and e 
are given and will also determine )8 up to a quadrant determina- 
tion that may be made by a consideration of signs. 

A third set of spherical coordinates may be obtained by passing 
a half plane through P and the x axis as in Fig. 18. Let the ray 
OD be the intersection of this half plane and the y, z plane, and 
define the clockwise direction on the y, z plane when it is viewed 
from the positive x direction. Then we measure the x azimuth 
angle 7 from the positive z axis to the ray OD. We similarly 
measure the x zenith angle \!/ from the positive x axis to OQ. The 
angles just defined satisfy the inequalities 



(33) ^ 7 < 2;r, TT. 
They also satisfy the relations 

(34) x = r cos $< y = r sin \f/ cos 7, 2 = r sin ^ sin 7, 
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and 7 and ^ may be determined when a and are given by the 
formulas 

/oi-\ i A tan e 

(35) cos ^ = sm a cos e, tan 7 = 

v ' ' COS a 

and a consideration of signs. 

z 




FIG. 18. 

6. The matrices of planar rotations. We have called a rota- 
tion of axes in which one coordinate axis remains unaltered a 
planar rotation. If the unaltered axis is the z axis, we call the 

y 




FIG. 19. 

rotation a yaw. We call the rotation a pitch if the unaltered 
axis is the x axis and a roll if the unaltered axis is the y axis. 
We shall measure all three of these angles in a clockwise positive 
direction as viewed from the positive direction on the unaltered 
axis on the coordinate plane perpendicular to this axis. 
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In Fig. 19 a yaw through an angle f has been performed that 
carries x, y, z to x f , y', z f . It should be evident that e is unaltered 
and that 

a = a' + f, = e'. 

We will normally be given a', e' and will desire a, e. Then we 
have 

RULE I. To remove the effect of a yaw through f increase the 
azimuth a' by f and leave e' unaltered. 

The result above is evident geometrically from Fig. 19. It 
should also be noted that the equations of rotation are 



(36) 



= F 




(cos f sin f 
sin f cos f 
001; 



where we apply this result to the unit vector to see that 

cos 



rsm a cos c cos f sin f 

(37) 1 cos a cos e ) = I sin f cos f 
sin e / \ 




sin e 



This formula implies that sin e = sin e' and so e == e', cos e = 
cos e'. Hence, sin a = sin a' cos f + cos a' sin f = sin (a +f), 
cos a = sin f sin a 1 + cos f cos a = cos (a' + f ), and so a = 



a 




FIG. 20. 

In a similar fashion we use Fig. 20 to measure a roll angle rj 
about the y axis. Then /3 = /3' + 17, 0' = 0, and we have the fol- 
lowing rule: 
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RULE II. To remove the effect of a roll through our angle r? 
increase the y azimuth /?' by 77, and leave the y zenith 6' unaltered. 
As before, the equations of rotation are 



(38) (y = 7, 



where we use 




r cos 77 sin 
= |01 

n 77 cos 



(sin cos )8\ /cos 77 sin rj\ /sin 0' cos /3'^ 
cos 1 = f 1 1 f cos 0' 

sin sin /?/ \sin 77 cos 777 \sin 0' sin $'/ 

to verify that = 0', cos j8 = cos 77 cos 0' sin 77 sin 0' = 
cos (77 + j8 ; ), sin /?' = sin 77 cos ' + cos 77 sin /?' = sin (77 + ft'), 

We finally use Fig. 21 to measure a pitch angle such that 



y' 




FIG. 21. 



7 = T 7 + , ^ = ^ ; . Then we have the following rule: 

RULE III. To remove the effect of a pitch through an angle 

increase the x azimuth y' by and leave the x zenith \j/ unaltered. 
The equations of rotation in this case are 



(40) 
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where we verify the correctness of our formulas by the computa- 
tion 

/ cos ^ \ /I \ / cos t 

(41) I sin \l/ cos 7 1 = I cos sin M I sin $' cos 7' 
\sin \[/ sin 77 \0 sin cos / \sin ^ sin 

This formula yields \l/ = ^', cos 7 = cos cos 7' sin $ sin 7' 
= cos ( + 7')) si n 7 = sin cos 7' + cos $ sin 7' = sin ( + 7'). 
6. Rotations as products of planar rotations. In Sec. 3 of 
Chap. 7 we gave a geometrical argument showing that every 
rotation of axes may be expressed as a product of three planar 
rotations. There are many expressions of a rotation of axes in 
terms of three planar rotations, and in particular we may write 





(42) ?/ = L U' , L = MI M2 MS = YfPtR,. 



Then we have expressed the rotation matrix L as a product of a 
yaw matrix, a pitch matrix, and a roll matrix. There actually 
exist physical instruments for the measurement of f , , 77. We 
shall be interested here primarily in observing how the elements 
of the matrix L are related to the angles , 77, f . 
We first form the product 

(cos 77 sin 77 

sin sin 77 cos sin cos 77 
cos sin 77 sin cos cos Tjy 
The definition of formula (42) implies that 
(44) L = 



cos f 


sin f 


\ 


(cos 77 





sin 77 


sin f 


cos f 





sin 





sin 


TJ 


cos 


sin ^ 


\ cos 


y 








I/ 


cos 





sin 


f] 


sin 


cos I 


COS 





y \/ vy 

and it follows that 

\i = cos f cos 77 sin f sin sin 77, 

X 2 = sin f cos , 

Xa = cos f sin 77 + sin f sin cos 77, 
(45) m = sin f cos 77 cos f sin sin 77, 

M2 = cos f cos , 

Ma = sin 77 sin cos f sin cos 77, 

PI = cos sin 77, j> 2 = sin , v$ = cos cos 17. 
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When L is given, the angles , 77, f may be determined, apart 
from quadrant, by the use of the formulas 

(46) sin = */ 2 , tan 77 = > tan f = 

J>3 M2 

In the case where all rotation angles lie between ir/2 and 7r/2, 
formula (46) completely determines , 17, f . 

It should be clear that if we write L as a product of three 
planar rotations in a different order, then the yaw pitch and roll 
angles are different; for example, if L = Y^P^R,, = P\RnY v , then, 
in general, it will not be true that f = v. 

7. Stabilization of coordinates. The problem of stabilizing 
coordinates arises in physical situations where it is required to 
find the azimuths a and elevations e of many points P relative to a 
fixed coordinate system under conditions where the azimuth a' 
and elevation e' of each point P relative to a rotated coordinate 
system, as well as the corresponding rotation, are measurable. 
It is customary to call a', e' the unstabilized coordinates of P and 
a, e the stabilized coordinates of P and we are thinking of a case 
where the rotation varies as P changes. 

Let us assume that the rotation of axes is given by an equation 

sin a cos sin a! cos e' 




cos a cos c I = FfP^/2, I cos a! cos e r 
sin / \ sin c' 

where we have written the effect on the general unit vector of a 
rotation which is the product of a yaw, a pitch, and a roll. Since 





(sin (a f ) cos e\ 
cos (a - f) cos c 1 = 
sin e / 

The required coordinates may then be obtained by the following 
procedure : 

a. Convert a', e' to j8', 0' by the use of formula (32). 

b. Compute 0" = + ij, 0" = 0'. 



138 SOLID ANALYTIC GEOMETRY [Chap. 8 

c. Convert 0", 6" to 7", $" by the use of the formulas 

cos \f/ = sin 6 cos 0, tan 7 = tan sin /3. 

d. Compute 7'" = 7" + , ^'" = *" 

e. Convert 7'", i//" to a!" , e'" by the use of the formulas 

sin e = sin \p sin 7, tan a = cot ^ sec 7. 

/. Then a = f + a'", e = e'". 

The formulas given in (c) and (d) are obtained by combining 
formulas (4), (31), and (35). The procedure above evidently 
requires a considerable amount of computation. 

A more direct formula is obtained by the use of formula (42) 
and matrix multiplication. It yields 

(50) sin (a f ) cos e = cos 77 sin a! cos e' sin 17 sin e' 

(51) cos (a f) cos e = sin sin rj sin c/ cos e 7 

+ cos cos a' cos e sin cos rj sin e' 

(52) sin = cos sin 77 sin a' cos e' + sin cos a' cos ' 

+ cos cos rj sin e'. 

The computation of a and e by the use of these formulas should 
require less computation than the procedure outlined above. 
However, there is a graphical (gnomonic chart) method for using 
the procedure above that makes it a very rapid one. 

If the rotation angles , r?, f are very small, we can use the 
approximations 

sin T/ == 77, sin = , cos rj == cos = 1, 

sin sin rj = 0. 
Then formulas (50), (51), and (52) yield 

(53) tan ( - f ) * MlL^^^)^ 
v ' cos a cos e f sin c' 

(54) sin e = (17 sin c/ + cos a/) cos e' + sin e'. 

We may then write sin e sin e' = 2 sin = cos = = 

z z 

(e ; ) cos c' and thus obtain the further approximation 

(55) e = e' + T] sin a' + cos a' 

for a rapid computation of c. However, even these approxima- 
tions are not so rapid as the graphical method we have referred 
to above. 
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8. Gnomonic projections. The gnomonic projection of a point 
P on the plane G whose rectangular coordinate equation is 
y = a > is the point Q of G whose azimuth and elevation are 
those of P. The plane G is tangent to the sphere of radius a at 
the point of zero azimuth and elevation, and Q is the point of 
intersection of the line x = t sin a cos e, y = t cos a cos c, z = 
t sin e with G. Then the coordinates of Q are 

(56) x = a tan a, y = a, z a tan sec a. 

The gnomonic projection o/ cm?/ surface S on G is the locus of 
all the points of G that are the gnomonic projections of the 
points of S. Consider in particular the cone consisting of all 
points in space having a fixed elevation e. The projections of 
these points are given by formula (56) for e fixed, and they 
satisfy z 2 = a 2 tan 2 e sec 2 a = tan 2 e(a 2 tan 2 a + a 2 ) = tan 2 e(x 2 
+ a 2 ). Then the gnomonic projection of the cone of points 
having elevation e = is the hyperbola 

(57) - * = a, y = a. 



The points of zero elevation project into the line z = 0, y = a. 

The points having a fixed azimuth a lie on a half plane that is 
a part of the plane x = y tan a. Then the locus of the gnomonic 
projections of these points is the vertical line 

(58) x a tan a, y = a. 

To obtain the gnomonic projections of those points which have 
a constant y zenith angle 0, we use formula (31) and put y a. 
Then 

(59) x = a tan cos /}, ?/ = a, 2 = a tan sin /3, 
so that 

(60) x 1 + z 2 = a 2 tan 2 0, y = a. 

Hence, the projection of the cone of all points having y zenith 
equal to lie on a circle of G with center at (0, a, 0) and radius 
a tan 0. The points having y azimuth ft all lie on the plane 
z = x tan j8, and this plane cuts the plane G in the line 

(61) z = x tan j8, y = a. 
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Here /3 is the angle from the coordinate axis z = 0, y = a in G 
to this line as measured in the usual counterclockwise direction. 
Note that G is parallel to the zx plane and that we are viewing 
G from the negative y direction. 

Let us finally consider the x zenith and azimuth. We use 
formula (34) to write 

(62) x = a cot \l/ sec 7, y = a, z = a tan 7. 

Then the gnomonic projections of all points whose x zenith is ^ 
satisfy the equation x 2 = a- cot 2 \[/ sec 2 7 = cot 2 ^(a 2 tan 2 7 
+ a 2 ) = cot 2 iK* 2 + a 2 ). Thus if ^ 7* ir/2, the locus of the 
gnomonic projections of the points having x zenith ty is the 
hyperbola 

(63) - * = a , y - a. 



The points with ^ = ?r/2 are the points whose projections from 
the line x = 0, ?/ = a, and these are the points with a = 0. 

Finally, the points having fixed x azimuth 7 lie on the horizontal 
line 

z = a tan 7, y = a. 

We have shown how we may determine curves on G that, when 
properly labeled, yield the six spherical coordinates of all points 
in space. 

Gnomonic projections may be made on other planes than the 
one we have selected, but are always made on planes perpendicu- 
lar to a radius of a sphere with center at the origin. The pro- 
jections on planes parallel to the z axis are all called equatorial 
projections, and those on planes perpendicular to the z axis are 
called polar projections. 

EXERCISES 



1. Derive the curves of gnomonic projection for points having a 
constant spherical coordinate a, c, /3, 0, 7, or ty on the following planes 
G, where a > 0: 

(a) y = -a (c) x = -a 

(b) x = a (d) z a 

2. Derive the equations of the curve of gnomonic projection of all 
points of elevation e on the plane perpendicular to the line x = y, z 0. 
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9. Gnomonic charts. A gnomonic chart is a map of a portion 
of a plane that is tangent to a sphere of radius a inches. The 
value of the constant a depends on the reading accuracy desired. 

We take the origin of a rectangular coordinate system on the 
chart to be the point C = (0, a, 0) and the chart to be the plane 
y = a. Then C may also be taken to be the origin of a translated 
(x', y'j z') coordinate system in which the tangent plane at C is 
the x', z' plane, i.e., the plane y r = 0. The x 1 axis is then the 
gnomonic projection of all points having elevation zero and is 
represented on the chart by a horizontal line. The z' axis is the 
gnomonic projection of all points having azimuth zero and is 1 
represented on the chart by a vertical line. All points having 
azimuth a A project into a vertical line and thus the azimuth 
of the projection P f on the chart of a point P in space is the same 
as the azimuth of P and may be read by the use of a scale on the 
x' axis in which the point marked A is at the distance x' = 
a tan A from C. The reading of the coordinate A may be facili- 
tated by the printing on the chart of a vertical grid. It is neces- 
sary for the finite portion of the tangent plane represented by 
the chart to include the portion where k ^ a g A-, where k 
is somewhat more than 45. Then the chart extends to a point 
where x' is greater than a inches. Since very large charts are 
difficult to handle, there are definite physical limitations on the 
use of gnomonic charts as accurate devices for the reading of 
spherical coordinates. 

All points having a fixed value of the coordinate 7 project into 
a horizontal line and therefore may be read by the use of a scale 
on the z' axis in which the point marked E is at a distance of 
z' = a tan E from C. The spherical coordinate /? could be read 
by the use of a circular protractor scale oh the chart and a rotating 
arm pivoted at C to read the coordinate ft = J5. The arm could 
contain a scale like the A scale or it could have a slider that would 
be used to locate the point P' accurately. The arm could then 
be rotated to either the A scale or the E scale for the reading of 
the coordinate 6. 

All points having the same elevation c lie on a cone intersecting 
the tangent plane in a hyperbola. A gnomonic chart then con- 
tains a quasihorizontal family of hyperbolas. The hyperbola 
that cuts the z' axis in the point labeled E has the rectangular 
coordinate equation z = tan E \/x 2 + a 2 , y = a, and all points 
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whose gnomonic projections lie on this hyperbola have elevation 
= E. The number of such curves which are actually drawn 
depends partly on the reading accuracy desired and partly on 
the physical limitations which the drawing of accurate curves 
and readability will impose. The chart must also contain a 
family of quasivertical hyperbolas whose equations are x = 
tan A -\/z 2 + a 2 , y = a, and all the points whose projections 
lie on the hyperbola cutting the x 1 axis in the point labeled A 
have the same x zenith \[/ = 90 A. 

We have now seen how a gnomonic chart may be used as a 
direct reading device for the six spherical coordinates a, |8, ^, 
e, 0, 7 of any point on the chart. Thus, if any two of these 
coordinates are given for any point P in space, a corresponding 
chart projection P' may be plotted and the remaining coordinates 
of P may be read. The chart may be used also as a device for 
computing the effect on a, e of a rotation of axes when the rota- 
tion has been expressed as a product of planar rotations. For the 
rules given in Sec. 5 may be applied and the operations of that 
section become chart motions. 

In using gnomonic charts, the limitation of a chart to a finite 
portion of the tangent plane requires that the charts be inter- 
preted also as tangent planes whose point of tarigency is one of 
the points (a, 0, 0), (-a, 0, 0), (0, -a, 0), (0, 0, a), (0,0, -a) 
as well as the point (0, a, 0) used in the description above. 
The symmetry of our definitions of the three pairs of spherical 
coordinates implies that the six spherical coordinates may be 
obtained by making the same readings as before. However, the 
correspondence between readings and coordinates will be dif- 
ferent for each case. 

The interested student should construct a chart with a = 10 
inches and curves drawn for every three degrees of A and E so as 
to get a better idea of the use of such a chart. 

EXERCISES 

1. Derive the interpretations of the six angles measured on a gnomonic 
chart for each of the positions of C. 

2. Interpret planar rotations as chart motions in the case where 
C = (0, a, 0). Give the interpretations also for C = (0, a, 0), 
C = (a, 0, 0), and C = (0, 0, a). 



CHAPTER 9 
ELEMENTS OF PROJECTIVE GEOMETRY 

1. Homogeneous coordinates. A point P of real n-dimen- 
sional Euclidean geometry was represented in Chap. 1 by a real 
n-dimensional vector, and we are accustomed to writing P = 
(xi, . . . , Xn). The numbers Xi, . . . , x n are unique if a 
coordinate system is specified, and we shall call them the non- 
homogeneous coordinates of P. 

It is sometimes convenient to represent P by an (n + 1)- 
dimensional vector 

P = (j/i, . . . , J/n+l), 

where y n +i is any real nonzero number and y\ = #;2/y+i. We then 
call the vector P a set of homogeneous coordinates of P. Evidently 
Q = (zi 9 . . . , Zn+i) is also a set of homogeneous coordinates of 
the same point P if and, only if Q = tP where t 5^ 0. It follows 
that P = Q if and only if P and Q are linearly dependent. 

We are thus led to the study of the geometry of points P 
represented by corresponding nonzero (n + l)-dimensional vec- 
tors P such that two points P and Q coincide if and only if P and 
Q are linearly dependent. This is a first postulate of the subject 
called projective geometry. The restriction y n +i 7* will be 
omitted, and we only assume that all (n + 1) -dimensional vec- 
tors considered are nonzero vectors. The points corresponding 
to vectors with y n +\ = are usually interpreted as points at 
infinity. 

In the case n = 3, we study all nonzero vectors (x, y, z, t). 
Then the nonhomogeneous coordinates of the points correspond- 
ing to those vectors with t ^ are the ratios 

x y z 
T T T 

2. Lines and planes. Let Pi, . . . , P m be any m distinct 
points whose nonhomogeneous coordinates are given by 

(1) Pi = (Xiy, . . . , Xnj) (j = 1, . , m). 

143 
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Then the vectors 

(2) Pi = (*i/, . . . , x nh 1) (j - 1, . . . , n) 

are m corresponding sets of homogeneous coordinates. Define 
P = (#1, . . . , rr n ) by a set of parametric equations 

(3) P = Pi + 2 (P 2 - Pi) + + UPm - Pi) 

formI independent real parameters 2 , . . . , m. The case 
m = 2 is the case of lines, and the case m = 3 is the case of planes. 
Then if P = (xi, . . . , x n , 1) and Xi, . . . , X m are defined by 



we see that 

(5) P = Xi/ 5 ! + + \mPrn. 

For Pj PI = (x\j x\i, . . . , Xnj Xni, 0) and formula (5) 
is equivalent to 

(6) P - P l = 6(^2 - Pi) + + UPm - Pi). 

The (n + l)st coordinates of the two sides of this vector equation 
are both zero, and the statement of the equality of the first n 
coordinates is precisely formula (3). 

Conversely, let P be defined by formula (5) where the param- 
eters Xi, . . . , Xj range independently over all real numbers. If 

* = Xi + - + \ n 
is zero, P represents a point at infinity. Define 

(7) j = * (j = 2, . . . , m) 
for all finite points, and see that 

(8) Q = (xi, . . . , x m , 1) = P = y 1 P! + 2 P 2 



Then tr l (\i - = t~ l (** + ' ' ' + X m ) = -(2 + + 
and 



(9) Q ~ Pl = Pl 

= {,(A - Pl) + ' ' ' + UPm - Pl). 
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Hence, formula (3) holds, with p = (x\, . . . , x n ). Note that 
if Ri = tJPi, then formula (5) is equivalent to P = n\R\ + 
+ Mm&n, where mRi = XiP t = /xrfJP,-. Then X; = ju^ t - with , ^ 0, 
and the parameters m range over all real numbers when the 
parameters X, do. Hence, formula (3) is equivalent to formula 
(5) where Pi is any vector of homogeneous coordinates of P t . 

We have now shown that the geometric configuration defined 
by the parametric equations of formula (3) is also defined by the 
equations of formula (5). However, formula (5) is homogeneous 
in its parameters and is a much more convenient form for the 
parametric equations of what are called linear spaces in geometry. 

In case m = 2 and PI and P 2 are distinct points, the correspond- 
ing vectors Pi and P 2 are linearly independent. Then the line 
through PI and P 2 is the set of points (including points at infinity) 
whose homogeneous coordinates are all linear combinations of 
Pi and P 2 . In the case n = 3, the homogeneous coordinate 
equations of a line become 

(10) (x, y, z, t) = a(xi, yi, Zi, ti) + 0(x 2 , y 2 , z 2 , t 2 ) 

for independent parameters a and 0. 

Three noncollinear points PI, P 2 , P 3 determine a plane. Then 
Pi and P 2 must be linearly independent and P 3 must not be a 
linear combination of Pi and P 2 . It follows that PI, P 2 , P 3 
determine a plane if and only if Pi, P 2 , P 3 are a set of three linearly 
independent vectors. Then the equations 

P = aPi + /3P 2 + 7 Pa 

are a set of parametric equations of the plane. In the case 
n = 3, these equations become 



(x, y, z, = a(xi, yi, 21, ti) + P(x 2 , y 2 , z 2 , * 2 ) 
We shall limit all further study to the case n = 3. 

EXERCISES 

1. Give a set of parametric equations in homogeneous coordinates of 
the lines defined by following pairs of points : 

(a) (1, -1, 1), (3, -3, 2) (e) (1, 2, -3), (-3, -2, 1) 

(6) (1, 2, 3), (0, 2, 3) (/) (4, 1, 6), (2, -1, 4) 

(c) (-2, 1, 1), (3, -2, 2) (g) (1, 0, 2), (2, -1, 1) 

(d) (3, -1, 2), (1, -1, -2) (h) (1, 1, 1), (2, 2, 5) 
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2. Give a set of parametric equations in homogeneous coordinates of 
the planes defined by the following triples of points: 

(a) (1, - 1,1), (3, -3, 2), (2, 1,0) 
(6) (1,2,3), (0,2,3), (-1,1,2) 

(c) (-2, 1,1), (3, -2,2), (0,1,0) 

(d) (3, -1, -2), (1, -1, -5), (-2, -6,3) 

3. Projective transformations. Let C be any nonsingular 
4X4 matrix 



(11) C = 



and define 

(12) Q = (*', y', z', t'} = PC 

for every point P w r ith corresponding homogeneous coordinate 
vector P = (x, y, z, t). Then every point P is mapped on 
(i.e., determines) a unique point Q, which we shall call the image 
of P under a projective transformation with matrix C. Con- 
versely, if the image Q is given, P is uniquely determined; for 
P = QC~\ 

It should be noted that a projective transformation may map 
finite points on points at infinity. For example, 

/O 1 

(13) (1,0,0,1)1 J 

\1 

It is natural to ask when two matrices B and C define the same 
projective transformation. This occurs when and only when 
Q = PB and R = PC define the same point Q for every P. 
This means that 

(14) PC = t(PB) 

for every P where Ms a nonzero real number which may con- 
ceivably depend on P. We shall prove the following theorem. 

Theorem 1. Two matrices B and C define the same protective 
transformation if and only if C is a nonzero scalar multiple of B. 

For if C = IB, it is evident that PC = t(PB). Conversely, 
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suppose that PC is a nonzero scalar multiple of PB for every 
P. Take P = (1, 0, 0, 0) and see that PB is the first row BI 
of B. Then the first row of C must be a scalar multiple t\B\ of 
the first row of B. We have a similar result for the other rows 
and see that C = DB, where D = diag{ti, t z , t s , I*} is a diagonal 
matrix. We next take P = (1, 1, 1, 1) and see that 

(15) PC = t l B l + t 2 B 2 + tzB 3 + tiB* = tPB 

= t(B l + B 2 + Bs + B,) 

where Bi is the ith row of A . Then 

(16) (ti - t,t 2 - t, t z - t,tt - t)B = 0. 

Since B is nonsingular, t\ = t 2 = 3 = t t, and we have proved 
that C = tB. 

4. Tetrahedral coordinates. If PI is any point, the corre- 
sponding vector Pi is a nonzero vector with the property that all 
nonzero scalar multiples of Pi define the same point PI. Let 
P 2 be a point distinct from PI. Then PI and P 2 are linearly 
independent, and the set of all linear combinations of Pi and P 2 
is the set of all vectors P defining points on the line joining PI 
toP 2 . 

We next let P 3 be a point not on the line joining PI and P 2 . 
Then Pi, P 2 , P 3 are linearly independent, and the set of all linear 
combinations of Pi, P 2 , P 3 is the set of all vectors P defining 
points P on the plane determined by PI, P 2 , PS. 

Suppose, finally, that P 4 is a point not on the plane determined 
by Pi, P 2 , P 3 . Then P 4 is not a linear combination of Pi, P 2 , P 3 . 
But then Pi, P 2 , P 3 , P 4 are linearly independent. It follows that 
the matrix 



(17) 



whose rows are the vectors Pi, P 2 , P 3 , P 4 , is a nonsmgular matrix. 
For |J3| = if and only if there exist numbers Xi, X 2 , X 3 , X 4 not 
all zero such that XiPi + X 2 P 2 + XaPs + X 4 P 4 = 0. 

Let PI, P 2 , P 3 , P 4 be four fixed points which are not coplanar 
and let C be the nonsingular matrix formed by the corresponding 
vectors Pi, P 2 , P 3 , P*. Then every vector P = (x, y, z, t) defines 
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a unique set of numbers (x', y', z', t') such that 
(18) P = z'Pi + y'P 2 + z'P 3 + t'P*. 

We shall call these numbers a set of tetrahedral coordinates of P 
and shall say that Pi, P 2 , PS, P are the vertices of a tetrahedron 
of reference. The existence of such coordinates is due to the 
fact that formula (18) is equivalent to 

(]Q\ (r 11 ? t\ (v* ?/ y 9 /Mft 

\i\j) \x, y, z, i) \x y y , z , i )>, 

and therefore to 

(20) (x 1 , y', z', t) = (x, y, z, <)C, 

where C = B~ l . 

We have been using the points whose homogeneous coordinates 
are (1, 0, 0, 1), (0, 1, 0, 1), (0, 0, 1, 1) and (0, 0, 0, 1) as a tetra- 
hedron of reference and we see that (x, y, z, t} = x(l, 0, 0, 1) + 
y(0, 1, 0, 1) + z(0, 0, 1, 1) + (t - x - y - z)(0, 0, 0, 1). How- 
ever, x, y, z, t are actually the tetrahedral coordinates of a point 
relative to the tetrahedron defined by the vectors (1, 0, 0, 0), 
(0, 1, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1). 

The equations of formula (20) have been interpreted in this 
section as the relations connecting two sets of coordinates of a 
fixed point. They were given in formula (12) of Sec. 3 as defining 
a point-to-point correspondence called a protective transformation. 
It is important to observe that these two interpretations are two 
geometric interpretations of the same algebraic phenomenon. 

We should also note that if L is a three-rowed orthogonal 
matrix and 

fen c-( L 0> ! 

\*tL) C/ I A 1 ]j 

\ !/ 

then the equations 

(22) (x, y, z, 1) = (x', y', z', 1) (^ j) 

define a projective transformation of coordinates with orthogonal 
matrix C. But these equations are clearly equivalent to the 
orthogonal transformation 

(23) (x, y, z) = (a/, y', z')L 
on nonhomogeneous coordinates. 
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5. The unit point. As we have seen, any set of four points 
PI, P 2 , PS, P\ that are not coplanar defines a corresponding 
matrix B of formula (17) and hence a set of tetrahedral coordi- 
nates. However, the same points also define a matrix 



(24) 



for any set of nonzero real numbers fa, . . . , < 4 . It follows that 
PI, . . . , P 4 do not completely specify B. The specifcation 
may be completed however by prescribing a fixed point (XQ, y^ z , 
to) as a point called a unit point and which is such that XQ = y 
= z ' = <o' = 1. For then we have 





(25) (x,y, z, t) = t(l, 1, 1, 1) 

which is equivalent to 

(26) t(ti, /2, 3, tt)B = (x , ?/o, < 
and therefore to 

(27) t(faj fay ts, 4 ) (#0, yo, ZQ, 

But this uniquely determines the numbers (fa, t%, fa, fa) apart from 
a proportionality factor t. 

Note that the unit point is not arbitrary since we must have 
fay 2, fa, <4 all not zero. But the equation t(fa, fa, fa, t*)B = 
(#o, yo, ZQ, fa), with one of ti, t^ fa, t\ zero, implies that (XQ, yo, ZQ, to) 
is a linear combination of three of the vectors Pi, P 2 , P 3 , P 4 and 
therefore PO is a point on a plane through three of the vertices of 
the tetrahedron of reference. The four planes determined by 
selecting three of the four vertices are called the faces of the 
tetrahedron of reference, and we have shown that a point may be 
selected as unit point if and only if it is not a point of a face of the 
tetrahedron. 

EXERCISES 

1. Let the vertices of a tetrahedron of reference be PI = (f, $, f), 
P 2 = (i i -*), P 3 = (-1, I *), P 4 = (- J, i *) and the unit point 
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be (2, 1, 1). Find the tetrahedral coordinates of the points whose 
homogeneous coordinates are given by the following vectors: 

(<*) (1,0,0,1) (e) (1, -1,2,0) 

(b) (1,1,1,1) (/) (4,1, -1,1) 

(c) (1,0,0,0) (jf) (-1,2, -1,1) 

(d) (0, 0, 0, 1) (A) (1, 2, 2, 1) 

2. Let the vectors of Exercise 1 be the tetrahedral coordinates of a set 
of points. Find the (nonhomogeneous) rectangular coordinates of the 
finite points of the set. 

6. Invariant points. A project! ve transformation with matrix 
C maps every point P on a unique image Q where Q = PC. 
Then we call P an invariant point if Q = P, that is, P is its own 
image. But this occurs if and only if PC = tP, that is, 

(28) P(tl - C) = 0. 

Since P must be a nonzero vector, the determinant 

(29) |*7 - C| = 0. 

Thus t must be a characteristic root of the matrix C, and we have 
shown that P is an invariant point if and only if P is a characteristic 
vector of C. 

If { is a simple root of xl C = 0, an argument like that used 
in the proof of Theorem 17 of Chap. 6 may be used to prove that 
there is only a single invariant point. When t is a double root of 
xl C 0, there are two corresponding linearly independent 
vectors and all linear combinations of them are characteristic 
vectors. Then there is a line of invariant points corresponding 
to the root t. The root t may be a triple root of xl C = and 
then there are three linearly independent characteristic vectors 
and a corresponding plane of invariant points. In the final case 
of a root of multiplicity four, C = tl and all points are invariant 
points. 

7. Quadric surfaces. An algebraic surface of degree n is 
defined by a polynomial equation f(x, y, z) = where f(x, y, z) 
is a polynomial of degree n. Then the equation 



(30) ^(x,y,,0= 

is an equation in homogeneous coordinates of the same surface. 
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Moreover, the polynomial F(x, y, z, t) is a homogeneous poly- 
nomial in x, y, z, t and this is the source of the " homogeneous " 
coordinate terminology. 

A projective transformation (x, y, z, t) = (#', y'j z', t')B re- 
places F(x, y, z, t) by another homogeneous polynomial and the 
corresponding equation <(x', y', z', t') = is an equation in 
tetrahedral coordinates of the given surface. In particular a 
quadric surface is defined by a polynomial equation F(x } y, z, t) 
= where 

(011 012 013 014\ /X] 
012 022 023 024 1 I V 
a 
013 023 033 034 
014 024 034 044/ V, 

Thus F(x, y, z, t) = PAP*, where A is a symmetric matrix. 
Then 

(32) 0(z', 7/, z 1 , t'} ^ QBAB*Q* 

and the matrix of this quadratic form is the matrix BAB*. 
It can be shown that the matrix B may be selected so that 

(di 
o o' S, o 
dt} 

where di = 1, 1, or zero. The corresponding quadratic form 
is then d\x' 2 + d^y' 2 + d&' 2 + d\t t<2> . We shall not prove this 
result here but refer the reader to the author's " Introduction to 
Algebraic Theories " for proof. The number of nonzero diagonal 
terms is the rank of the quadratic form F(x, y, z, f) and the 
number of positive diagonal elements di is its index. Both 
integers are invariants of the form. Note that \A$\ = |J5| 2 |A|. 

The determinant of the matrix A of the quadratic form 
F(x, y, z, t) is called a discriminant A of the corresponding quadric 
surface, and we have seen that a projective transformation of 
coordinates replaces A by |5| 2 A, where B is the matrix of the 
transformation. We call the surface nonsingular if A ^ and 
see that then |B| 2 A = did^ld* ^ and each di = 1 or 1. By 
permuting variables and multiplying by 1 if necessary, we see 
that every nonsingular quadric surface is defined with respect 
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to a properly chosen tetrahedral coordinate system by one and 
only one of the equations 

x z + y* + z 2 + t* = 0, 
(34) x 2 + y* + z* - t 2 = 0, 

8. Cross ratios of points. If the rows of the matrix 

jj _ (ai a 2 a n \ 

\bi b 2 ' b n / 

are not proportional, then it must be possible to select two col- 
umns of U such that the rows of 



(35) 



(en a 
\bi b 



are not proportional. Then these rows are linearly independent 
and the only solution of 






' = xi(a i} a y ) + 2 (6,-, &,-) = 



ii a/ 
>,- 6y 



is Xi = #2 = 0. It follows that the determinant 

(36) A = 
Let us now consider a matrix 

(37) A = 




whose rows are the homogeneous coordinates of four distinct 
collinear points. Then we have the following result: 

Lemma 1. The two-rowed determinants defined by two columns 
of A are all not zero or all zero. 

For let Ui = (ai, 61, Ci, di) and U 2 = (a 2 , b 2) c 2 , d 2 ) be the 
coordinates of any two distinct points on the line joining the four 
given points and so have 

/oo\ / j\ ___ \TJ I .77 
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A = 




di #2 0,3 (14 
62 63 64 



) 



and the two-rowed matrices in the jth and Mh columns of A are 
the two-rowed square submatrices of 



(40) 



f\i H\ 

X2 M2 
Xs M3 

^4 M4/ 

If we let D P(/ be the determinant formed by the pth and qth rows 
of this matrix, we have 




(41) 



x,, 



If (X p , Up) and (X <; , M) are proportional, then so are (x p , y p , z p , t P ) 
and (x q , y q , z q , t q ) contrary to the hypothesis that the four given 
points are distinct. It follows that the determinants 



(42) 



are all not zero and therefore D pg = if and only if the deter- 
minant 



(43) 



a, a k 



which is independent of p and g, is zero. This proves the lemma. 
The cross ratio of four distinct collinear points PI, P 2 , Pa, P* is 
the ratio 



(44) 



of the corresponding two-rowed minors D Pq obtained from a 
corresponding matrix 
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of homogeneous coordinates by selecting two columns of A in 
which the minors D pg are not zero. Then if Ui and C/ 2 are the 
vectors of homogeneous coordinates of any two points on the 
given line through Pi, P 2 , Pa, P 4 , we see that 



(45) 



k(Pi, P 2 , Pa, PO = 



AL 

Xi 
X 4 



X2 M2 



X2 

X 3 M3 



Then k(Pi, P 2 , PS, P 4 ) is independent of the particular columns of 
A selected in its computation. It is independent of the base 
vectors f/i and C/ 2 , since it is not defined in terms of these vectors 
but rather in terms of homogeneous coordinates of the given 
points. Finally, it should be noted that the replacement of 
Pi, P 2 , Ps, P 4 by *iPi, tj*i, ^sPa, * 4 P 4 will not alter k(Pi, P 2 , P 3 , P 4 ). 
Theorem 2. The cross ratio of four points is unaltered by a 
projective transformation. 



For let 
have 



(46) 



= PiC, Q 2 = P 2 C, <J 3 = P 3 C, Q 4 = P 4 C. Then we 




= AC = 




Since /r(Pi, P 2 , P 3 , P 4 ) is independent of the selection of base 
points, the cross ratio obtained with the use of the vectors 
UiCy U%C will have the same value as that obtained by the use of 
Ui, C/ 2 . Then k(Qi, Q 2 , Qs, QO = *(Pi, P 2 , P 3 , P 4 ) as desired. 

The converse of Theorem 2 is true, and we shall prove it. We 
state the result as the following 

Theorem 3. Let PI, P 2 , P 3 , P 4 and Qi, Q 2 , Q 3 , Q^ be two sets 
each of four distinct collinear points, and let 

k = Ar(Pi, P 2 , P 8 , P 4 ) = fr(Qi, Q 2 , Qs, <20. 

T/ien ffcere exists a projective transformation carrying Qi into Pi, 
/or i = 1, 2, 3, 4. 

To prove this result, we first notice that PI and P 2 may be 
used as base points and therefore P 3 = tj*i + 2 P 2 , where neces- 
sarily <i< 2 5^ 0. Then the vector tjP* is a homogeneous coordinate 
representation of P t and may be taken to be P>, for i =*= 1, 2. 
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Hence, we may take P 3 = PI + -Pa- Now P^ = siPi + szPa, 
and we see as before that Si2 7* 0. Then we may replace P 4 by 
S2 -1 P^ and thus obtain P 4 = kPi + P%, where A 1 is a nonzero real 
number. This yields the relation 



(47) 



A = 




and we use formula (45) to see that 

1 



(48) 



k(Pi, P t , P., P.) = 



1 1 



= k. 



By the argument of Sec. 4, there exists a nonsingular matrix G 
whose first two rows are Pi and Pa and then 



(49) 



A = 




may be chosen so 




In a similar fashion the vectors Qi, 
that the corresponding matrix 



B = 



where H is a nonsingular matrix whose first two rows are Q\ and 
Q 2 . The value of k is the same and B = AC, where C = G~ 1 H 
is a matrix defining the desired projective transformation. This 
proves the theorem. 

The form of the matrix A given by formula (47) may be used 
to determine the effect of permuting the points Pi, PI, P&, Pt. 
There are 4! = 24 possible permutations and thus possibly 24 
distinct values. However, there, are actually only six formally 
distinct values. These are the values 



(50) k, 



1 

F 



1 



1 



1 



k- 1 



f56 



SOLID ANALYTIC GEOMETRY 



[Chap. 9 



obtained from k by the two operations of inverting k and sub- 
tracting k from 1. That these are values follows from the fact 
that 



(51) 
and 
(52) 
so that 

(53) 



Pi, Ps, 





1 




1 






1 1 




k 1 


1 




1 




1 


*' 




k 1 




1 1 




1 




1 1 




1 

T 6 





k 1 


= 1 -*, 


1 1 


k 1 




1 





, PI, p 2 , p 4 ) = 



*(P3, P s , Pi, 



k - 1 



But from each of these six values we obtain three other equal 
values by permutation, since it should be evident from formula 
(45) that 



(54) 



8 , P,, P 4 ) = k(P t , Pi, P 4 , P,) = k(P,, P 4 , P 1( P 2 ) 



ILLUSTRATIVE EXAMPLE 

Compute the cross ratio of PI = (1, 1, 2), P 2 = (-3, 1, 1), P 3 = 
(-T, I,T), ^4 = (-i 1,1) and the cross ratio of Q, = (-|,T,|), 
Q 2 = (-1, 1, 0), Q 3 = (0, 2, 1), Q 4 = (-i, |, |), and show that the 
first four points can be carried into a proper arrangement of the second 
set of four points by a projective transformation. 



Solution 



We form the matrices 



A = 




B = 
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Then we use the first and second columns to compute 



*(Pi, P,, P, f P) = 



1 1 

-7 5 




-3 1 

-1 3 


1 1 
-1 3 




-3 1 

-7 5 



and 



-3 711-1 1 




2 

Then k(Q\, Q-i, Q 3 Q 4 ) = 3, and there exists a projective transformation 
that will carry Pi into Q 3 , P 2 into Q 2 , Pa into Qi and P 4 into Q 4 . 

EXERCISES 

1. Compute k(Pi, P 2 , P 3 , P 4 ) and k(Qi, Q 2 , Q 3 , <M in the following 
cases, and find an arrangement of the second set of points into which the 
first set can be carried by a projective transformation. 

(a) P! = (2, 1, -3), P 2 = (-1, 2, 0), P 3 = (-1, -5, -3), P 4 = 

(4, -3, -3) 
Q l = (1, -2, -1), Q 2 = (-1, 0, 1), Q 3 = (-1, 1, 1), Q 4 = (-4, 

3,4) 

(6) P 1 = (1, 1,0), P 2 = (-1,1, 1),P 3 = (-1,3,2),P 4 = (2,0, -1) 
Q! = (2, -1, 1), Q, = (1, 1, -1), Q 3 = (5, -1, 1), Q 4 = (-1, 
2 2) 

(c) P! = (1, -1, 1, 1), P 2 = (1, 0, -1, 2), P 3 = (0, -1, 2, -1), P 4 

= (5, -1, -3,_9) 

Q, = (2, 1, -1, 0), Q 2 = (6, 5, 3, 4), Q 3 = (1, 1, 1, 1), Q 4 = (1, 
0, -2, -1) _ 

(d) P! = (0, 1, 0, 1), P 2 = (1, 1, 1, 1), P, = (2, 5, 2, 5), P 4 = (1, 3, 

Q, = (1, 1, 0, 0), Q 2 = (2, 3, 1, 0), Q 3 = (0, 1, 1, 0), Q 4 = (7, 

(e) P! = (0, 1, 0, 4), P 2 = (2, 3, -2, 2), P, = (1, 1, -1, -1), P 4 

= (1,2, -1,3) 
Q, = (0, 1, 0, 4), Q 2 = (2, 3, -2, 2), Q 3 = (2, 4, -2, 6), Q 4 = 

(2, 2, 2, 2) 
(/) P! = (2, 0, 0, 1), 'P 2 = (1, 2, 3, 4), P 3 = (-1, 2, 3, 3), P 4 = (5, 

2,3,6) 
Q, = (2, 0, 0, 1), Q 2 = (5, 2, 3, 6), Q 3 = (-3, 2, 3, 5), Q 4 = (9, 

2, 3, 8) 

2. Use the method of proof of Theorem 3 to show that a projective 
transformation can be found which will carry any five points PI, . . . , 
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PB, no four of which are coplanar, into any other five noncoplanar points. 
HINT: Show that P lt . . . , P 5 can be carried into (1, 0, 0, 0), (0, 1, 0, 0), 
(0,0,1,0), (0,0,0,1), (1,1,1,1). 

9. Plane coordinates and duality. The coefficients of the 
homogeneous coordinate equation ax + by + cz + dt = form 
a nonzero vector (a, fr, c, d) that may be regarded as a set of 
coordinates of the corresponding plane. Clearly two nonzero 
vectors define the same plane if and only if each is a nonzero real 
scalar multiple of the other, and thus we have a new coordinate 
system for our projective geometry. We may then express all 
our projective geometric properties in terms of plane coordinate 
vectors, rather than point coordinate vectors. The reader 
should carry this out for the material of Sec. 2. 

A point, represented in point coordinates by (x, y, z, t), is said 
to be incident with a plane, represented in plane coordinates by 
(a, 6, c, d), when the point lies on the plane (the plane contains 
the point). This occurs if and only if the inner product of the 
two corresponding vectors vanishes, i.e., if and only if 

(x, y, z, t) - (a, 6, c, d) = ax + by + cz + dt = 0. 

Since the incidence relation is unaltered by the interchange of the 
point and plane vectors, we have a property of our geometry that 
we may state as the following principle: 

DUALITY PRINCIPLE. // the words point and plane are inter- 
changed in any true statement about the incidence of points and 
planes, the result is a true statement. 

We may also formulate this property by saying that if a geo- 
metric configuration formed out of points and planes has certain 
incidences the configuration obtained by the interchange of 
points and planes will have the same incidences. 

The Duality Principle may be extended to linear subspaces of 
an n-dimensional space. The precise results may be found in 
more advanced texts on projective geometry. 

EXERCISES 

1. State the dual of the theorem that three points not on a line deter- 
mine a plane. 

2. Find some other incidence relation, and state its dual. 
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